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Abstract 

These two accompanying papers are concerned with entanglement for sys¬ 
tems of identical massive bosons and the relationship to spin squeezing and 
other quantum correlation effects. The main focus is on two mode entangle¬ 
ment, but multi-mode entanglement is also considered. The bosons may be 
atoms or molecules as in cold quantum gases. The previous paper I dealt with 
the general features of quantum entanglement and its specific definition in the 
case of systems of identical bosons. Entanglement is a property shared between 
two (or more) quantum sub-systems. In defining entanglement for systems of 
identical massive particles, it was concluded that the single particle states or 
modes are the most appropriate choice for sub-systems that are distinguishable, 
that the general quantum states must comply both with the symmetrisation 
principle and the super-selection rules (SSR) that forbid quantum superposi¬ 
tions of states with differing total particle number (global SSR compliance). 
Further, it was concluded that (in the separable states) quantum superpositions 
of sub-system states with differing sub-system particle number (local SSR com¬ 
pliance) also do not occur. The present paper II determines possible tests for 
entanglement based on the treatment of entanglement set out in paper I. 

Several inequalities involving variances and mean values of operators have 
been previously proposed as tests for entanglement between two sub-systems. 
These inequalities generally involve mode annihilation and creation operators 
and include the inequalities that define spin squeezing. In this paper, spin 
squeezing criteria for two mode systems are examined, and spin squeezing is 
also considered for principle spin operator components where the covariance 
matrix is diagonal. The proof, which is based on our SSR compliant approach 
shows that the presence of spin squeezing in any one of the spin components 
requires entanglement of the relevant pair of modes. A simple Bloch vector test 
for entanglement is also derived. Thus we show that spin squeezing bdecomes 
a rigorous test for entanglement in a system of massive bosons, when viewed as 
a test for entanglement between two modes. 

In addition, other previously proposed tests for entanglement involving spin 
operators are considered, including those based on the sum of the variances for 
two spin components. All of the tests are still valid when the present concept 
of entanglement based on the symmetrisation and super-selection rule criteria 
is applied. These tests also apply in cases of multi-mode entanglement, though 
with restrictions in the case of sub-systems each consisting of pairs of modes. 
Tests involving quantum correlation functions are also considered and for global 
SSR compliant states these are shown to be equivalent to tests involving spin 
operators. A new weak correlation test is derived for entanglement based on 
local SSR compliance for separable states, complementing the stronger corre¬ 
lation test obtained previously when this is ignored. The Bloch vector test is 
equivalent to one case of this weak correlation test. Quadrature squeezing for 
single modes is also examined but not found to yield a useful entanglement test, 
whereas two mode quadrature squeezing proves to be a valid entanglement test, 
though not as useful as the Bloch vector test. The various entanglement tests 
are considered for well-known entangled states, such as binomial states, rela- 
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tive phase eigenstates and NOON states - sometimes the new tests are satisfied 
whilst than those obtained in other papers are not. 

The present paper II then outlines the theory for a simple two mode inter¬ 
ferometer showing that such an interferometer can be used to measure the mean 
values and covariance matrix for the spin operators involved in entanglement 
tests for the two mode bosonic system. The treatment is also generalised to 
cover multi-mode interferometry. The interferometer involves a pulsed classical 
field characterised by a phase variable and an area variable defined by the time 
integral of the field amplitude, and leads to a coupling between the two modes. 
For simplicity the centre frequency was chosen to be resonant with the inter¬ 
mode transition frequency. Measuring the mean and variance of the population 
difference between the two modes for the output state of the interferometer for 
various choices of interferometer variables is shown to enable the mean values 
and covariance matrix for the spin operators for the input quantum state of the 
two mode system to be determined. The paper concludes with a discussion of 
several key experimental papers on spin squeezing. 


PACS Numbers 03.65 Ud, 03.67 Bg, 03.67 Mn, 03.75 Gg 
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1 Introduction 


The previous paper I dealt with the general features of quantum entanglement 
and its specific definition in the case of systems of identical bosons. In dehn- 
ing entanglement for systems of identical massive particles, it was concluded 
that the single particle states or modes are the most appropriate choice for 
sub-systems that are distinguishable. Further, it was conclude that the general 
quantum states must comply both with the symmetrisation principle and the 
super-selection rules (SSR) forbidding quantum superpositions of states with 
differing total particle number (global SSR compliance). As a consequence, it 
was then reasoned that in the separable states quantum superpositions of sub¬ 
system states with differing sub-system particle number (local SSR compliance) 
do not occur [T] . Other approaches - such as sub-systems consisting of labelled 
indistinguishable particles and entanglement due to symmetrisation [2] or al¬ 
lowing for non-entangled separable but non-local states [3]- were found to be 
unsuitable. The local (and global) SSR compliant definition of entanglement 
used here was justified on the basis of there being no non-relativistic quantum 
processes available to create SSR non-compliant states and alternatively on the 
absence of a phase reference [3]. 

Paper I can be summarised as follows. Section 2 covered the key definitions 
of entangled states, the relationship to hidden variable theory and some of the 
key paradoxes associated with quantum entanglement such as EPR and Bell 
inequalities. A detailed discussion on why the symmetrisation principle and the 
super-selection rule is invoked for entanglement in identical particle systems was 
discussed in Section 3. Challenges to the necessity of the super-selection rule 
were outlined, with arguments against such challenges dealt with in Appendices 
D and E. Two key mathematical inequalities were derived in Appendix B and 
details about the spin EPR paradox set out in Appendix C. The final Section 4 
summarised the key features of quantum entanglement discussed in the paper. 

The present paper II focuses on tests for entanglement in two mode systems 
of identical bosons, with particular emphasis on spin squeezing and correlation 
tests and how the quantities involved in these tests can be measured via two 
mode interferometry. Two mode bosonic systems are of particular interest be¬ 
cause cold atomic gases cooled well below the Bose-Einstein condensation (BEC) 
transition temperature can be prepared where essentially only two modes are 
occupied ( 0 , 0 ). This can be achieved for cases involving a single hyperfine 
components using a double well trap potential or for two hyperfine components 
using a single well. At higher temperatures more than two modes may be occu¬ 
pied, so multi-mode systems are also of importance and the two mode treatment 
is extended to this situation. 

As well as their relevance for entanglement tests, states that are spin squeezed 
have important applications in quantum metrology. That squeezed states can 
improve interferometry via the quantum noise in quadrature variables being 
reduced to below the standard quantum limit has been known since the pio- 
neeering work of Caves [7] on optical systems. The extension to spin squeezing 
in systems of massive bosons originates with the work of Kitagawa and Ueda 
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[5], who considered systems of two state atoms. As this review is focused on 
spin squeezing as an entanglement test rather than the use of spin squeezing 
in quantum metrology, the latter subject will not be covered here. In quantum 


metrology involving spin operators the quantity ^ I (which in¬ 

volves the variance and mean value of orthogonal spin operators) is a measure of 
the uncertainty A9 in measuring the interferometer phase. The interest in spin 
squeezing lies in the possibility of im prov ement over the standard quantum limit 
where A9 = \/\fN (see Subsection 3.7). As we will see, for squeezed states 


ASl) < 


Sz) 1/2 so we could have A9 l/\/l 


^ 1/ViV 

which is less than the standard quantum limit. In SubSection |3.8| we give an 
example of a highly squeezed state where A9 ~ Vln A /N which is near the 
Heisenberg limit. Suffice to say that increasing the number of particles in the 
squeezed state has the effect of improving the sensitivity of the interferometer. 
Aspects of quantum metrology are covered in a number of papers (see [5], [TTi]'). 
based on concepts such as quantum Fisher information, Cramers-Rao bound 
m, [12], quantum phase eigenstates. 

The proof of the key conclusion that spin squeezing in any spin component 
is a sufficiency test for entanglement [T| is set out in this paper, as is that for 
a new Bloch vector test. A previous proof |13j that spin squeezing in the z 


spin component W (AS 


Sx) I/2demonstrates entanglement based on 


< 

treating identical bosonic atoms as distinguishable sub-systems has therefore 
now been superceded. It is seen that correlation tests for entanglement of quan¬ 
tum states complying with the global particle number SSR can be expressed in 
terms of inequalities involving powers of spin operators. Section sets out the 
definitions of spin squeezing and in the following Section]^ it is shown that spin 
squeezing is a signature of entanglement, both for the original spin operators 
with entanglement of the original modes, for the principle spin operators with 
entanglement of the two new modes and finally for several multi-mode cases. 
Details of the latter are set out in Appendices El and A number of other 
correlation, spin operator and quadrature operator tests for entanglement pro¬ 
posed by other authors are considered in Sections EE and with details of 
these treatments set out in Appendices [MjandflS] Some tests also apply in cases 
of multi-mode entanglement, though with restrictions in the case of sub-systems 
each consisting of pairs of modes. A new weak correlation test is derived and 
for one case is equivalent to the Bloch vector test. 

In Section it is shown that a simple two mode interferometer can be used 
to measure the mean values and covariance matrix for the spin operators in¬ 
volved in entanglement tests, with details covered in Appendices El and [TS] 
The treatment is also generalised to cover multi-mode interferometry. Actual 
experiments aimed at detecting entanglement via spin squeezing tests are ex¬ 
amined in Section The final Section summarises and discusses the key 
results regarding entanglement tests. Appendices and IR provide details re¬ 
garding certain important states whose features are discussed in this paper - the 











’’separable but non-local ” states and the relative phase eigenstate. 
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2 Spin Squeezing 

The basic concept of spin squeezing was first introduced by Kitagawa and Ueda 
[8] for general spin systems. These include cases based on two mode systems, 
such as may occur both for optical fields and for Bose-Einstein condensates. 
Though focused on systems of massive identical bosons, the treatment in this 
paper also applies to photons though details will differ. 

2.1 Spin Operators, Bloch Vector and Covariance Matrix 

2.1.1 Spin Operators 

For two mode systems with mode annihilation operators d, b associated with 
the two single particle states \4>a)^ \4>b)^ where the non-zero bosonic com¬ 
mutation rules are [e, e^] = 1 (e = d or 5), Schwinger spin angular momentum 
operators = x, y, z) are defined as 

Sx = (b^a + a^b)/2 Sy = (b^a — a^b)/2i = (b^b — a^a)/2 (1) 

and which satisfy the commutation rules [5"^ , 5"^ ] = for angular mo¬ 

mentum operators. For bosons the square of the angular momentum operators 
is given by S‘^ + Sy+ 3'^ = {N/2){N/2 + 1), where N = {b^b + a^a) is the boson 
total number operator, those for the separate modes being dg = e^e (e = d 
or b). The Schwinger spin operators are the second quantization form of sym¬ 
metrized one body operators S;^ = (</>a(OI + l</'a(*)) (</>b(*)l)/2 ; % = 

(</'b(*)|-|0a(*)) (</'a(*)l)/2 

, where the sum i is over the identical bosonic particles. In the case of the two 
mode EM field the spin angular momentum operators are related to the Stokes 
parameters. 

2.1.2 Bloch Vector and Covariance Matrix 

If the density operator for the overall system is 'p then expectation values of the 
three spin operators = Tr{pS^) = x,y,z) define the Bloch vector. Spin 

squeezing is related to the fluctuation operators AS^ = in terms of 

which a real, symmetric covariance matrix C{S^,Sy) p = x,y,z) is defined 
[H], H] via 

C(S^,Sy) = {(^AS^ASy^^(^ASyAS^))/2 

= (s^Sy + S,^l2-(s^)(Sy) ( 2 ) 

and whose diagonal elements C{S^,S^) = gives the variance for the 

fluctuation operators. The covariance matrix is also positive definite. The 
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variances for the spin operators satisfy the three Heisenberg uncertainty prin¬ 
ciple reations (^ASy'^'j > 3 I <^5^^ p; (^ASy'^'^ (^AS^'^'j > || P; 

(^ASz"^^ (^ASx'^^ > j\ spin squeezing is defined via conditions such 

as (^ASx^'j < || I with (^ASy"^^ > (^Sz^ |, for Sx being squeezed com¬ 

pared to Sy and so on. Spin squeezing in these components is relevant to tests 
for entanglement of the modes a and b, as will be shown later. Spin squeezing 
in rotated components is also important, in particular in the so-called principal 
components for which the covariance matrix is diagonal. 

2.2 Spin Operators - Multi-Mode Case 

As well as spin operators for the simple case of two modes we can also de¬ 
fine spin operators in multimode cases involving two sub-systems A and B. 
For example, there may be two types of bosonic particle involved, each compo¬ 
nent distinguished from the other by having different hyperfine internal states 
|A), \B). Each component may be associated with a complete orthonormal set 
of spatial mode functions and so there will be two sets of modes 

\4>ai), I'/'bi)) where in the |r) representation we have (r = (/'ai(i’) 1 ^) 

(r \(j)yf) = \B). Mode orthogonality between A and B modes arises from 

{A\B) = 0 rather from the spatial mode functions being orthogonal. 

We can define spin operators for the combined multimode A and B sub¬ 
systems m via 

Sx = ^/dr ($];(r)$a(r)-b $j^(r)${,(r)) 

Sy = ^/dr (§|;(r)$a(r) - $l(r)$b(r)) 

= i/dr ($t(r)$,(r)-$t(r)$^(r)) (3) 

where the field operators satisfy the non-zero commutation rules 

='^cdd(r-r') c,d = a,b (4) 

It is then easy to show that the standard spin angular momentum conmmutation 
rules are satisfied. [5"^ , ] = ie^y\S\. 

For convenience we can expand the field operators in terms of an orthonormal 
set of spatial mode functions 4>i{r). We can choose the spatial mode functions to 
be the same 4 >ai{^) ~ —4‘ii'^) (these might be momentum eigenfunctions) 

and then the field annihilation operators for each component are 

^a(r) = Eoi ^bir) = J2bi ^,(r) (5) 

i i 

These expansions are consistent with the field operator commutation rules 0 
based on the usual non-zero mode operator commutation rules [a^, a]] = [bi, 6 |] = 

Sij. 
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By substituting for the field operators we can then express the spin operators 
in terms of mode operators as 



and it is then easy to confirm that the standard spin angular momentum con- 
mmutation rules are satisfied. [5'.t , 5'^ ] = ie^^xSx- We now have both field 
and mode expressions for spin operators in multimode cases involving two sub¬ 
systems A and B. 

Finally, the total number of particles is given by 

N = /dr ($j;(r)$b(r)-h$l(r)$a(r)) 

= H(bfbi + alai^ 

= m + Na (7) 

in an obvious notation. 

2.3 New Spin Operators and Principal Spin Fluctuations 

The covariance matrix has real, non-negative eigenvalues and can be diago- 
nalised via an orthogonal rotation matrix M(—a, —/3, — 7 ) that defines new spin 
angular momentum operators = x, y, z) via 

= ( 8 ) 

and where 

C(Jj,j;) = Y,M^x{-a,-l3,-^)C{Sx,Se)M^e{-a,-l3,-^) 
xe 

= (9) 

is the covariance matrix for the new spin angular momentum operators = 
X, y, z), and which is diagonal with the diagonal elements and 

giving the so-called principal spin fluctuations. The matrix M(a,/ 3 , 7 ) 
is parameterised in terms of three Euler angles a, j3, 7 and is given in |lfij (see 
Eq. (4.43)). 

The Bloch vector and spin fluctuations are illustrated in Figure 1. In Fig 1 
the Bloch vector and spin fluctuation ellipsoid is shown in terms of the original 
spin operators {f = x,y,z) 
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Figure 1. Bloch vector and spin fluctuations shown for original spin operators. 


These rules also apply to multimode spin operators as defined in SubSection 

o 
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2.4 Spin Squeezing Definitions 

We will begin by considering the case of the spin operators in the most general 
case. We will also specifically consider spin squeezing for the two new modes 
and for the multi-mode situation. 


2.4.1 Heisenberg Uncertainty Principle and Spin Squeezing 


Since the spin operators also satisfy Heisenberg uncertainty principle relation¬ 
ships 


A5,2) 

(a5/) 

> 

4I 



A5/) 

(a5,2^ 

> 

4I 

(5.) P 


ASz^) 

(a^^v 

> 

aI 


(10) 


spin squeezing will now 

be defined for the spin operators 

via condtions such as 

(aV) 

< 


(a>i 

and ^A5j,^y 

'>51 

(a>i 


(ASy^) 

< 

2I 

( 5 .) 

and 


( 5 .) 1 


( a 5 /) 

< 

2I 

{Sy)\ 

and (^ASx^^^ 


1 (Sy) 1 

(11) 

for 53;being squeezed compared to Sy, 

and so on. 





Note also that the Heisenberg uncertainty principle proof (based on 




-I- iXASisj yASa — iXASp j ^ > 0 for all real A) also establishes the 
general result for all quantum states 




ASb^) >1(5, 


( 12 ) 


where a, /3 and 7 are x, y and z in cyclic order. 

Since the two new mode spin operators defined in Eq. ([^ satisy the standard 
angular momentum operator commutation rules, the usual Heisenberg Uncer¬ 


tainty rules analogous to (101 apply, so that spin squeezing can also exist in the 
two mode cases involving the new spin operators Jx, Jy and Jz as well. These 
uncertainty principle features also apply to multimode spin operators as defined 
in Subsection 12.21 

It should be noted that finding spin squeezing for one principal spin operator 


Jy with respect to another Jx 


any of the old spin operators Sx, Sy and 


does not mean that there is spin squeezing for 
Sz In the case of the relative phase 


eigenstate (see SubSection 3.8) Jy is squeezed with respect to Jx - however none 


of the old spin components are spin squeezed. 
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2.4.2 Alternative Spin Squeezing Criteria 

Other criteria for spin squeezing are also used, for example in the article by 
Wineland et al m . To focus on spin squeezing for compared to any orthog¬ 
onal spin operators we can combine the second and third Heisenberg uncertainty 
principle relationships to give 

+ (ASy^)) > 1 (l P + I ( 5 ,) p) (13) 

Then we may define two new spin operators via 

S'_L 1 = cos 6 Sx + sinO Sy S ±2 = — sin 9 Sx + cos 9 Sy (14) 

where 9 corresponds to a rotation angle in the xy^ plane, and which satisfy the 
standard angular momentun commutation rules [S'_l 1 , S'.l 2] = iSz^ 2 , Sz] = 
iS^i, [S'^ = iS± 2 - It is straightforward to show that (^ASx^^ + (^ASy'^'^ = 

(^ASxi^) + {aSx2^) and | (^^1) p + l (§^2) P = | (s,) |2 + | (^Sy) p so that 
(a5.2) [{aSxi^) + {aSx2^)) > i (i P + I (5^2) P) (15) 

so that spin squeezing for S'zCompared to any two orthogonal spin operators 
such as S± 1 or 2 would be defined as 



and 


(a5xi^) + (a5x 2^) > iy(|(5xi)p + |(5x2)p) (16) 


For spin squeezing in (^ASz'^J we require the spin squeezing parameter ^ to 
satisfy an inequality 



The last step is an approximation for an N particle state based on the assump¬ 
tion that the Bloch vector lies in the xy plane and close to the Bloch sphere, 
this situation being the most conducive to detecting the fluctuation (ASz"^)- 


In this situation -y/P + I (*^-*- 2 ) 1^) approximately N/2. The con¬ 
dition < 1 is sometimes taken as the condition for spin squeezing [TS] . 
but it should be noted that this is approximate and Eq. (161 gives the correct 
expression. 
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2.4.3 Planar Spin Sqneezing 

A special case of recent interest is that referred to as planar squeezing |19| in 
which the Bloch vector for a suitable choice of spin operators lies in a plane 
and along one of the axes. If this plane is chosen to be the xy plane and the 
X axis is chosen then = 0 and = 0, resulting in only one Heisen¬ 

berg uncertainty principle relationship where the right side is non-zero, namely 


(A5, 


gives 


((as,^ 


■’y 


> i 

A5 2 


Combining this with ( AS'a;^ \ ( AS”, 




V) > 0 


A5, 


> il(5. 


2. So the total spin fluctu- 
= (ASy^\ + (aSx'^') will be 


ation in the xy plane defined as (^ASpara 
squeezed compared to the spin fluctuation perpendicular to the xy plane given 
by (AJperp^\ = (aS^‘^\ if 


AS 2 

‘-^^para 


< 


5, 


and (ASperp^ 


>i|(S, 


(18) 


By minimising (AS^ 


-'para 


whilst satisfying the constraints (Sz) = (Sy) = 


0 a spin squeezed state is found that satisfies (181 with (^ASpara" 
(^ASperp^) ^ J4/3, 
on the Bloch sphere. 


J2/3, 


J for large J = N/2 [12]. The Bloch vector is 


2.4.4 Spin Squeezing in Multi-Mode Cases 

Since the multi-mode spin operators defined in Eq. satisy the standard angu¬ 
lar momentum operator commutation rules, the usual Heisenberg Uncertainty 


rules analogous to (101 apply, so that spin squeezing can also exist in the multi¬ 


mode case as well. Thus 


A^,2^ 

< 


(&> 

1 and 

(ASy^) 

>5' 

(s; 

A§,2) 

< 


(5.) 

1 and 

(aSz^) 



AS.^) 

< 


(S.) 

1 and 

(A§,2\ 

.>i| 



(19) 


for S'a;being squeezed compared to 5A and so on. 

Sinmilar alternative criteria to (16) can also be obtained, for example for Sz 


being squeezed compared to any two orthogonal spin operators such as S± i or 
2 defined similarly to (14). 
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2.5 Rotation Operators and New Modes 

2.5.1 Rotation Operators 


The new spin operators are also related to the original spin operators via a 
unitary rotation operator i?(a,/?, 7 ) parameterised in terms of Euler angles so 


that ^ ^ 

= R{a,l3,"f)S^R{a,P,j)~^ ( 20 ) 


where 

R{a,l3,j) = Rz{a)Ry(l3)Rz{'j) (21) 


with = exp{i4>S^) describing a rotation about the ^ axis anticlockwise 

through an angle (j). Details for the rotation operators and matrices are set out 
in [B] . Note that Eq. (20) specifies a rotation of the vector spin operator rather 
than a rotation of the axes, so (^ = x, y, z) are the components of the rotated 
vector spin operator with respect to the original axes. 


2.5.2 New Mode Operators 

We can also see that the new spin operators are related to new mode operators 
c and d via 

Jx = {S^c + c^d)/2 Jy = {d^c—c^d)/2i = {d^ d— c^c)/2 (22) 

where 

c = R{a, (3,^)dR{a, d = R{a, (3,j)hR{a, (23) 

For the bosonic case a straight-forward calculation gives the new mode op¬ 
erators as 

c = exp(ii7) ^cos(^) exp(^ia)a-hsin(^) exp(-^ia)S^ 
d = exp(-ii7) ^-sin(^) exp(iia)a-hcos(^) exp(-im)S^ 

(24) 

and it is easy to then check that c and d satisfy the expected non-zero bosonic 
commutation rules are [e, e^] = 1 (e = c or d) and that the total boson number 
operator is N = (d’' d + c'c). As N is invariant under unitary rotation operators 
it follows that Jl + Jl + Jl = (N/2){N/2 + 1). 
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2.5.3 New Modes 

The new mode operators correspond to new single particle states \(j)c)^ \4>d) 
where 


\(t>c) = exp(-^i 7 ) ( cos(^) exp(-^ia) |(/)„) + sin(^) exp(ita) |^f,) 


l</>d) = exp(-i 7 ) -sin(^) exp(--m) + cos(^) exp(-ia) 10^) 


(25) 

These are two orthonormal quantum superpositions of the original single particle 
states \4>a}i \4>b)i such represent an alternative choice of modes that could 

be realised experimentally. 

Eqs. (24) can be inverted to give the old mode operators via 


1 


1 


1 


= exp(--ia) cos(-) exp(--t 7 ) c - sin(-) exp(+-i 7 )d 


b = exp(+(-ia) ( sin(^) exp(^J 7 ) c +cos(^) exp{-]-i-/) d 


(26) 


2.6 Old and New Modes - Coherence Terms 


The general non-entangled state for modes d and b is given by 

'P = '^Pr'Pr®?r 


(27) 


and as a consequence of the requirement that and are physical states for 
modes d and b satisying the super-selection rule, it follows that 

{{ar)c = Tr{pi{dr) = 0 {id^r)^ = Tr(pi{d^r)=0 

{(br)^ = TrCpl(br)=Q ((?)™)^=Tr(pi(5t)-)=0 

(28) 

Thus coherence terms are zero. 

For our two-mode case we have also seen that the original choice of modes 
with annihilation operators d and b may be replaced by new modes with anni¬ 
hilation operators c and d. Since the new modes are associated with new spin 
operators (^ = x, y, z) for which the covariance matrix is diagonal and where 
the diagonal elements give the variances, it is therefore also relevant to consider 
entanglement for the case where the sub-systems are modes c and d, rather than 
d and b. Consequently we also consider general non-entangled st ates for modes 
c and d in which the density operator is of the same form as (27), but with 


p)^ —?> 'pd and p^ —>■ p^. Results analogous to (28) apply in this case, but with 


c and 6 —>■ d. 
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2.7 Quantum Correlation Functions and Spin Measure¬ 
ments 

Finally, we note that the principal spin fluctuations can be related to quantum 
correlation functions. For example, it is easy to show that 

(^ASx'^'^ = i + 2 (b'^a'^at^ + + (a'^'a)^ 

(((^^a) +((a^^) +2^(?a^ (29) 

showing that is related to various first and second order quantum 

correlation functions. These can be measured experimentally and are given 
theoretically in terms of phase space integrals involving distribution functions 
to represent the density operator and phase space variables to represent the 
mode annihilation, creation operators. 
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3 Spin Squeezing Test for Entanglement 

With the general non-entangled state now required to be such that the density 
operators for the individual sub-systems must represent quantum states that 
conform to the super-selection rule, the consequential link between entanglement 
in two mode bosonic systems and spin squeezing can now be established. We first 
consider spin squeezing for the original spin operators Sx, Sy, Sz and entangled 
states of the original modes a, b. and then for the principal spin operators Jx, 
Jy, Jz and entangled states of the related new modes c, d. We show jT] that spin 
squeezing in any spin component is a sujficiency test for entanglement of the 
two modes involved. Examples of entangled states that are not spin squeezed 
and states that are not entangled nor spin squeezed for one choice of mode 
sub-systems, but are entangled and spin squeezed for another choice are then 
presented. 

3.1 Spin Squeezing and Entanglement - Old Modes 

Firstly, the mean for a Hermitian operator 12 in a mixed state 

'p = Y.Pr'pr ( 30 ) 

R 

is the average of means for separate components 

{^) = Y.PR{n) (31) 

R 


where (^12^ —Tr^pj^^). 

Secondly, the variance for a Hermitian operator 12 in a mixed state is always 
never less than the the average of the variances for the separate components 
(see pOj l 

(A12fl2^^ (32) 


where ^A12 = Tr(p ISPI ^) with AH = 12—^12^ and ^AH = Tr{j)^ AH^. ^) 

with A1 
and 12^ 


with A12fl = 12 — (12y . To prove this result we have using (311 both for 12 


- E {(F). - F)) + E (s): - (e F.)' 

R R \ R / 

= + ( 33 ) 
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The variance result ( [3^ follows because the sum of the last two terms is always 
> 0 using the result (175) in Appendix 2 of paper 1, with Cn = (ft j , 


VCr = 



which are real and positive. 


3.1.1 Mean Values for S^, Sy and 

Next, we find the mean values of the spin operators for the product state = 
Pr^Pr 


and 




(('^7^ (a)fl - = 0 

(34) 


(35) 


for SSR compliant sub-system states using (28 1 , and thus using (31) the overall 
mean values for the separable state is 


= 0 (^,/)=0 (36) 

and 

( 5 ,) =^Pfli((|S^6)^-(a^a)^)) (37) 

Hence if either ^ 0 or (^Sy'^ 7 ^ 0 the state must be entangled. This may 

be called the Bloch vector test. This result will also have later significance. 


3.1.2 Variances for Sx and Sy 


Next we calculate ( AS 1) ,( ASy'j and {Sx) , {Sy 


(Sz) for the case 

R \ “/ R \ / R' \ " / R I R 

11 ^ \ V 

of the separable state (|27[) where = Pn® Pr- From Eqs. (Ilj) we find that 


Si = i((?)2(a)2 + 5t^^t+2tS55t + (5)2(^t)2) 

si = —\{{h‘^)'^{a)'^ — b‘^b'aa^—a^abb‘^ + {b)‘^{a^Y) (38) 

y 4 


so that on taking the trace with and using Eqs. (28) we get after applying 
the commutation rules [e,e^] = 1 (e = a or &) 
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are 


Hence using (34) for x'j and (^Sy'j we see finally that the variances 




(40) 


and therefore from Eq. (321 


> 


> 


I] Pr + 2 

R 

Y. Pr (41) 


Now using (371 for iS z) we see that 


5.)| <^PKi|((?6)^-(atS)^))| < J((^^^)j^ + (“^«)i?)) (42) 


R R 

and thus for any non-entangled state of modes a and b 


ASD- 


1 , 


S: 


> 


> 


Y Pr -Ypr 

R R 


> 0 


(43) 


Similar final steps show that (^ASyJ — ^1 ^S'zy | > 0 for all non-entangled state 
of modes a and b. 

This shows that for the general non-entangled state with modes a and b as 
the sub-systems, the variances for two of the principal spin fluctuations 

and (^ASp^ are 6ot/i greater than i| z'j |, and hence there is no spin squeezing 

for Sx compared to Sy (or vice versa). Note that as | (5'y ) | =0, the quantity 


/ / / ^ \ 




v'(l(«x.) 

I^ + I 

(*-) 

p) 


in Eq. (16) is the same as that in Eq. (11) which is used here. 


It is easy to see from (361 that 


(A5^)-i|(5,)|>0 




> 0 


(44) 
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for any non-entangled state of modes d and b. This completes the set of in¬ 
equalities for the variances of Sx and Sy. These last inequalities are of course 

trivially true and amount to no more than showing that the variances 
and ( AS y \ are not negative. 


3.1.3 Variance for Sz 

For the other principal spin fluctuation we find that for separable states 


AS‘. 


1 


fl 4 
so that using 


J (b^b- J\ T((dtd- (d^d) J (d^d- (d^d)^))^ 


(45) 


(^ASl"^ >Y,PR\i(^(b^b -^^ + (|(d^d-(d^d)j")^ (46) 


From Eq. (36) it follows that 


1 , 


ASi)-^\ s 


> 


R ' 


-I- ( (d’^d — (d’^d)^) 


R 


> 0 

? 2 \ 1 
' z) 2 

that the variances are non-negative. 


(47) 


Similarly ( AS — 5 l(S'y)| >0. Again, these results are trivial and just show 


3.1.4 No Spin Squeezing for Separable States 

So overall, we have for the general non-entangled state of modes d and b 



> 


(5.) 

1 and 

(ASy^) 


(a) 1 


A§,2) 

> 


(5.) 

1 and 

(aSz^) 

.>i| 

(^x)l 



> 


(S.) 

1 and 

(aSx^) 

.>i| 

(^y)l 

(48) 


The first result tells us that for any non-entangled state of modes d and b the 
spin operator Sx is not squeezed compared to Sy (or vice-versa). The same is 
also true for the other pairs of spin operators, as we will now see. 


3.1.5 Spin Squeezing Tests for Entanglement 

The key value of these results is the spin squeezing test for entanglement. We see 
that from the first inequality in (48) for separable states, that if for a quantum 


23 




state we find that 


AS”; 




or (asI)<1\(S, 


(49) 


then the state must be entangled for modes a and b. Thus we only need to 
have spin squeezing in either of Sx with respect to Sy or vice-versa to demon¬ 
strate entanglement. Note that one cannot have both (aS'^') < 

S 

Because ( S. 


and( ASy) < ^ 


Note that one cannot have both (^AS^J < || ^5'; 
I etc. due to the Heisenberg uncertainty principle. 


= 5, 




= 0 the second and third results in (48) merely 

show that (^ASx'^'j > 0, > 0 and > 0 for SSR compliant 

non-entangled states, it may be thought that no conclusion follows regarding 
the spin squeezing involving for entangled states. This is not the case. If for 
a given state we find that 


or 


AS: 


{as: 


< 


< 


5= 


(a5: 


A,5: 


< 


< 


s. 


(50) 


(51) 


then the state must be entangled. This is because if any of these situations 
apply then at least one oi (Sx) oi ( Sy) must be non-zero. But as we have 


seen from (36) both of the quantities are zero in non-entangled states. Thus we 


only need to have spin squeezing in either of S^ with respect to Sy or vice-versa 
or spin squeezing in either of Sz with respect to Sx or vice-versa to demonstrate 
entanglement. 

Hence the general conclusion stated in [T] , that spin squeezing in any spin op¬ 
erator Sx, Sy, Sz shows that the state must be entangled for modes a and b.. The 
presence of spin squeezing is a conclusive test for entanglement. Note that the re¬ 
verse is not true - there are many entangled states that are not spin squeezed. A 
notable example is the particular binomial state |<i)) = ((a-|- h)^\/2)^/ |0) 


for which (Sx) = N/2, (Sy) = (Sz) =0 and 


(ASx 


ASy^ 


= (ASz 


= A/4, 

p \ / p \ / p \ / p \ / p 

= 0 (see i) The spin fluctuations in Sy and Sz correspond to the 
standard quantum limit. 

This is a key result for two mode entanglement. All spin squeezed states are 
entangled. We emphasise again that the converse is not true. Not all entangled 
two mode states are spin squeezed. This important distinction is not always 
recognised - entanglement and spin squeezing are two distinct features of a two 
mode quantum state that do not always occur together. 


For the two orthogonal spin operator components (14) in the xy plane S'_li 
and S' _Li it is then straightforward to show that 


If 


ASii 




(52) 
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or 


(53) 


If (A§i2)<i|(5,)| 

that is, if S _li is squeezed compared to S _l 2 or vice versa - then the state must 
be entangled. Spin squeezing in any of the spin operator component in the xy 
plane will demonstrate entanglement. 


3.1.6 Significance of Spin Squeezing Test 


The spin squeezing test for two mode systems was based on the general form 
(27) for all separable states together with the requirement that the sub-system 
density operators and were compliant with the local particle number SSR. 
From the point of view of a supporter for applying the local particle number SSR 
if the result of an experiment is that spin squeezing has occurred, the immediate 
conclusion is that the state is entangled. On the other hand from the point of 
view of a sceptic about being required to apply the local particle number SSR for 
the sub-system states, such a sceptic would draw different conclusions from an 
experiment that demonstrated spin squeezing. They would immediately point 
out that in this case spin squeezing is not a test for entanglement. However, as 
we will now see these conclusions are still of some interest. 

To discuss this it is convenient to divide possible mathematical forms for the 
density operator into categories. Considering all general two mode quantum 
states that are compliant with the global particle number SSR, we may first 
divide such quantum states into three categories, as set out in Table 1. 


REGION 

OVERALL 

SUB-SYSTEM 

CATEGORY 


QUANTUM STATE 

QUANTUM STATE 






A 

P = '^PrPr^Pr 

R 

Both and p^ are local 

* Separable 



particle number SSR compliant 






B 

P = '^PrPr^Pr 

R 

Neither 'p^ nor p^ is local 

* Separable but 



particle number SSR compliant 

non-local (3j; 




* Entangled [l] 





C 

P ^^Pr'Pr®'Pr 

R 

N/A 

* Entangled 


Table I. Categories of two mode quantum states. 

The regions referred to are shown in Figure 2. All authors would regard 
the quantum states in Region A as being separable and those in Region C as 
being entangled - it is only those in Region B where the category is disputed. 
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Those such as [T] (local SSR supporters) who require local particle number 
SSR compliance for each sub-system state would classify the overall state as 
entangled, those who do not require this (local SSR sceptics) such as [3] would 
classify the overall state as separable but non-local. Note that no further sub¬ 
classification is needed. 

In Subsection 3.29 of paper I we show that if states of the form (561 
are globally SSR compliant, then both the sub-system states pj^ and pj^ are 
local particle number SSR compliant in general. However, we point out that 
there are special matched choices for both p^^ and along with the Pr, where 
neither p^ nor is local particle number SSR compliant even though p is 
global particle number SSR compliant. But the case where just one of Pr or 
Pr is non SSR compliant does not occur, so Region B does not need to be 
sub-divided along these lines. 



Figure 2. Categories of two mode quantum states that are global particle 
number SSR compliant. The regions A, B and C are described in Table I. 

Now let us consider quantum states for which (^ASx'^^ > ^ \ I 

> ^1 (Sz) |. Such states are clearly not spin squeezed. Firstly, we 

/ p \ / p 
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know that all states in Region A satisfy these inequalities. However, some states 
in Region B and some states in Region C may also satisfy these inequalities. 
In Figure 3 the quantum states in Region B that satisfy these inequalities are 
depicted as lying in Region D, those in Region C that do so are depicted as 
lying in Region F. Hence, if we find that the quantum state is such that spin 
squeezing does occur (as in the test of (59)) we can definitely say that it does 
not lie in Regions A, D or F. It must therefore be located in Regions E or L. 
The question is - Does this determine whether the state is entangled or not 
according to the supporters of applying the local SSR as in the definition of 
entanglement used in the present paper? The answer is that it does. This is 
because the quantum state must be located within either of Regions B or C, 
since these regions include E and L respectively. In both cases it would be 
entangled according to the definition used here [1] (see Table 1). 



Figure 3. Categories of two mode quantum states satisfying inequalities 

> || I and (^ASy'^^ > I These regions are A, D, F. 

Set unions are denoted +. Referring to Figure 2, B = D (B E and C = F (B L. 


However, the sceptics of applying the local SSR would draw a different con¬ 
clusion from the experiment that demonstrated spin squeezing (as in the test 
of (59l). They would agree that the mathematics shows that a state in Region 
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A could not demonstrate spin squeezing. Nor by assumption could states in 
Regions D or F. This means that the state must lie in either Region L or Region 
E. So from the point of view of the sceptic, either the state is entangled (if it 
lies in Region L) or the sub-system states in all separable states (Region E) 
do not comply with the local particle number SSR. The sceptic’s conclusion is 
clearly interesting - in the first case the quantum state is entangled, and hence 
may demonstrate other non-classical features, and in the second case the pos¬ 
sibility exists of finding sub-systems in states that have the unexpected feature 
in non-relativistic many body physics of having coherences between states with 
differing particle number. If there was a second experimental test that could 
show that the state was not entangled, then this would demonstrate the exis¬ 
tence of quantum states (sub-systems are themselves possible quantum systems) 
in which the particle number SSR breaks down. 

The second experiment would seem to require a test for entanglement which 
is necessary as well as being sufficient - the latter alone being usually the case 
for entanglement tests. Such criteria and measurements are a challenge, but not 
impossible even though we have not met this challenge in these two papers. 
Thus, in principle there could be a pair of experiments that give evidence of 
entanglement, or failure of the Super Selection Rule. For such investigations 
to be possible, the use of entanglement criteria that do invoke the local super¬ 
selection rules is also required. Such tests are the focus of these two papers, 
though here our primary reason is because we consider applying the local particle 
number SSR is required by the physics of non-relativistic quantum many body 
systems involving massive particles. 


3.1.7 Inequality for \(^S ^ 

Of the results for a non-entangled physical state for modes d and b we will later 


find it particularly important to consider the first of (48) 


{asI)>^\{s^)\ and {asI)>^\{s 
This is because we can show that for any quantum state 


i(n,-n.)\| <i(|(nfc)| + |(n,)|) =i(iV 


(54) 


(55) 


and the mean number of bosons in 


there is an inequality involving | (^S 
the two mode system. Note that there are some entangled states (see SubSection 
3.61 for which (ASI) and (ASI) are both greater than ^ \ (S z) |, since all 


that has been proven is that for all non-entangled states we must have both 

{ASl)>^^\{Sz)\and {asI)>1\{s 

Hence we may conclude that spin squeezing in any of the original spin fluc¬ 
tuations Sx , Sy or Sz requires the quantum state to be entangled for the modes 
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a and b as the sub-systems. Similarly, we may conclude that spin squeezing 
in any of the principal spin fluctuations Jx , Jy or requires the quantum 
state to be entangled for the modes c and d as the sub-systems, these modes 
being associated with the principal spin fluctuations via Eq. (22). Although 


finding spin squeezing tells us that the state is entangled, there are however no 
simple relationships between the measures of entanglement and those of spin 
squeezing, so the linkage is essentially a qualitative one. For general quantum 
states, measures of entanglement for the specific situation of two sub-systems 
(bi-partite entanglement) are reviewed in [21) . 


3.2 Spin Squeezing and Entanglement - New Modes 


It is also of some interest to consider spin squeezing for the new spin operators 
Jx, Jy, Jz with the new modes c and d as the sub-systems, where these spin 
operators are associated with a diagonal covariance matrix. The definition of 


spin squeezing in this case is set out analogous to that in Eq.(19). In this case 
the general non-entangled state for the new modes is given by 


"R Pr ' 


' Pr 


(56) 


with the and representing physical states for modes c and d, and where 


results analogous to Eqs. (28) apply. The same treatment applies as for spin 
operators Sx, Sy, Sz with the modes d and b as the sub-systems and leads to 
the result for a non-entangled state of modes c and d 


AJ: 




and 




(57) 


showing that neither Jx or Jy is spin squeezed for the general non-entangled 
state for modes c and d given in Eq. (24). We also have 


J.) = Y.Pr{J 


R 


= 0 


R 




= 0 


(58) 


so all the results analogous to Eqs. (48) also follow. Following similar argue- 


ments as in SubSection |3.1| we may also conclude that spin squeezing in any of 
the original spin fluctuations requires the quantum state to be entangled when 
the original modes c and d are the sub-systems. Thus the entanglement test is 


or 


or 


If 

(aj^) 


(J.) 1 

or 

("v 

<5' 

(b) 1 

(59) 

If 

i^n) 

■<i| 

(^.)l 

or 

(ajI) 

<i| 

(^.)l 

(60) 

If 

i^n) 

<i| 

{jy)\ 

or 

(Aj“) 

<i| 

(^.)l 

(61) 
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then we have an entangled state for the original modes c and d. 

The result ( |58[ ) also provides a Bloch vector entanglement test - if either 
^0 or \S^ ^ 0 the state must be entangled. 

Hence we have seen that spin squeezing - either of the new or original spin 
operators requires entanglement of the new or original modes. Which spin 
operators to consider depends on which pairs of modes are being tested for 
entanglement. 

3.3 Spin Squeezing and Entanglement - Multi-Mode Case 

As we have seen the multi-mode case involves a set of n modes with annihilation 
operators di for bosons with hyperfine component A, and another set of n modes 
with annihilation operators bi for bosons with hyperfine component B. Since 
entanglement implies a clear choice of what sub-systems are to be entangled, 
there are numerous choices possible here for the present multi-mode case. Case 
1 involves two sub-systems, one consisting of all the di modes as sub-system A 
and the other consisting of all the bi modes as sub-system B. Case 2 involves 
2n sub-systems, the Ai th containing the mode di and the Bi th containing 
the mode bi. Case 3 involves n sub-systems, the ith containiing the two modes 
di and bi. These three cases relate to entanglement causing interactions in 
differing circumstances. Case 1 might apply to cases where separable states can 
be created with all the di modes coupled together to produce states and the bi 
modes coupled together to produce states . Case 2 might apply cases where 
separable states can be created with the di and all the bi modes independent 
of each together to produce states p^^*^ 0 p^^^Case 3 might apply cases where 
separable states can be created with the di and the matching bi modes coupled 
together to produce states Cases 2 and 3 will be discussed further in 

Subsection |4.4| dealing with the entanglement test introduced by Sorensen et al 

m- 

The density operators for separable states in the three cases will be of the 
form 


bep 

= '^PrPr 



R 


Pr 

= 'Pr^'Pr Case 1 


Pr 



Pr 

'^ab(l) „ ^ab(2) „ „ '^ab(i) „ '^ab('. 


AC 

Case 


(62) 

Case 2 (63) 
3 (64) 

(65) 


Discussion of whether there is a spin squeezing test for Case 1 in the multi¬ 


mode case involves a generalisation of the theory set out in SubSection 3.1 


There is a Bloch vector entanglement test, in that if either of xj or yj is 
non-zero, then the state is entangled. We also find that spin squeezing in any 
spin component requires the state to be entangled, thus generalising the spin 
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squeezing test to the multi-mode case, for two sub-systems consisti ng of all the 
modes di and al the modes bi The details are covered in Appendix |ll[ 

For Case 2 a spin squeezing test for entanglement also be obtained. The 
test is again that spin squeezing in any spin component , Sy or confirms 
entanglement of the 2n sub-systems consisting of single modes di and bi. Fur¬ 
thermore, there is also a Bloch vector entanglement test, in that if either of 


xj or \^S yj is non-zero, then the state is entangled. As these systems can 
have quantum states with large numbers N of bosonic particles, it can be said 
that entanglement in an N particle system has occurred if spin squeezing is 
found. The proof of these tests is set out in SubSection |12.1| of Appendix 

For Case 3 there is also a spin squeezing test for entanglement, but it is 
restricted. Here the test is that spin squeezing in confirms entanglement of 
the n sub-systems consisting of pairs of modes di and bi, but the test is restricted 
to the situation where exactly one boson occupies each mode pair. No spin 
squeezing test was found for the other spin operators, nor was a Bloch vector 
entanglement test obtained. The proof of this result is set out in SubSection |12.2| 
of Appendix |12[ That no general spin squeezing test for entanglement exists 
can be shown by a counter-example. If all the N bosons occupied one mode 
pair di and bi, and the quantum state for this pair corresponded to the 

relative phase eigenstate with phase Op = 0 (see SubSection 3.8) then although 
the overall state is separable, spin squeezing in Sy compared to Sz occurs (with 
= l + ^\nN, = (i - ii) and (^S ^ f • Thus there 

is a situation where a non-entangled state for sub-systems consisting of mode 
pairs is spin squeezed, so spin squeezing does not always confirm entanglement. 

As in the previous two mode cases, having established in multi-mode cases 
that spin squeezing requires entanglement a further question then is: Does 
entanglement automatically lead to spin squeezing? The answer is no, since 
cases where the quantum state is entangled but not spin squeezed can be found 
- an example is given in the previous paragraph . Thus in general, spin squeezing 
and entanglement are not equivalent.- they do not occur together for all states. 
Spin squeezing is a sufficient condition for entanglement, it is not a necessary 
condition. 


3.4 Bloch Vector Entanglement Test 


We have seen for the general non- 
d and b that 

-entangled states 

of modes c and d or of modes 

(Jx) 

= 0 

(t.) 

= 0 

(66) 

(s.) 

= 0 

(s.) 

= 0 

(67) 
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Hence the two mode Bloch vector entanglement tests 


J.) ^ 0 
^ 0 


Sy\ ^ 0 


( 68 ) 


for modes c and d or of modes d and b. The same Bloch vector test also applies 
in the multi-mode case for Case 1, where there are just two sub-systems each 
consisting of all the modes di or all the modes bi and in Case 2, where there are 
2n sub-systems consisting of all the modes di and all the modes bi. 

From Eqs. (22) and ([^ these results are equivalent to 


dc^'^ = 0 ~ 

bd^^ = 0 


= 0 


(69) 

(70) 


Hence we find further tests for entangled states of modes c and d or of modes 
d and b 


dc^ 

bd'^ 


> 0 
> 0 


> 0 


> 0 


(71) 

(72) 


As we will see in Section these tests are particular cases with m = n = 1 
of the simpler entanglement test in Eq. (165) that applies for the situation in 


the present paper where non-entangled states are required to satisfy the super¬ 
selection rule. 


3.5 Entanglement Test for Number Difference and Sum 

There is also a further spin squeezing test involving the operator S^, which is 
equal to half the number difference ^(db — da). We note that simultaneous 
eigenststes of da and db exist, which are also eigenstates of the total number 
operator. Eor such states the variances (Ada), (Ad^), (^ASf'j and 
are all zero, which does not suggest that useful general inequalities for these 
variances would be found. However, a useful entanglement test - which does 
not require SSR compliance can be found. Eor the variance of Sz in a separable 
state we have 


AS': 






R 




(73) 


R 
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( 74 ) 


For such a separable state we also find 

^ II PrU^^Dr + {^K)r) 

R 

This leads to the useful if somewhat qualitative test that if we have a state with 
a large fluctuation in the total boson number and a small fluctuation in the 
number difference, then it must be an entangled state. If it was separable and 
is large, then must also be large. There is also the converse test 

- if we have a state with a small fluctuation in the total boson number and a 
large fluctuation in the number difference, then it must be an entangled state. 


3.6 Entangled States that are Non Spin-Squeezed - NOON 
State 

One such example is the generalised N boson NOON state defined as 

p = |$)($| 

|<i)) = cos( 


(dt)^ 


(6t)^ 

^ |O)+sin0h^|O) 

\/M 


= cost 


y,-y)+sm0 


N N 


( 75 ) 


which is an entangled state for modes d and b in all cases except where cos 0 or 
sin 0 is zero. In the last form the state is expressed in terms of the eigenstates 
for (_^)^ and 5^, as detailed in [5]. 

A straight-forward calculation gives 


S 



N 

~4 


(Syj=0 (§,)=-yCos20 

= j (as!) = ^(1-cos^20) ( 76 ) 


for > I, so that using the criteria for spin squeezing given in Eq. (11) we see 
that as (^AS'^'^ ~ 5I I — hence spin squeezing does not occur 

for this entangled state. 


3.7 Non-Entangled States that are Non Spin Squeezed - 
Binomial State 

Of course from the previous section any non entangled state is definitely not 
spin squeezed. A specific example illustrating this is the N boson binomial state 
given by 

P = l^)($l 

l < 5 >) = ^ 10 ) ( 77 ) 
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where c and d are given by Eqs. (241 with Euler angles a = —tt + y, /3 = —20 
and 7 = — tt, we find that 


c = — COS0 exp(-ix)d — sm0 exp(—-ix) b = —di 

d = sin 9 exp{-ix)d — cos 6 exp{—-ix) b =— 0,2 


(78) 


where the mode operators oi and 0,2 are as defined in [6] (see Eqs. (53) therein). 
The new spin angular momentum operators (^ = x, y, z) are the same as those 
defined in [6] (see Eqs. (64) therein) and in [6] it has been shown (see Eq. (60) 
therein) for the same binomial state as in ([77]) that 



(see Eqs. (162) and (176) therein). Hence the binomial state is not spin 
squeezed since = 1| (^Jz'j \- It is of course a minimum 

uncertainty state with spin fluctuations at the standard quantum limit. Here 


J = l/y/N. Clearly, it is a non-entangled state for modes c 


and d , being the product of a number state for mode c with the vacuum state 
for mode d. 

Note that from the point of view of the original modes d and b, this is an 
entangled state, so the question is: Is it a spin squeezed state with respect to 
the original spin operators (^ = x, y, z) ? The Bloch vector and variances for 
this binomial state are given in (see Eq. (163) in the main paper and Eq. 
(410) in the Appendix). The results include: 


(s = -y cos 20 

(^AS'l^ = y (cos^ 20 cos^ X + sin^ x) ^AS'y^ = y (cos^ 20 sin^ x + cos^ x) 

(80) 

This gives (^AS'l^ (^ASy'^ — || ^ = ,^N‘^{cos‘^ 20 — 1)^ cos^ xsin^ X > 0 

as required for the Heisenberg uncertainty principle. With x = 0 we have 
(^AS^'^ = cos^ 20 and (^ASy'^ = whilst | = ^|cos20|. As 

(^AS'^'j < 5I I there is spin squeezing in for this entangled state of 

modes d and b, though not of course for the new spin operator J ^ since this 
state is non-entangled for modes c and d. This example illustrates the need to 
carefully define spin squeezing and entanglement in terms of related sets of spin 
operators and modes. The same state is entangled with respect to one choice 
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of modes - and spin squeezing occurs, whilst it is non-entangled with respect to 
another set of modes - and no spin squeezing occurs. 

To summarise - with a physically based definition of non-entangled states for 
bosonic systems with two modes (related to the principal spin operators that 
have a diagonal covariance matrix) being the sub-systems and with a criterion 
for spin squeezing that focuses on these principal spin fluctuations, it seen that 
whilst non-entangled states are never spin squeezed and therefore although en¬ 
tanglement is a necessary condition for spin squeezing, it is not a sufficient one. 
There are entangled states that are not spin squeezed. Furthermore, as there 
is no simple quantitative links between measures of spin squeezing and those 
for entanglement, the mere presence of spin squeezing only demonstrates the 
qualitative result that the quantum state is entangled. Nevertheless, for high 
precision measurements based on spin operators where the primary emphasis 
is on how much spin squeezing can be achieved, knowing that entangled states 
are needed provides an impetus for studying such states and how they might be 
produced. 

3.8 Entangled States that are Spin Squeezed - Relative 
Phase Eigenstate 

As an example of an entangled state that is spin squeezed we consider the 
relative phase eigenstate \ ^;,0p) for a two mode system in which there are N 
bosons. For modes with annihilation operators a, h the relative phase eigenstate 
is defined as 



(81) 


where the relative phase Op = p{2tt/{N+\)) withp = —N/2, —N/2+1 ,..., -|-A^/2, 
is an eigenvalue of the relative phase Hermitian operator of the type introduced 
by Barnett and Pegg [25] (see jS] and references therein). Note that the eigen¬ 
values form a quasi-continuum over the range —tt to -I-tt, with a small separation 
between neighboring phases of 0{1/N). The relative phase state is consistent 
with the super-selection rule and is an entangled state for modes a, b. The 
Bloch vector for spin operators Sx, Sy, Sz is given by (see i) 



(82) 


but the covariance matrix (see Eq. (178) in [B]) is non-diagonal. 

3.8.1 New Spin Operators 

It is more instructive to consider spin squeezing in terms of new spin operators 
Jx, Jy, Jz for which the covariance matrix is diagonal. The new spin operators 
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are related to the original spin operators via 

X = 

Jy = sin 9p Sx + cos Op Sy 

Jz = — cos Op Sx + sin Op Sy (83) 

corresponding to the transformation in Eq. Q with Euler angles a = —tt + Op, 
j3 = — t : 12 and 7 = — tt . 


3.8.2 Bloch Vector and Covariance Matrix 


The Bloch vector and covariance matrix for spin operators Jx, Jy, Jz are given 
by (see Eqs. (180), (181) in [6] - note that the C{ Jy, Jy) element is incorrect in 
Eq. (181)) 


and 


j.x)=0 (Jy)=0 


C{Jx,Jx,) C{Jx,Jy) C{Jx,Jz) 
C{Jy,Jx) C{Jy,Jy) C {Jy , Jz) 
C{Jz,Jy) C{Jz,Jy) C{Jz,Jz) 
0 

0 
0 


J.) = -N^ 


(84) 


7 + |lniV 0 


_ 21! W2 

1^6 64 ^ 


1 


(85) 


With (AJ2) = ViV^, = 1 + ilnV and (aJ^) = (i - f^) and 


the only non-zero Bloch vector component being \ J zj = —N | it is easy to see 

that ^AJ^^ (^AJy^ > j\ (^J P as required by the Heisenberg Uncertainty 

Principle. The principal spin fluctuations in both Jx and Jz are comparable 
to the length of the Bloch vector and no spin squeezing occurs in either of 
these components. However, spin squeezing occurs in that Jy is squeezed with 

respect to Jx - (^AJy'j only increases as ginIV whilst 1| (^Jz^ \ increases as 
■^N for large N. Hence the relative phase state satisfies the test in Eq. (49) to 

~ V\^/N 


demonstate entanglement for modes c, d. Here y (^Aj‘^yJ/\ ^ J 

which indicates that the Heisenberg limit is being reached. 

Note that none of the old spin components are spin squeezed. As shown 
in Ref. [B] (^ASl^ = cos^ Op N^, sin^^p and 

(^ASl^ = , along with (^Sx'^ = f cosO. 


Sy)= -N I sinOp 


S, 


0. All variances are of order whilst the non-zero means are only of order 
N. Hence spin squeezing in one of the principal spin operators does not imply 
spin squeezing in any of the original spin operators. This is relevant to spin 
squeezing tests for entanglement of the original modes. 
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3.8.3 New Modes Operators 


To confirm that the relative phase state is in fact an entangled state for modes 
c, d as well as for the orig inal modes a, 6, we note that the new mode operators 
c, d are given in in Eq. (24) with Euler angles a = —tt + 9p, (d = —7r/2 and 
7 = —TT. The old mode operators are given in Eq. (26) and with these Euler 
angles we have 



b = 



( 86 ) 


This enables us to write the phase state in terms of new mode operators c, d as 


N 

Y’ 



\/ivTT \V 2 


-1 


N N/2 NI4:-k/2 N/4+k/2 

E E E 

k= -N/2 r= -JV/4+fe/2 s= -N/4-k/2 


,, ( i'|jV/4-/c/2+r 

y/{N/2-k)\ ^{N/2 + k)r 

{N/2-ky. {N/2 + ky. 

^ {N/4 - k/2 - r)!(fV/4 - k/2 + r)! (N/4 + fc/2 - s)!(iV/4 + k/2 + s)! 

x(^t)Af/2-(r+s) f^^^N/2+(r+s) |q^ 


(87) 


We see that the expression does not depend explicitly on the relative phase 0p 
when written in terms of the new mode unnormalised Fock states (cl')'^/^“l''+®) 
This pure state is a linear combination of product states of the form \ N/2 — m)^® 
|7V/2 + m)^ for various m - each of which is an N boson state and an eigenstate 
for J 2 with eigenvalue m, and therefore is an entangled state for modes c, d 
which is compatible with the global super-selection rule. Note that there cannot 
just be a single term m involved, otherwise the variance for J 2 would be zero 
instead of fj) We will return to the relative phase state again in 

Subsection 14.11 
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4 Other Spin Operator Tests for Entanglement 


In this Section we examine a nnmber of previously stated entanglement tests 
involving spin operators. It turns out that many of the tests do confirm entan¬ 
glement for massive bosons according to the SSR and symmetrisation principle 
compliant definition as it is defined here, though not always for the reasons 
given in their original proofs. Importantly, in some cases for massive bosons the 
tests can be made more general. 

There are a number of inequalities involving the spin operators that have pre¬ 
viously been derived for testing whether a state for a system of identical bosons 
is entangled. These are not always associated with criteria for spin squeezing - 
which involve the variances and mean values of the spin operators. Also, some 
of these inequalities are not based on the requirement that the density operators 
'PR’ 'Pr the expression for a non-entangled state conform to the super-selection 
rule that prohibits quantum superpositions of single mode states with differing 
numbers of bosons (which was invoked because they represent possible quantum 
states for the separate modes - local particle number SSR compliance). Only 
generic quantum properties of the sub-system density operators were 

used in the derivations. In contrast, our results are based in effect on a stricter 
criterion as to what constitutes a separable state., so of course we obtain new 
entanglement tests. However, entanglement tests which were based on not re- 
quiring SSR compliance for p^, p^ will also confirm entanglement when SSR 


compliance is required. This outcome occurs in the SubSection |4.I| in the case 
of the Hillery spin variance entanglement test. It also occurs in SubSection |4.2| 


for the entanglement test in (117) involving spin operators for two mode sub¬ 
systems, in SubSection 4.3 for the entanglement test in (129) involving mean 


values of powers of local spin operators, and in two entanglement tests (152) 


(154) in SubSection 4.5 that involve variances of two mode spin operators. 


Other entanglement tests have been proposed whose proofs were based on 
forms of the density operator for non-entangled states that are not consistent 
with the symmetrisation principle. The sub-systems were regarded as labelled 
individual particles, and strictly speaking, this should only apply to systems 
of distinguishable part icles. These incl ude the spin squeezing in the total spin 

In that SubSection we show that 


4.4 


operator Sz test ( 133[ ) in SubSection 
the original proof in |13) can be modified to treat identical particles but now 
with distinguishable pairs of modes as the sub-systems, but the proof requires 
that the separable states are restricted to one boson per mode pair. However, 
in SubSection |3.1| we have already shown that for two mode systems in which 
SSR compliance applies spin squeezing in Sz demonstrates two mode entan¬ 
glement. Also, in SubSection |3.3| we showed that in multi-mode cases modes 
associated with two different internal (hyperfine) components, that spin squeez¬ 
ing in Sz also shows entanglement occurs in two situations - one where there 
are two sub-systems each just consisting of modes associated with the same 
internal component (Case 1), the second where each mode counts as a separate 
sub-system (Case 2). Thus the spin squeezing in Sz test does demonstrate en- 
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tanglement for identical massive bosons, though not for the reasons given in the 
original proof. These new SSR compliant proofs now confirm the spin squeezing 
in Sz as a valid test for entanglement in two component or two mode BECs. 


4.1 Hillery et al 2006 

4.1.1 Hillery Spin Variance Entanglement Test 


An entanglement test in which local particle number SSR compliance is ignored 
is presented in the paper by Hillery and Zubairy [23] entitled ’’Entanglement 
conditions for two-mode states”. The paper actually dealt with EM field modes, 
and the density operators for photon modes allowed for coherences be¬ 

tween states with differing photon numbers. A discussion of SSR for the case of 
photons is presented in Paper I, in SubSection 3.2. Hence the conditions in 
Eq. (28) were not applied. 


We will now derive the Hillery spin variance inequalities involving ( AS 


(^AS 


3.1 


but now 


y / by applying a similar treatment to that in SubSection 
ignoring local particle number SSR compliance. It is found that for the original 
spin operators Sx, Sy, Sz and modes a and b 


-j(((P)’)„<(S)bB + ((b")„<(S*)"» (88) 

where terms such as and previously shown to be zero have 

been retained. Note that in [23j the operators Sx, Sy, Sz constructed from the 
EM held mode operators as in Eq. Q would be related to Stokes parameters 
Hence 


= + + 1)) (89) 

where the terms ((a^)^) cancel out. This is the same as before. 

However, 
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is now non-zero, since the previously zero terms have again been retained. Hence 


so that we now have 

= + 1 ) + " {^^)r {^)r 

= ^ (Hh + - I I'l ^')r I') 


But from the Schwarz inequality - which is based on ^(al — (a^))(a — (a))) > 0 
for any state 


(a)R^ < (a^a) 


R 


< 


Pb) 


/ R 


SO that 




(93) 

(94) 


and thus from Eq. (321 it follows that for a general non entangled state 


(a52) + ^ASl) > Y, Pr^ iim)^ + (n,)^) (95) 

R 

However, half the expectation value of the number operator is 

\ {^) ^ \ fl + (^a)fl) (96) 

R 


SO for a non-entangled state 

l^ASl) + [asI) >\(n) (97) 

This inequality for non-entangled states is given in [23] (see their Eq. (3)). 
The above proof was based on not invoking the SSR requirements for separable 
states that we apply in this paper. 


4.1.2 


Validity of Hillery Test for Local SSR Compliant Non-Entangled 
States 


However, it is interesting that the inequality (97) can be more readily derived 
from the definition of entangled states used in the present paper - which is based 
on local particle number SSR compliance for separable states. We would then 


find that ( S 


AS: 


= s 


= 0 and hence 






+ {Pa)j 




40 




instead of Eq.(92). Since the term always positive we find 

after applying Eq. (32) that 




(99) 


which is the same as in Eq. (97). Hence, finding that ^AS^y + (^ASy^ < 

i would show that the state was entangled, irrespective of whether or not 
entanglement is defined in terms of non-physical unentangled states. 

Thus, the Hillery spin variance entanglement test |23j is that if 


AS: 




( 100 ) 


then the state is an entangled state of modes a and b. This test is still used 
in recent papers, for example (53], which deal with the entanglement of 
sub-systems each consisting of single modes a, b for a double well situation (in 
these papers S x —t Jab^ ^ y ^ z ~’^ab)- 


4.1.3 Non-Applicable Entanglement Test Involving | zj \ 

Previously we had found for a general non-entangled state that is based on 
physically valid density operators 


(ASI 


> 

0 


'aK 


> 

0 

(101) 

SO that the sum of the variances 

satisfies the inequality 


(a52) 

+ (a5A>| 


>1 

(102) 


This is another correct inequality required for a non-entangled state as defined 
in the present paper. It follows that if only physical states p^, p^ are allowed, 

the related entanglement test involving ( AS w + ( AS" y j would be 


(a5; 




5 , 


(103) 


For any quantum state we have 


= ^1 iinb) - (Ua)) I < ^ (ins) + (fia)) = ^ (n) (104) 


which means that it is now required that ^AS'^\ -|- (^ASy^ be less than a 
quantity that is smaller than in the criterion in (97). 
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However, it should be noted (see (0) that all states, entangled or otherwise, 
satisfy the inequality 




(105) 


so the inequality in (1021 - though true, is of no use in establishing whether a 


state is entangled in the terms of the meaning of entanglement in the present 
paper. There are no quantum states, entangled or otherwise that satisfy the 


proposed entanglement test given in Eq. (103). This general result was stated 
by Hillery et al [53]. To show this we have 


(( 


)’(^ 


.)) 


AS” ,: — lAAS',,) ( AS'a; + zAA5'u ) ) > 0 (106) 

A5^^ + A^5,^ + A^^A5^^ > 0 (107) 

for all real A. The condition that this function of A is never negat ive g ives the 


Heisenberg uncertainty principle ( AS' 


AS I 


>i\ 


s, 


P and (1051 follows 

from taking A = 1 and A = — 1. Even spin squeezed states with (AS^ 


< 


^1 (S^ ) I still have ( ASl) + {ASl 


> 


so it is never found that 


\ “ 

AS^^ + (^ASy'j < I ^S 2 ^ I and hence this latter inequality cannot used as a 
test for entanglement. 

Fortunately - as we have seen, showing that spin squeezing occurs via either 
l^AS"^ < ^1 i^S I or (^ASy'^ < || z'j I is sufficient to establish that the 

state is an entangled state for modes a,b, with analogous results if principle 
spin oper ators are considered. Applying the Hillery et al entanglement test in 
Eq. (100) involving ^ ^ valid entanglement test, but is usually less 

stringent than the spin squeezing test involving either (^AS^^ < || | or 

(^ASy'^ < || |. For the Hillery et al entanglement test to be satisfied at 

least one of (^AS'^'j or (^ASy'j is required to be less than ^ whereas for 

the spin squeezing test to apply at least one of (^AS'lj or (^ASy'^ must be less 

than I- The quantity (^Sz'j \ is likely to be smaller than ^ - 

for example the Bloch vector may lie close to the xy plane, so a greater degree 
of reduction in spin fluctuation is needed to satisfy the spin squeezing test for 
entanglement. 

However, this is not always the case as the example of the relative phase state 
discussed in SubSection 13.81 shows. The results in the current SubSection can 
easily be modified to apply to new spin operators J^, Jy, Jz , with entanglement 
being considered for new modes c and d. The Hillery et al |23j entanglement 
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test then becomes 


(108) 


(aj) 


AJ; 




In the case of the relative phase eigenstate we have from Eq. (85) that (^AJ ^ 

l^AJy'^ = + i + l in A R:! for large N. This clearly exceeds ^ (^n'^ = 

|iV, so the Hillery et al [ 
have seen in SubSection 


3.8 


test for entanglement fails. On the other hand, as we 
« ^ A, so the spin squeezing 






test is satisfied for this entangled state of modes c and d. 


4.1.4 Hillery Entanglement Test - Multi-Mode Case 

It turns out that the Hillery spin variance test can also be applied in multi- 
mode situations, where the spin operators are defined as in SubSection |2.2[ As 


explained in SubSection 3.3 three cases occur in regard to specifying the sub¬ 
systems. For Case 1, where there are two sub-systems each consis ting of all the 
modes di or all the modes 6 ^. the Hillery spin variance test as in (100) applies. 


The proof is set out in Appendix |13[ and again does not require the sub-system 
density operators to be local SSR compliant. Also, for Case 2 where there 
are 2n subsystems consisting of all modes di and all modes bi the Hillery spin 
variance test as in (100) applies. The proof is set out in Appendix[^ and again 
does not require the sub-system density operators to be local SSR compliant. 
However, for Case 3 where there are n sub-systems consisting of all mode pairs 
di and bi the Hillery spin variance test does not apply. The reason is explained 
in Appendix [l^ Basically, it is because specihc sub-system density operators 
(see (|64|) could be entangled states of the modes di and bi all of which do 


satisfy the Hillery test involving \NiJ for this Ah sub-system. If we choose a 


special separable state of the form (641 with just one term (no sum over i?), it is 
easy to see that the Hillery test will be satished for the full system. However, the 
full system state involving these sub-systems is still a separable state, showing 
that satisfying the Hillery spin variance test does not always require the state 
to be entangled. 


4.2 He et al 2012 

In two papers dealing with EPR entanglement He et al [5S], [21] a four mode 
system associated with a double well potential is considered. In the left well 
I there are two localised modes with annihilation operators oi, bi and in the 
right well 2 there are two localised modes with annihilation operators 02 , 62 . 
The modes in each well are associated with two different internal states A and 
B. Note that we use a different notation to |2S], [M]. This four mode system 
provides for the possibility of entanglement of two subsystems each consisting 
of pairs of modes. We can therefore still consider bipartite entanglement how¬ 
ever. With four modes there are three different choices of such sub-systems but 
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perhaps the most interesting from the point of view of entanglement of spatially 
separated modes - and hence implications for EPR entanglement - would be 
to have the two left well modes oi, bi as sub-system 1 and the two right well 
modes 0 , 2 , 62 as sub-system 2. This is an example of the general Case 3 consid¬ 
ered for multi-modes in SubSection |3.3[ Consistent with the requirement that 
the sub-system density operators conform to the symmetrisation 

principle and the super-selection rule, these density operators will not in general 
represent separable states for their single mode sub-systems Oi, bi or 0 , 2 , &2 - 
and may even be entangled states. As a result when considering non-entangled 
states for the pair of sub-systems 1 and 2 we now have 

{(a%r) = Tr{rn^^\a%r) ^0 * = 1,2 (109) 

\ / ab{i) 


in general. In this case where the sub-systems are pairs of modes the spin 
squeezing entanglement tests as in Eqs.(59) - (611 for sub-systems consisting of 
single modes cannot be applied, as explained for Case 3 in SubSection 3.3 unless 
there is only one boson in each sub-system. Nevertheless, there are still tests of 
bipartite entanglement involving spin operators. We next examine entanglement 
tests in Refs. [2S] , [53] to see if any changes occur when we invoke the definition 
of entanglement based on SSR compliance. 


4.2.1 Spin Operator Tests for Entanglement 

There are numerous choices for defining spin operators, but the most useful 
would be the local spin operators for each well [24] defined by 

5 ^ = (pfii +a%)/2 Si = {V[ai - a\bi)/2i Sl = {VJ)i - a\ai)/2 

si = (^02-I-0262)/2 S'y = (^02 - a 2 & 2 )/ 2 * 5 '^ = ( 5|62 - a2«2)/2 

( 110 ) 


These satisfy the usual angular momentum commutation rules and those or the 
different wells commute. The squares of the local vector spin operators are 
related to the total number operators Ni = b\bi -|-ajai and N 2 = b\b 2 + ala 2 
as = ( Ni/2)iNi/2 + 1) and = ( ^ 2 / 2 )(iV 2/2 + 1). The total 

O' O' 

spin operators are 

Sa = Sl + Sl a = x,y,z ( 111 ) 

and these satisfy the usual angular momentum commutation rules. Hence there 
may be cases of spin squeezing, but these do not in general provide entanglement 
tests. 

Eor the local spin operators we have in general 





rr(p“'(") sl) ^ 0 


a = x,y,z 

( 112 ) 
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based on (1091, and applying (31) we see that in general 
states. Thus the Bloch vector 


for separable 

fest for entanglement does not apply, 
there is no spin squeezing test either. Following a similar 


Furthermore 

approach as in Section we can obtain the following inequalities for separable 
states of the sub-systems 1 and 2 


(AS|)-i|(s.)l 

> (((aS;)")- f(s;) n+Y.^K (((aS")’)- sI(»’)d 

R R 

(113) 

(AS»>-i|(s.>l 

> (((aSJ)") - h (si) I) + Cr({(as;)’) - ll (s?) I) 

R R 

(114) 


Similar inequalities can be obtained for other pairs of spin operators. In neither 
case can we state that the right sides are always non-negative, For example, 
each may be a spin squeezed state for S ^ versus S y and each may 
be a spin squeezed state for S ^ versus S' ^. In this case the right side of the first 
inequality is a negative quantity, so we cannot conclude that the total S^, is 
not squeezed versus S y for all separable states. As the and can be 
chosen independently we see that separable states for the sub-systems 1 and 2 
may be spin squeezed, so the presence of spin squeezing in a total spin operator 
is not a test for bipartite entanglement in this four mode system. This does not 
of course preclude tests for bipartite entanglement involving spin operators, as 
we will now see. 

In Subsection 2.8 of paper 1 it was shown that | fls 

for a non-entangled state of general sub-systems A and B, so with —>■ = 

Sij. — iSy and —?> = S'^ — iSy = (5*^)^ to give 

< (^SlSlSlSl'^ (115) 


for a non-entangled state of sub-systems 1 and 2. For the non-entangled state 
of these two sub-systems we have 


sI^-) = Yp« 


(s^ 

ab{l) 


ib(2) 


(116) 


which in general is non-zero from Eq.(112). 


Hence a valid entanglement test involving spin operators for sub-systems 1 
and 2 - each consisting of two modes localised in each well exists, so if 



(117) 
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then the two sub-systems are entangled. A similar conclusion is stated in 
where the criterion was predicted to be satisfied for four mode two well BEC 
systems. This test for entanglement involves the local spin operators, though it 
is not then the same as spin squeezing criteria. It is referred to there as spin 
entanglement. Other similar tests may be obtained via different choices of ^Ia 
and flB- 


4.3 Raymer et al 2003 

In a paper also dealing with bipartite entanglement where the sub-systems each 
consist of two modes, Raymer et al [27] derive entanglement tests involving 
spin operators for the sub-systems defined in (IIO). With Hermitian operators 
^IajAa and for the two sub-systems we consider 

U = ahA +13^3 V = qAa- I3^b (118) 

where a,/3 are real. Then with p = = Pr ® using 

(32) it can first be shown that 


ai/A >Y,Pr (^Ur) > E 


R 


(119) 


where At/R = C/-<^C/^ ,AVr = V-(VJ with (u) ^ = TriUp^), (v) ^ = 
TriVpR). _ 


R 


Substituting for U and V from (118) and using Pn = Pn® Pr we can then 
evaluate the various terms as follows. 


W 


U 


(( 


u 


R 


a~ {VLa/ (3 (Og ) ^-|-2 q!/3 ( Oa ) (O; 


\“^/r 


R 


= a 


AC/I 


with a similar result for (AEJ 



Here for sub-system A we define = 

R _ . /^ \ R 


Tr{n\p^), = Tt^^a'Pr) and (a^)^ = Tr(A\p^), (Aa)^ = Tr(AApi) 

with analogous expressions for sub-system B. 
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We thus have 

(aU^) > 


AW 


> oi^ ^ Pr ^ ^ Pr (^^^br'^ 


r2 \ 


R 
B 
R 
' B 


( 121 ) 


where AHar = — ('^a') , A^br = Hb — (^b') , AAar = Aa — (Aa) 


/ B 


and AAbr = Ab - ( A 


Adding the two results gives 

(^AU^'^ + (jAy2^ 




^Ir)^ 


aa 


2 \ 
br}^ 


( 122 ) 


a general variance inequality for separable states. 

This last result can be developed further based on the commutation rules 


[flA,AA]=*0A [nB,AB]=i0B 
The Schwarz inequalities - valid for all real and Ab 

(A^Iar — jA^AA^ib) Pb {AflAR + IXaAAar) 

(AIIbb — iAbAAbb) Pb {^^br + IXbAAbr) 

lead to the following inequalities 


l2 \ 


(An^B/ +-^A (0 a) + A^ ( AA^b) ^ 0 


(AVIbr 


2 \ 


/ B \ IB 


R 


+ Ab ( 0B ) + -^B ( 


=): ^ 

R 


l2 \ 


-brJb ^ 


so by taking Aa,b = 1 or —1 we have 






B 




0A 


0B 


B 


> 0 
> 0 


> 0 


(123) 


(124) 


(125) 


(126) 


The Heisenberg Uncertainty principle results l^AA\j^ > I ^0A^ P/4 

etc also follow from (125). 
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Noting that I/©a) | >|^Pr(0a) | = | (©a) | and ^ Pr |/0b\ | > 


R 


R 


I ^ Pr I = since the modulus of a sum is never greater than 

the sum of the moduli, we finally arrive at the final inequality for separable 
states 

A(af]A + + (A(aA^ - ^Ab)^) > (©a) | + /3'| (©b) | (127) 

This leads to the following test for entanglement 

'A{anA + /3hBf) + (A{aAA-PABf)<a^\(eA)\ +I3 ^\(Qb)\ (128) 


which is usually based on choices where = (3 =1. 

We now choose ^a = S]., Qb = S^, A a = Sy and Ar = as in Eq. (llOl 
along with a = /3 = 1 Here ©a = Si and Q\b = -S^.Here sub-systems A = 1, 
B = 2 consist of modes di, bi and 02 , &2 respectively. Hence if we have 

(a(§i + Sir) + (a(§i - Sir) < I ( 51 ) I + I ( 52 ) I (129) 

then bipartite entanglement is established. Note that this test did not require lo¬ 
cal particle number SSR compliance, but still will apply if this is invoked. Other 
tests involving a cyclic interchange of x, y, z can also be established, as can other 
tests where the signs within the left terms are replaced by (—, -b), (-b, -b), (— , —) 
via appropriate choices of a, (3. . These tests involve mean values of powers of 


local spin operators. Similar to tests in SubSections 4.1 4.2 this test also does 
not require SSR compliance. 


4.4 Sorensen et al 2001 

4.4.1 Sorensen Spin Squeezing Entanglement Test 

In a paper entitled ” Many-particle entanglement with Bose-Einstein conden¬ 
sates” Sorensen et al m consider the implications for spin squeezing for non- 
entangled states of the form 

P = '^^^Pr Pb® Pr® p\® ■■■ ( 130 ) 

R 

where Pr is a density operator for particle i. As discussed previously, a density 
operator of this general form is not consistent with the symmetrisation principle 
- having separate density operators Pr for specific particles i in an identical par¬ 
ticle system (such as for a BEG) is not compatible with the indistinguishability 
of such particles. It is modes that are distinguishable, not identical particles, 
so the basis for applying their results to systems of identical bosons is flawed. 
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However, they derive an inequality for the spin variance ( AS” ‘ 


A5?> > I 


+ (5 


(131) 


that applies in the case of non-entangled states. Key steps in their derivation 
are stated in the Appendix to m , but as the justification of these steps is 
not obvious for completeness the full derivation is given in Appendix!^ of the 


present paper. This inequality (1311 establishes that if 


r = 


(a^; 


< 


N 


(132) 


then the state is entangled, so that if the condition for spin squeezing analogous 


to that in Eq. (16) is satisfied, then entanglement is required if spin squeezing 
for S' 2 to occur. Spin squeezing is then a test for entanglement in terms of their 
definition of an entangled state. 

If the Bloch vector is close to the Bloch sphere, for example with ( S x) =0 


and ( S y) = N/2 then the condition (1321 is equivalent to 


AS; 


< 2l(^^ 


(133) 


which is the condition for squeezing in S 2 compared to S a,. Spin squeezing is 
then a test for entanglement in terms of their definition of an entangled state. 


Note that the condition (1331 requires the Bloch vector to be in the xy plane 


and close to the Bloch sphere of radius N/2. By comparison with (16) we see 
that the Sorensen spin squeezing test is that if there is squeezing in Sz with 
respect to any spin component in the xy plane and the Bloch vector is close to 
the Bloch sphere, then the state is entangled. 

As explained above, the proof of Sorensen et al really applies only when the 
individual spins are distinguishable. It is possible however to modify the work of 
Sorensen et al [13] to apply to a system of identical bosons in accordance with 
the symmetrization and super-selection rules if the index i is re-interpreted as 
specifying diffferent modes, for example modes localised on optical lattice sites 
i = l,2,..,n or distinct free space momentum states listed i = l,2,..,n. On 
each lattice site or for each momentum state there would be two modes a, 6 - 
for example associated with two different internal states - so the sub-system 
den sity operator then applies to the two modes on site i. However the proof 
of (131) requires the to be restricted to the case where there is exactly 
one identical boson on each site or in a momentum state. Such a localisation 
process in position or momentum space has the effect of enabling the identical 
bosons to be treated as if they are distinguishable. Details are given in the next 
Subsections. A similar modification has been carried out by Hyllus et al [29]. 
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However, as we have seen in SubSection |3.2| it does in fact turn out for 
two mode systems of identical bosons that showing that is spin squeezed 
compared to S x or S y is sufficient to prove that the quantum state is entangled. 
There are no restrictions either on the mean number of bosons occupying each 
mode. The proof is based on applying the requirement of local particle number 
SSR compliance to the separable states in the present case of massive bosons 
and treating modes (not particles) as sub-systems. In SubSection 3.3 we have 
also shown that the same result applies to multi-mode situations in cases where 
the sub-systems are all single modes (Case2) or where there are two sub-systems 
each containing all modes for a single component (Casel). So the spin squeezing 
test is still valid for many particle BEC, though the justification is not as in the 
proof of Sorensen et al m - which was derived for systems of distinguishable 
particles, with each sub-system being a single two state particle. 


4.4.2 Revising Sorensen Spin Squeezing Entanglement Test - Lo¬ 
calised Modes 

The work of Sorensen et al really applies only when the individual spins are 
distinguishable. It is possible however to modify the work of Sorensen et al m 
to apply to a system of identical bosons in accordance with the symmetrisation 
and super-selection rules if the index i is re-interpreted as specifying diffferent 
modes, for example modes localised on optical lattice sites or in different mo¬ 
mentum states i = l,2,..,n. Another example would be single two state ions 
with each ion being trapped in a different spatial region. The revised approach 
draws on the results established for multi-mode cases in Appendices and 

With two single particle states a, b available on each site (these could be two dif¬ 
ferent internal atomic states or two distinct spatial modes localised on the site) 
the modes would then be labelled with a = a, b. The mode orthogonality 
and completeness relations would then be 

= 1 ( 134 ) 

Oii 

With the particles now labelled K = 1,2,3, ...one can define spin operators in 
first quantization via 

K i 

K i 

= EE(l<(>^>*(^))('(’b.(^)|-ka,(l^))('^a,(l^)l)/2 (135) 

K i 

In second quantization if the annihilation, creation operators for the modes \4>ai) 
,\4>bi) 8-re di, bi and d|, b\ respectively, then the Schwinger spin operators are 


12 
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just 




S,, 


J2{^a^+a%)/2 = J2Sl 

i i 

J2i^a^-a%)/2^ = J2Sl 

i i 

J2{b%-ala.)/2 = Y,Sl 

i i 


( 136 ) 


It is easy to confirm that the overall spin operators Sa and the spin operators 
S'), for the separate pairs of modes (o^' for short) satisfy the 

same commutation rules as Sorensen et al m have for the overall spin operators 
and those for the separate particles. With this modification the non-entangled 


state in Eq. (1301 could be interpreted as being a non-entangled state where the 


subsystems are actually pairs of modes \4>ai) ■\4'bi) density operators 

'pji would then refer to a subsystem consisting of these pairs of modes. This 
corresponds to Case 3 discussed in SubSection |3.3[ It is to be noted that entan¬ 
glement of pairs of modes is different to entanglement of all separate modes - 
Case 2 discussed in SubSection |3.3| It is an example of a special kind of multi- 
mode entanglement - since the modes may themselves be entangled 

we may have ’’entanglement of entanglement”. In terms of the present paper 
the density operators would be restricted by the super-selection rule to sta¬ 
tistical mixtures of states with specific total numbers Ni of bosons in the pair 
of modes |(/),jj) ,|^i>bi). In terms of Fock states \nai ), [n-bi) for this pair of modes 
the allowed quantum states for the sub-system will be 


Ni 




(137) 


SO at this stage the general mixed physical state for the two mode system could 
be 


oo Ni Ni 

?R= E E E^f «rifc)a.aL 

Ni=0 k=0l=0 


\N,-k),^{N,-l\,^ (138) 


This state has no coherences between states of the two mode subsystem with 
differing total boson number Ni for the pair of modes. However this is still an 
entangled states for the two modes \cj)^^) ,|0f,j), so the overall state in Eq. (1381 


is not a separable state if the subsystems were to consist of all the distinct 
modes. 


4.4.3 Revising Sorensen Spin Squeezing Entanglement Test - Sepa¬ 
rable State of Single Modes 

It is possible however to link spin squeezing and entanglement in the case where 


the sub-systems consist of all the distinct modes (Case2 in SubSection 3.3). 
To obtain a fully non-entangled state of all the modes \ 4 >ai) J'/’b*) the density 
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operator 'p^ must then be a product of density operators for modes \4>ai) 

PH = rn^^PH (139) 

giving the full density operator as 

P = '^Pr (Tr ® Pfl ) 0 (pn^ ® ® (p’r (140) 

R 

as is required for a general non-entangled state all 2N modes. Furthermore, 
as previously the density operators for the individual modes must represent 
possible physical states for such modes, so the super-selection rule for atom 
number will apply and we have 

= rr(p-(a,)")=0 ((«!)") =Tr(p“*(hl)")=0 

= Tr{p^R{hr) = 0 ((6|)-)^^ = Tr(/^*(&l)™) = 0 

(141) 


The question is whether this reformulation will lead to a useful inequality 
for the spin variances such as This issue is dealt with in Appendix 

|12| and it is found that we can indeed show for the general fully non-entangled 


state (1401 that 




^ 2 ' 




and 


(As;)>i| 


(142) 


or S y then the state must 


This shows that if there is spin squeezing in either S 
be entangled. Note that this result depends on the general non-entangled state 
being non-entangled for all modes and that the density operator for each mode 
a-i or bi being a physical state with no coherences between mode Fock states with 
differing atom numbers. In terms of the revised interpretation of the density 
operator to refer to a multi-mode system with modes {(fai) il'/’&i) th® statement 
that spin squeezing for systems of identical massive bosons requires all the modes 


to be entangled is correct. However superposition states of the form (137) that 


are consistent with the super-selection rule applying to pure states of a two mode 
system are precluded, and such states ought to be allowed if entanglement of 
pairs of modes rather than of separate modes is to be considered. 

In addition, we can show that if either or is non-zero, then 

the state must be entangled - the Bloeh vector test. Finally, if it is found that 
if there is spin squeezing in S'z then the state must be entangled. Thus spin 
squeezing in any spin component confirms entanglement of the 2n individual 
modes. 


4.4.4 Revising Sorensen Spin Squeezing Entanglement Test - Sepa¬ 
rable State of Pairs of Modes with One Boson Occupancy 


It is also possible however to link spin squeezing and entanglement in the case 
where the subsystems consist of pairs of modes (CaseS in SubSection 3.3), but 
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only if further restrictions are applied. The general non-entangled state of the 


pairs of modes would actually be of the form (see (64), here the ab dropped for 
simplicity) 

" ' ^ " (143) 


p=y^pRPR' 


R 


'Pr^ 


'Pr' 


where the are now of the form given in Eq. (1381 and no longer are density 
operators for the ith identical particle. Unlike in (141) we now have expectation 


values ((oi)”) - = Trlfpf^fdi)'^) etc that are non-zero, so considerations of the 
link between spin squeezing and entanglement - now entanglement of pairs of 
modes, willl be different. 

If the density operators p)^ associated with the pair of modes a^, hi are all 
restricted to be associated with one boson states then this density operator is 
of the form 


Pr 


PL(|l).a(lL®|0).5 {OU + Plbi\l)^a 

+Pba(|0)m (l|ia ® |l)i& (^U) + Pbh(|0)i 


, ® |0)ib (l|ib) 
la ® Wib (llib) 


(144) 


where the p^^ are density matrix elements. With this restriction the pair of 

modes at, bi behave like distinguishable two state particles, essentially the case 
that Sorensen et al [13] implicitly considered. The expectation values for the 
spin operators S'*, S* and S* associated with the ith pair of modes are then 


SI 


2 (,Pab 
2 ^Phb Paa) 


Pba) ~ ~2i ~ 


(145) 


If in addition Hermitiancy, positivity, unit trace Tr{ffjf) = 1 and Trifp'jff' < 
1 are used (see Appendix 14) then we can show that p),^ and are real and 
positive, pi,, = {pU* and p*„^p^(, - |p*„j2 > 0. The condition Tr(p)^) = 1 leads 
to p)j„ -I- p),j = 1, from which Trlfp^f}'^ < 1 follows using the previous positivity 
results. These results enable the matrix elements in (144) to be parameterised 
in the form 


Pftb = cos a^ 

Pab = V C(i cos^ ai sin^ /?, exp{+icj)i) p\,^ = \J sin^ cos^ Ui sin^ / 3 j exp(— 

(146) 


where a^, Pi and pi are real. In terms of these quantities we then have 


SI) ^ = ^sin2Q;j sin^/3j cospi i^Si^ ^ = ^sin2ai sirP P, sin(()j 


R 


1 


SI) = - cos 2ai 

R 2 


(147) 
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(148) 


and then a key inequality 


SI 


+ 




1 

4 


1 

4 


sin^ 2ai (1 


sin4ft)< 


1 

4 


follows. This result depends on the density operators being for one boson 
states, as in (144|. The same steps as in Sorensen et al m (see Appendix [I^ I 


leads to the result 




S'. 


+ (S 


(149) 


for non-entangled pair of modes ai, bi. Thus when the interpretation is changed 
so that are the separate sub-systems are these pairs of modes and the sub¬ 
systems are in one boson states, it follows that spin squeezing requires entan¬ 
glement of all the mode pairs. 

A similar proof extending the test of Sorensen et al m to appply to systems 
of identical bosons is given by Hyllus et al [23] based on a particle entanglement 
approach. In their approach bosons in differing external modes (analogous to 
differing i here) are treated as distinguishable, and the symmetrization principle 
is ignored for such bosons. 


4.5 Benatti et al 2011 


In earlier work Toth and Gunhe [18j derived several spin operator based inequal¬ 
ities for separable states for two mode particle systems based on the assumption 


that the particles were distinguishable. As in Eq.( 130), the density operator was 


not required to satisfy the symmetrisation principle. Tests for entanglement in¬ 
volving the mean values and variances for two mode spin operators resulted. 
Subsequently, Benatti et al [30] considered whether these tests would still ap¬ 
ply if the particles were indistinguishable. Their work involves considering states 
with N bosons. 

For separable states they found (see Eq.(lO)) that for three orthogonal spin 
operators J„i, J „2 and J „3 


?2 


nl 


+ ( '^2 


+ {j, 


^ \^N{N + 2) 


n3 / — 


(150) 


from which it might be concluded that if the left side exceeded iV(iV-|-2)/4 then 
the state must be entangled. However, since = 5^ + Sy+ = 

N{N + 2)/4 the left side is always equal to N{N + 2)/4 for all states with N 
bosons, so no entanglement test results. This outcome is for similar reasons 
as for the failed entanglement test -I- (^ASy'^ < \ (^S \ discussed in 

Subsection 14.11 

For separable states they also found (see Eq.(ll)) that for three orthogonal 
spin operators J„i, J „2 and J „3 


A 7^ 




N 

> — 
- 2 


(151) 
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so that if 


(152) 




AJ^2 


A 7^ 


N 

<T 


then the state mu st be entangled. This test is an extended form of the Hillery 
spin variance test (100). To prove this result we note that + 


(aj^3) = {asi)+{asi)+{as^^) =.(^sl)+{s^^)+{s^,)-{s,y-{syy- 

Sz\ = N{N + 2)/4 — (Sx') — (Sy\ — (Sz') for all states with N bosons. 


For separable states we have = 0 so that + l^AJ'^ 2 '^ + 

l^AJ^y = N{N + 2)/4 — (^Sz'^ ■ As the eigenvalues for Sz lie between —N/2 

and +A/2 we have (^Sz'j < N'^/i. Thus + (^AJ^y + (^AJ^y > A 

as required. The test in ( fT^ is quite useful in that it applies to any three 
orthogonal spin operators, though it would be harder to satisfy compared to the 
Hillery spin variance test because of the additional .(^ASy term. 

For separable states they also found (see Eq.(13)) that for three orthogonal 
spin operators J„i, J „2 and Jn 3 


?2 


(A-l)(( 


A.^i) 


AJ, 
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n2 


)-( 


y2 \ ^ NiN-2) 




SO that if 


(A - 1 ) {{Azy +(aj^^)) - {zy < 


N{N -2) 


(153) 


(154) 


then the state must be entangled. To p rove this result for nl = it, n2 = 

Z and n3 = we use the result (41) for separable states that (^AS^ + 

(aA^) > + + ERPHi{a^^(b^yZ = + 

Pr straightforward to show that 5'f = (6'^6+a'l'a)^/4— 

b^b.ata, so that ^y = A^/4 — J2rPr ({^^Zr Hence for separable 

states iN-l)yAZ) + {ASy)-{sy > (A-l)A/2+(A-l) n>~ 

AV4 + ErPr {(^^r y^y^- Thus (A - 1)((A52) + (a^^^) - (^Z) > 

+ N J2rPr y^y second term on the right side is 

always positive the required inequality follows. 

Finally, they considered another inequality (see Eq. (12)) found to apply for 
separable states involving distinguishable particles in Ref. |18j . 

((^i) + (^2))-y-(A-l)(AJ^3)<0 (155) 
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so the question is whether an entanglement test ~ ~ 

(TV — 1) > 0 applies for the case of indistinguishable particles. For the 

case where nl = ^, n2 = ^ and n3 = — TV/2 — (TV — 

1) = (( 52 ) + + ( 52 )) - N/2 - (A) (52) + (A - 1) = 

A(A + 2)/4 - A/2 - A (aV 4 - ((a^a)^ + (^ " 1) 

As (S,y < AV4 we see that (^(s^ + - N/2 - (A - 1) < 

A^^ Pfj ((a’l’a)^ ), which is certainly > 0 and not < 0 as required. 

However, perhaps an entanglement test such that if is could be shown that 

{{sl) + {§ 1 )) - N/2 - (A - 1) (a^^) >nJ2Pr (156) 

R 


always applies then it could be included that the state is entangled. Unfortu¬ 
nately the right side could be too large for the left side to always exceed the 
right side for some separable states. Noting that -I- (b'^b^ = A for the 

A bosons states being considered we find that the right side is maximised when 
(ala)^ = (b^b'^ = N/2 for all Pj^, giving a maximum for the right side of 

A^/2 - and this can occur for some separable states. To show that the state is 
entangled the left side must exceed this value, otherwise the state might be one 
of the separable states. However, the left side is at most of order A^ from the 
firsr two terms and the negative terms on ly m ake the left side smaller. Hence 


there is no entanglement test of the form (156).(ala)^ = (b^bj 


/ R 


Hence Benatti et al [30] have demonstrated two further entanglement tests 


(152) and (154) for two mode systems of identical particle that involve spin 


operators. Again, these tests do not involve invoking the local particle number 
SSR for separable states. 


4.6 Sorensen and Molmer 2001 

In a paper entitled ’’Entanglement and Extreme Spin Squeezing” Sorensen and 
Molmer |31j first consider the limits imposed by the Heisenberg uncertainty 
principle on the variance considered as a function of | ^S' 2 ^ | for states 

with A two mode bosons where the spin operators are chosen such that (^S y = 

S y'j =0. Note that such spin operators can always be chosen so that the Bloch 
vector does lie along the z axis, even if the spin operators are not principal spin 
operators. Their treatment is based on combining the result from the Schwarz 
inequality 




+ \'S'2) ^T(T-I-1) 


(157) 
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where J = -/V/2, and the Heisenberg uncertainty principle 






(158) 


where ^ > 1. In fact two inequalities can be obtained 



(159) 


< ;.<^U(J+l)-(5.) +W J(J+l)-(5,) -^(5 


(160) 


which restricts the region in a (^ASlj versus | (^S zj \ plane that applies for 
states that are consistent with the Heisenberg uncertainty principle. Note that 
in the first inequality the minimum value for occurs for ^ = 1, and in 

the second inequality the maximum value for also occurs for ^ = 1 - the 

minimum HUP case. The first of these two inequalities is given as Eq. (3) in 
m- For states in which S'a: is squeezed relative to Sy the points in the^AS 

versus | ( S z ) | plane must also satisfy 


(AS^ 


< 


1 , 


(161) 


Note that as is a spin angular momentum component we always have 
S z ) I < d, which places an overall restriction on | (S z ) |. However, for 


^ > 1 there are values of | (S z ) | which are excluded via the Heisenberg un¬ 


certainty principle, since the quantity ^J(J -|- 1) — (S, 
becomes negative. This effect is seen in Figure 4. 




then 
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The question is: Is it possible to find values for and | | in 

which all three inequalities are satisfied? The answer is yes. Results showing 
the regions in the (AS ^ ) versus | (S' 2 ) | plane corresponding to the three 


inequalities are shown in Figures 4 and 5 for the cases where J = 1000 and with 
^ = 1.0 and ^ = 10.0 respectively. The quantities for which the regions are shown 
are the scaled variance and mean ( AS / J and | ( S 2 ) \/J, with ( AS f 


as a function of | ^S 2 y | via (1591, (160) and (161). The spin squeezing region 
is always consistent with the second Heisenberg inequality (160) and for large 


given 


J = 1000 there is a large region of overlap with the first inequality (159). For 


small J and large ^ the region of overlap becomes much smaller, as the result in 
Figure 6 for J = 1 and with ^ = 10.0 shows. As the derivation of the Heisenberg 
principle inequalities is not obvious, this is set out in Appendix fTS] 


(a) 


MITA n iJ far SQIJMrT 



OflTAJX rJ far HlTMBaXMrrJ* 1000 
0 «- 


(b) 




Figure 4. Regions in the < AS^ > versus | < 5'2 > | plane (shown 
states that satisfy (a) the spin squeezing inequality Eq. (161) (b) 


shaded) for 
the smaller 
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Heisenberg uncertainty principle inequality Eq. (1591 and (c) the larger HUP 
inequality Eq. (1601. The case shown is for J = 1000 and HUP factor ^ = 1. 
Both < > and | < > | are in units of J. The spin operators are chosen 

so that < Sx >=< Sy >= 0. 
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(c) 


DU.TA ix-tJ MTMAxajMrr J-torn 



Figure 5. As in Figure 4, but with J = 1000 and HUP factor ^ = 10.0. 
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Ml TAIX 


hrSOUMrr 




Figure 6. As in Figure 4, but with J = 1 and HUP factor ^ = 10.0. 


Sorensen and Molmer 1311 also determine the minimum for ( AS" J ) = i St 


as a function of 


5 , 


= 5 


for various choices of J, subject to the constraints 
= 0. The results show again that there is a region in the 

versus | \ plane which is compatible with spin squeezing. 

So although these considerations show that the Heisenberg uncertainty prin¬ 
ciple does not rule out extreme spin squeezing, nothing is yet directly determined 
about whether the spin squeezed states are entangled states for modes a, 6, 
where the Sa are given as in Eq. Q. The discussion in [ST] regarding entan¬ 
glement is also based on using a density operator for non-entangled states as in 
Eq. (130) which only applies to distinguishable particles (see SubSection 4.4). 
Sorensen m also showed that for higher J the amount of squeezing attainable 
could be greater. This fact enables a conclusion to be drawn from the measured 
spin variance about the minimum number of particles that participate in the 
non-separable component of an entangled state |32j . 
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5 Correlation Tests for Entanglement 

In Subsection 2.4 of the accompanying paper I it was shown that for separable 
states the inequality | (^TIa 0 0 applies, so that if 

I 0 0 (162) 

then the state is entangled. This is a general correlation test. 

As will be seen the correlation tests can be re-expressed in terms of spin 
operators when dealing with SSR compliant states. 


5.1 Dalton et al 2014 

5.1.1 Weak Correlation Test for Local SSR Compliant Non-Entangled 
States 

For a non-entangled state based on SSR compliant for modes a and b 

where the SSR is satisfied we have with Qa = (a)"* and = (&)" 


((a)- (Vr) =EPr ((a)™ (b^r)^ = EPr ((a)™)fl ((^^)”)^ = 0 (163) 


since from Eqs. analogous to (28) ((d)—)^ = ((^^)") = 0- Hence for a SSR 

compliant non-entangled state as defined in the present paper the inequality 
becomes 

0 < ((dt)-(d)- (6^)”(6)”') (164) 


which is trivially true and applies for any state, entangled or not. 

Since ((d)— (6'^)"^ is zero for non-entangled states it follows that it is merely 
necessary to show that this quantity is non-zero to establish that the state is 
entangled. Hence an entanglement test [1] in the case of sub-systems consisting 
of single modes d and b becomes 


|((d)-(?)'^)p>o 


(165) 


for a non-entangled state based on SSR eompliant p^. Note that for globally 
compliant states ((d)— (6^)"^ = 0 unless n = to, so only that case is of interest. 
This is a useful weak correlation test for entanglement in terms of the definition 
of entanglement in the present paper. A related but different test is that of 
Hillery et al [23] - discussed in the next SubSection. 

For the case where n = to = 1 the weak correlation test is 

I ((d6'^^ p > 0 (166) 


which is equivalent to 



7 ^ 0 and/or 



0, the Bloch vector test. 


62 



5.2 Hillery et al 2006, 2009 

5.2.1 Hillery Strong Correlation Entanglement Test 

In a later paper entitled ’’Detecting entanglement with non-Hermitian opera¬ 
tors” Hillery et al [33] apply other inequalities for determining entanglement 
derived in the earlier paper |23j but now also to systems of massive identical 
bosons, while still retaining .density operators that contain coherences 

between states with differing boson numbers. In particular, for a non-entangled 
state the following family of inequalities - originally derived in [23] . is invoked. 


(a)’" (Vr) I" < ((a^nar (b^riby 


(167) 


This is just a special case of (1621 with = (a)™ and Ds = (5)". Thus if 

(a)™ (Vr) p > ({a^riar (168) 


then the state is entangled.The Hillery et al [23] entanglement test (1681 is a 


valid test for entanglement and is actually a more stringent test than merely 
showing that | ({a)'^ (^^)"^ P > 0, since the quantity | ^(a)’” (fo^)"^ P is now 
required to be larger. In a paper by He et al |24| (see SubSection |5.3| ) the Hillery 
et al [23] entanglement test | ^(a)™ (&^)"^ P > ^(a'^)'"(a)™ (&'*')"(&)"'^ is applied 
for the case where A and B each consist of one mode localised in each well of a 
double well potential. This test whilst applica ble co uld be replaced by the more 
easily satisfied test | ^(q )™ (&^)"'^ P > 0 (see (^165|). However, as will be seen 
below in SubSection |5.3[ the Hillery et al [53j entanglement criterion is needed 
if the sub-systems each consist of pairs of modes, as treated in [26] , [24] . 

Note that if n m the left side is zero for states that are globally SSR 
compliant. In this case we can always substitute for two mode systems 


(a?) 

(St)"(S) 


(6t)"(6) 


{S, - iSyr 

(SSt)"= (1+^-5. 



(169) 


to write both the Hillery and the weak correlation test in terms of spin operators. 

A particular case for n = m = 1 is the test | p > (fhanf) for an 

entangled state. To put this result in context, for a general quantum state and 
any operator fl we have and ~ ~ {^))) ~ 

hence leading to the Schwarz inequality | p = | ^ Taking 
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ri = a5'f leads to the inequality < (no(n{, + l)), whilst choosing 

n = ba^ leads to the inequality < ((uq + I)??;,) for all quantum 

states. In both cases the right side of the inequality is greater than {naUb), so 
if it was found that | p > {naUb) (though of course still < {ria (ub + 1)) 

and < {{fia + l)u&)) then it could be concluded that the state was entangled. 
However, as we will see the left side | (a&W P actually works out to be zero if 


physical states for are involved in defining non-entangled states, so that 

for a non-entangled state defined as in the present paper the true inequality 
replacing | ^ < {ria fib) is just 0 < {fia fib)-, which is trivially true for any 

quantum state. 


For the case where n = m = 1 we can write the test (1681 in terms of spin 


operators using ab^ = Sx — iSy as 


2 /^ \ 2 1 
Sx) +(5„> > T 


N^) - 


st) 


(170) 


which when combined with the general result Sx^ + Sy^ + Sf^ = (A^/2)(iV/2-|-l) 
leads to the test 

+ (171) 


This is the same as the Hillery spin variance test (100) , so the Hillery first 


order correlation test does not add a further test for demonstrating non-SSR 
compliant entanglement. The Hillery correlation test for n = 2 leads to complex 
conditions involving higher powers of spin operators. 


5.2.2 Applications of Correlation Tests for Entanglement 

As an example of applying these tests consider the mixed two mode coherent 
states described in Appen dix [T^ whose density operator for the two mode a, 
b system is given in Eq. (336). We can now examine the Hillery et al [55] 


entanglement test in Eq.(168) and the entanglement test in Eq.(165) for the 


case where m = n = 1. It is straight-forward to show that 


(a’^a) {b^b)^ = jal' 


(172) 


so that I p = ^(a'^a) . A non-entangled state defined in terms of 

the SSR requirement for the separate modes satishes | p = 0, whilst for 

a non-entangled state in which the SSR requirement for separate modes is not 
specifically required merely satishes | p < ^(a'^a) {b''^b)'j. Hence the test 
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for entanglement of modes A, B in the present paper | p > 0 is satisfied, 

whilst the Hillery et al [33] test | is not. 

In terms of the definition of non-entangled states in the present paper, the 


mixture of two mode coherent states given in Eq.(336) is not a separable state 


not a separable state. As discussed in Paper 1 (see Section 3.4.3) this is be¬ 
cause a coherent state gives rise to a non-zero coherence ((d) 7 ^ 0 ) and thus 
cannot represent a physical state for the SSR compliant states involving iden¬ 
tical massive bosons (as in BECs). However, in terms of the definition of non- 
entangled states in other papers such as those of Hillery et al [53|, [33] the 
mixture of two mode coherent states would be a non-entangled state. It is thus 
a useful state for providing an example of the different outcomes of definitions 
where the local SSR is applied or not. 

A further example of applying correlation tests is provided by the NOON 
state defined in ( |75[ ) where here we consider modes A, B. All matrix elements of 

the form ((d 1 ')™(d)™ ( 6 t)"( 6 )"\ are zero for all m, n because both terms contain 


one mode with zero bosons. Matrix elements of the form ^(d)™ (&^)"y are all 
zero unless m = n = N and in this case 


{{sf(bY) = ((i 

= cose sine {N,0\{S-)^ \0,N) 

= cosO sin e VnV 2 N - 2V3N - 6V4:N - 12...//V(173) 


which is non-zero in general. Hence | ^(d)'^ p > 0 as required for both 

the weak and strong correlation tests, confirming that the NOON state is en¬ 
tangled. Carrying out this entanglement test experimentally for large N would 
involve measuring expectation values of high powers of the spin operators 
and Sy, which is difficult at present. 


5.3 He et al 2012 

For the four mode system associated with a double well described in SubSection 


I {{a^r (b]r) P < {(.airia^r i^ribyr) m) 

that apply for two non-entangled sub-systems A and B can now be usefully 

applied, since in this case the quantities ^(d^)™ (SJ)"^ are in general no longer 
zero for separable states. Thus there is an entanglement test for two sub-systems 
consisting of pairs of modes. If 

\{{d.r(b]r)\^ > {{airia^n^nbyr) 

foranyofij = 1,2 (175) 


4.2 (see [24]), the inequalities derived by Hillery et al [33] (see SubSection 5.2) 
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then the quantum state for two sub-systems A and B - each consisting of two 
modes localised in each well - is entangled. Again only the case where m = n is 
relevant for states that are global SSR compliant. 


6 Quadrature Tests for Entanglement 


In this Section we discuss tests for two mode entanglement involving so called 
quadrature operators - position and momentum being particular examples of 
such operators. These tests are distinct from those involving spin operators or 
correlation tests - the latter have been shown to be closely related to spin opera¬ 
tor tests. The issue of measurement of the quadrature variances involved in these 
tests for the case of two mode systems involving identical massive bosons will be 
briefly discussed in Section Again we have a situation where tests derived in 
which local particle number SSR compliance for separable states is ignored are 
still valid when it is taken into account. However, when local particle number 
SSR compliance for separable states is actually included new entanglement tests 
arise. The two mode quadrature squeezing test in (203) is an example, though 


this test is not very useful as it could be replaced by the Bloch vector test. The 


quadrature correlation test in (197) also applies and is equivalent to the Bloch 
vector test. However the non-existent quadrature variance test in (185) is an 


example where there is no generalisation of the previous entanglement test (see 


(177)) that applied when the SSR were irrelevant. 


6.1 Duan et al 2000 

6.1.1 Two Distinguishable Particles 

A further inequality aimed at providing a signature for entanglement is set out in 
the papers by Duan et al [55] , Toth et al [33] . Duan et al [55] considered a general 
situation where the system consisted of two disinguishable sub-systems A and B, 
for which position and momentum Hermitian operators xa,Pa and Xb,Pb were 
involved that satisfied the standard commutation rules [xajPa] = = * 

in units where h = 1. These sub-systems were quite general and could be two 
distinguishable quantum particles A and B, but other situations can also be 
treated. An inequality was obtained for a two sub-system non-entangled state 
involving the variances for the commuting observables xa + xb and pa—Pb 

(A(ia+ ®b)^) + (A(pa> 2 (176) 

which could be used to establish a variance test for entangled states of the A 
and B sub-systems, so that if 

(A(x^ + xb)2) + (A(p^-Pb)') <2 (177) 

then the sub-systems are entangled. For the case of distinguishable particles 
such states are possible - consider for example any simultaneous eigenstate of the 
commuting observables xa + xb and pA —PB- For such a state (A(xa + xbY') 
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and {A{pa — PbY) are both zero, so the simultaneous eigenstates are entangled 
states of particles A, B. For simplicity we only set out the case for which a = 1 
in [28]. T he p roof given in |28| considered separable states of the general form 
as in Eq.(]^ for two sub-systems but where and are possible states 
for sub-systems A, B. Consequently, a first quantization case involving one 
particle states could be involved, where super-selection rules were not relevant. 
As explained in the Introduction, the two distinguishable quantum particles are 
each equivalent to a whole set of single particle states (momentum eigenstates, 
harmonic oscillator states, ..) that each quantum particle can occupy, and 
because both p;^ and p^ represent states for one particle we have [ua^'Pr] = 
[nB,p^] = 0. Because p represent a state for the two particles [ua + nB,^] = 0, 
the SSR are still true, though irrelevant in the case of distinguishable quantum 
particles A and B. 

Another inequality that can be established is 

{A{xA-XBf) + {A{pA+PB)^)>‘2 (178) 

which could also be used to establish a variance test for entangled states of the 
A and B sub-systems, so that if 

{A{xa-xb)‘^) + {A{pa+pb)^) <“2 (179) 

then the sub-systems are entangled. 


6.1.2 Two Mode Systems of Identical Bosons 


However, we can also consider cases of systems of identical bosons with two 
modes A, B rather than two distinguishable quantum particles A and B. In 
this case both the sub-systems may involve arbitrary numbers of particles, so it 
is of interest to see what implications follow from the physical sub-system states 
Pb and p^ now being required to satisfy the local particle number SSR, and all 
quantum states p satisfying the global particle number SSR. It is well-known 
that in two mode boson systems quadrature operators can be defined via 


XA 


xb 


^(a+at) fo = ^(a-at) 


(180) 


which have the same commutation rules as the position and momentum oper¬ 
ators for distinguishable particles. Thus \xa,Pa] = ^b,Pb\ = * as for cases 
where A, B were distinguishable particles. 

Since the proof of Eq. (1761 in |28| did not involve invoking the SSR, then if 
the inequality in Eq.( 1771 is satisfied, then the state would be an entangled state 
for modes A, B. as well as for distinguishable particles A, B. The situation 
would then be similar to that for the Hillery et al [23], |32| tests - the SSR 
compliant sub-system states are just a particular case of the set of all sub¬ 
system states. However, in regard to spin squeezing and correlation tests for 
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entanglement, new tests were found when the SSR were explicitly considered 
and it is possible that this could occur here. This turns out not to be the case. 


As we will see, the inequality (1761 is replaced by an equation that is satisfied 


by all quantum states for two mode systems of identical bosons where the global 
particle number SSR applies. This equation is the same irrespective of whether 
the state is separable or entangled. To see this we evaluate {A(xa + xb)‘^) + 
{A{pa —Pb)"^) for states that are global SSR compliant. 

Firstly, 

{{xa+Xb)) = {{pa-Pb))=Q (181) 

since (a) = = 0 for SSR compliant states. 

Secondly, 

1 ( (aa^) + (ab^\ + tha^\ + (W\ 

Vi) 

using (a^) = ((at)^) = = (at^ = = 0 for global SSR 

compliant states. Hence using the commutation rules , introducing the number 
operator N and the spin operator Sx and using (181) we find that 


{A{xA + XBf) = {{xa + Xb^) 

= 1 + a) + (b^ b'^ + (b^ b 

= 1 + (n)+2(s, 


Similarly 

{A{pa - PB?) = 1 + (n)-2 (Sx) 

Thus we have for all globally SSR compliant states 

{A{xa + xb^) + {A{pa - PBf) =2 + 2 (n^ 


(182) 

(183) 

(184) 


Since ( N) > 0 for all quantum states we see that the Duan et al inequality 


(176) for separable states is still satisfied, but because (184) applies for all 


states irrespective of whether or not they are separable, we see that there is no 
quadrature variance entanglement test of the form 


( A ( a ;.4 + xb)"^} + {A{pa — Pb)'^} <2 + 2 


(185) 


for the case of two mode systems of identical massive bosons. The situa tion is 
similar to the non-existent test ( AS » ) + ( AS ?, > < I (<5 z ) I in Section 


The situation contrasts that in Section 


4.3 


4.1.3 


where a test ^A(S'^ + S”^)^ 
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(^A{Sy — Sy)^j < I I establishes entanglement between two sub-systems 

(1 and 2) - but in this case each consisting of two modes. 

We can also show for all globally SSR compliant states that 


{A{xa - xbY) 


1 + (n)-2 (s,) 


(186) 


{A{pa+PbY) 


1 + (n^+2(s,) 


(187) 


and hence 


{A{xa - xbT) + {A{pa + PBf) = 2 + 2 


(188) 


but again no entanglement test results. 


The universal result (1841 for the quadrature variance sum may seem para¬ 


doxical in view of the operators {xa + xb) and {pA — Pb) commuting. Math¬ 
ematically, this would imply that they would then have a complete set of si¬ 
multaneous eigenvectors \Xa,b, Pa,b) such that {xa + xb) \Xa,b, Pa,b) = 
^A.B \Xa,b, Pa,b) and {pA — Pb) \Xa,b, Pa,b) = Pa,b \Xa,b, Pa,b)- For 


these eigenstates (A{xa + xb)‘^) = {A{pa —Pb)"^) =0 which contradicts (1841 
for such states. However, no such eigenstates exist that are globally SSR 
compliant. For SSR compliant states \Xa,b, Pa,b) must be an eigenstate of 
N and for eigenvalue N we see that {xa + xb) \Xa,b, Pa,b)ai is a linear 
combination of eigenstates of N with eigenvalues ± 1. Hence (Ja Pxb) 
\Xa,b, Pa,b) 7 ^ Xa,b \Xa,b, Pa,b )so simultaneous eigenstates that are 
SSR compliant do not exist and there is no paradox. As pointed out above, 
this issue does not arise for the case of two distinguishable particles where the 
operators xa,xb,Pa and pb are not related to mode annihilation and creation 
operators - as in the present case. 

We can also derive inequalities for separable states involving xa,Pa and 


xb,Pb based on the approach in Section 4.3 Starting with Eq. (127) we choose 


flA = X A, ^B = X B, Aa = PA an d Ag = Pb- Here 0 a = 1a and Q\b = 1b 
F or separable states we have from 


{A{axA + I3 xb)‘^) + {A{apA - Ppb)'^) > 


(189) 


With the choice of = /?^ = 1 we then find the following inequalities for 
separable states 


{A{xa + Xb 

)^) + {A{pa-Pb)^) 

> 

2 


{A{xa - Xb 

)2) -t (A(pa +Pb)^) 

> 

2 


{A{xa + Xb 

)^) + {A.{pa +Pb)‘^) 

> 

2 


{A{xa - Xb 

f) + {AipA-pB)^) 

> 

2 

(190) 


depending on the choice of a and /3. With a = /3 = 1 the first result is obtained 
and is the same as in (176). This result is consistent with (184). However using 
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(182), (187), (186) and (183) we have for global SSR compliant states - separable 
and non-separable 


(A(x^ + xbY) + (A(pa - PbY) 
{A{xa - XBf) + {A{pa + PbY) 
{A{xA + XBf‘) + {A{pA+PBf) 
{A(xa - xb)"^) + {A{pa - PbY) 


2 + 2^iV^ 

2 + 2(^iV^ 

2 + 2 + 4 

2 + 2(^n)-4(s,) (191) 


The implications for the first two equalities have been discussed above. In the 
case of the (+, +) and (—, —) cases, we note that for states with eigenvalue N for 


N the eigenvalues for Sx lie in the range —N/2 to +N/2 and hence ±2 (^Sx^ 
is always > 0. Thus (1901 will apply for both separable and entangled states. 


Hence for global SSR compliant states none of 
test. 


(190) lead to an entanglement 


6.1.3 Non SSR Compliant States 


On the other hand if neither the sub-system nor the overall system states are 
required to be SSR compliant - though they may be - we find that for separable 
states 

{A(xA±XBf)p+{A{pATPBf)p > 2+2 (^n'^^+ 2{(h'b^ + (a'^V'^)-2\ i^)p+(b'’) 

(192) 

so entanglement based on ignoring local particle number SSR in the separable 
states is now shown if 


{A(xa ± xb)'^)p+{A{pa + PbT) p < 2+2 {^Nj + 2 ( 


iS) )-2|(S),+(?) 


' p 

(193) 

However, even if local particle number SSR compliance is ignored for the sub¬ 
system states (as in Ref m). global particle number SSR compliance is still 
required for the overall quantum state. This applies to both the separable states 
and to states that are being tested for entanglement. In this case the quantities 

are all zero, so the entanglement 

test in (193) would become the same as the hypothetical entanglement test (185). 

For the sceptic (see Section [3.1.6 ) who wishes to completely disregard the 
SSR (both locally and globally) and proposes to use tests based on quadrature 
variances such as (193) to establish entanglement, the challenge will be to find 


a way of measuring the allegedly non-zero quantities .. {^1 ■ This 

would require some sort of system with a well-defined phase reference. Such a 
measurement is not possible with the beam splitter interferometer discussed in 
this paper, and the lack of such a detector system would preclude establishing 
SSR neglected entanglement for systems of identical bosons. Essentially the 
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same problem arises in testing whether states that are non-SSR compliant exist 
in single mode systems of massive bosons. 

As mentioned previously, the result in Eq. (176) was established in Ref. [28] 


without requiring the sub-system states to be compliant with the lo¬ 

cal particle number SSR or the density operator p for the state being tested 
to comply with the global particle number SSR, as would be the case for 
physical sub-system and system states of identical bosons. However, in Ref. 
pS] it was pointed out that two mode squeezed vacuum states of the form 
1$) = exp(—— ah)) |0) satisfy the entanglement test. However, such 
stand alone two-mode states are not allowed quantum states for massive iden¬ 
tical boson systems, as they are not compliant with the global particle number 
SSR. To create states with correlated pairs of bosons in modes a and b processes 
such as the dissociation of a bosonic molecular BEC in a mode M into pair of 
bosonic atoms in modes a and b can indeed occur, but would involve interaction 
Hamiltonians such as H = K(a^b'^M + M^ab). The state produced would be an 
entangled state of the atoms plus molecules which would be compliant with the 
global total boson number SSR - taking into account the boson particle content 
of the molecule via N = 2nM + na + Ub- It would not be a state of the form 
1$) = exp(—~ |0)- 


6.2 Reid 1989 

Another test involves the general quadrature operators defined as in 


], for 


which those in (180) are special cases 


= -^(aexp{—i6)+a^ exp{+i9)) 

v2 

= -^{bexp{-i(l)) + b^ exp{+i(j))) 

v2 


(194) 


These operators are Hermitian. The conjugate operators are 

Pa = exp{-i9) - exp(-ht6»)) = X^+^/^ 

Pt = ^(^exp(-i0)-Pexp(-hi(/))) = (195) 

y li 

where [XlPt] = [Xt,Pt]=^■ 

Noting that for any state we have (^{X^ — XXf)'^'^ > 0 for all real A estab¬ 
lishes the Cauchy inequality for all quantum states 


'^ab — 




{x‘r‘){{x*y 


< 1 


(196) 


The quantity is a correlation coefficient. For SSR compliant separable states 




EH 


whilst for all globally SSR compliant 
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states = (”a) + 5 > 5 ((^^)^) = i^b) + 5 > j- Hence for 

SSR compliant separable states the correlation coefficient is zero. A quadrature 
correlation test for entanglement based on locally SSR compliant sub-system 
states is then 


CJ ^ 0 


(197) 


However, it is not difficult to show that for states that are globally SSR com¬ 
pliant 


{Kxt) = 


Sx ) cos(6* — (j)) + (Sy \ sin(0 — (j>) 


(198) 


so that the entanglement test based on locally SSR compliant sub-system states 
is equivalent to finding one of (Sx) or ( Sy) to be non-zero. This is the same 


as the previous Bloch vector test in Eq. (681 or the weak correlation test in 


Eq.(166|. 
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6.3 Two Mode Quadrature Squeezing 


From Eq. (194) we can define two mode quadrature operators as 

Xe(+) = + X®) = i(aexp(-i0) + 6'^ exp(+i0) + exp(+i0) + &exp(-i0)) 

V 2 2 

Ps{+) = -^{Pa + Pb) = ;^(aexp(-i0) - 6'^exp(+i0) - a'^ exp(+i0) + 6exp(-i0)) 

V2 


— Xe+TT/2{+) 


where we have [ Xg{+),Pg{+)] = i. Note that ^o(+) = {xa + xb)/V ^ and 

As we 


(199) 


6 . 1.2 


^o(+) = {pa +Pb)IV^ unlike the operators considered in Section 
have seen there is no entanglement test for systems of identical bosons of the 
form (A(a;^ + + —PbY) < 2+2 The Heisenberg uncertainty 

principle gives ^AXg(+)^ ^APg (+)^ > 1/4, so a state is squeezed in Xg{-\-) if 

AAg(+)\ < 1/2, and similarly for squeezing in Pg{+). 


We can show that for separable states both ( AXg (+) / > 1/2 and ( APg (+)) > 


1/2, so two mode quadrature squeezing in either Xg{-\-) or Pe{+) is a test for 
two mode entanglement. Firstly, for SSR compliant sub-system states 




(-??)„) = 0 


( 200 ) 


since (a) ^ = (b) =0. Secondly, 


1 


R 


> - 


R 

\ + \iN 
1 
2 


where for local SSR compliant states other terms involving (a^)^ , 
(a)fl > ^^^)r = ^^)r . are all zero. Hence 


( 201 ) 

a&\ = 

/ R 


^Xl{+)) = IXI{+)) - (Xg{+)) >- 


( 202 ) 
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which establishes the result. Since Pe{+) = -^e+ 7 i-/ 2 (+) we also have (^APg{+)J = 

5 5 ^ separable state. Hence the two mode quadrature squeezing 

test. If 


(ax 2(+)) <1 or {APi(+)) < i (203) 


then the state is entangled. Obviously Xg{+) and Pd{+) cannot both be 
squeezed for the same state. 

We can also define additional two mode quadrature operators as 


M-) 


Psi-) 


1 

71 

1 

71 

Xg+n/zi-) 


iK - 

iP!-Pb) 


-(aexp(—i0) — 6'^exp(+i0) +d’^exp(+i0) — &exp(—i0)) 

7(aexp(— 10 ) + exp(+i0) — d’^ exp(+i0) — 5exp(—i0)) 
2i 


(204) 


where we also have [ Xg{—),Pe{—)] = i. Again (^AXg{—)'j (^APg(—)'j > 1/4, 
so a state is squeezed in Xg{—) if (^AXg{—)'j < 1/2, and similarly for squeezing 
in Pe{-). 

A similar proof shows that for separable states both ^AA|(—)^ ~ k P 
i > 1/2 and ^AP2(_)^ = I + i (7) > 1/2, so if 

(aa 2(-)) <1 or (aP^H) < I (205) 


then the state is entangled. Hence any one of Xg{+),Pg{+), Xg(—), Pe(—) being 
squeezed will demonstate two mode entanglement. 

The question then arises - Can two of the four two mode quadrature op¬ 
erators be squeezed? For simplicity we only discuss 0 = 0 cases in detail. 
Obviously pairs such as Ao(+))7)(+) or Ao(—),Po(—) cannot. Next, we con¬ 
sider Ao(-I-) and Po(—)■ We note that for all global SSR compliant states 
(aAq (+)^ + (^^Po(.~)') = 5 ((^(^A + xb)"^) + {A{pa - Pb)"^)) = 1 + {n ^ 

using (184), so that if Ao(-I-) is squeezed ('AAg(-l-)') < 1 then ^APg(—)*) > 

1 
2 

spite of the operators commuting. A similar conclusion applies to Ao(—) and 
Po(-I-). For the pair Ao(-I-) and A‘o(—) we have ^AAq (+)^ -|- ^AAg(—)^ = 

1 ((A(a;^ + xbY) + {A{xa — xb)'^)) = 1 + (^n'^ using (1821 and (186), so the 

same situation as for Ao(+) and Po(—) applies, and thus Ao(-I-) and Ao(—) 
cannot both be squeezed. A similar conclusion applies to Po(—) and Po(-I-). In 
general, only one of Xg{+), Pe{+), Xe{—), Pg{—) can be squeezed. 


^A), showing that both Ao(-I-) and Po(—) cannot both be squeezed - in 
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Further questions are: What quantities need to be measured in order to test 
whether two mode quadrature squeezing occurs and how useful would it be to 
detect entanglement? It is straight-forward to show from (199) and (204) that 
for states that are global SSR compliant 


(^e(+)) 

= 

0 {Xei-)) 

= 0 


(206) 

(ax|(+)) 

= 

(i?(+)) = i ( 

(S) 

+ 2 


(ax|(-)) 

= 

(«'-)) - 5 ( 

(«) 

- 2 ( 5 ,)) 

(207) 

since terms such as (a^), 

(oS) 

etc are all zero 

for SSR compliant 

states. As 

explained in Section 

12.2 

both 

(n) + 2 (Sx) 

and 

(n) - 2 (Sx) 

are always 


non-negative, but the entanglement test would require 

1 for squeezing in 

1 for squeezing in Xg{—) 


{n) + 2( S, 
(fv) - 2 ( S, 


< 

< 


(208) 


This shows that the two mode quadrature squeezing test involves measuring 
and so that once again measurements of boson number and the 

mean value of a spin operator are involved. Similar conclusions apply for Pe{+) 
and Pg{—). However, since the test requires (Sx) to be non-zero it would 


simpler to use the Bloch vector test (see (68)) which merely requires showing 
that one of ( S'a,) or { Sy) to be non-zero. 


In most cases the inequalities in (|208|) will not be satisfied, since both (Nj 


and (Sx) are 0{N). However, for the binomial state in (77) with 9 = 37r/4 


and y = 0 we have for |<i)) = |0) /VN\ the results = X and 

= —N/2 (see (163) in Ref. [6]). Hence spin squeezing in Xe{-\-) occurs, 
confirming that this particular binomial state is entangled. Note that the test 
does not confirm entanglement for almost all other binomial states (those where 

Sx) is different from ±A^/2), though these are actually entangled. 
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7 Interferometry in Bosonic Systems 

In this section we discuss how interferometers in two mode bosonic systems 
operate. This topic has of course been discussed many times before, but for 
completeness we present it here. Our approach is essentially the same as in ear¬ 
lier papers, for example that of Yurke et al [3S]. Before discussing interferometry 
in two mode bosonic systems, we first need to set out the general Hamiltonian 
for the two mode systems that could be of interest. The two modes may be 
associated with two distinct single boson spatial states, such as in a double well 
potential in which case the coupling between the two modes is associated with 
quantum tunneling. Or they may be associated with two different atomic inter¬ 
nal hyperfine states in a single well, which may be coupled via classical fields in 
the form of very short pulses, for which the time dependent amplitude iSvA(t), 
the centre frequency is luq and the phase is f. Since this coupling process is 
much easier to control than quantum tunneling, we will mainly focus on the 
case of two modes associated with different hyperfine states, though the ap¬ 
proach might also be applied to the case of two spatial modes. The free atoms 
occupying the two modes are associated with energies hiUa, hujb, the transition 
frequency ujba = uib — to a being close to resonance with loq- It is envisaged that 
a large number N of bosonic atoms occupy the two modes. The bosonic atoms 
may also interact with each other via spin conserving, zero range interatomic 
potentials. We will show that measurements on the mean and variance for the 
population difference determine the mean values and covariance matrix for the 
spin operators involved in entanglement tests. 

For interferometry involving multi-mode systems, a straightforward general¬ 
isation of the two mode case is possible, based on the reasonable assumption the 
interferometer process couples the modes in a pairwise manner. This is based 
on the operation of selection rules, as will be explained below. 

However, although in the present section we show that two mode interfer¬ 
ometers can be used to measure the mean values and covariance matrix for the 
spin operators involved in entanglement tests for systems of massive bosons, 
the issue of how to measure mean values and variances for the quadrature 
operators involved in other entanglement tests for massive bosons is still to be 
established. Such a measurement is not possible with the beam splitter in¬ 
terferometer discussed in this paper. An interferometer involving some sort 
of phase reference would seem to be needed. Proposals based on homodyne 
measurements have been made by Olsen et al [37], [38], but these are based on 
hypothetical reference systems with large numbers of bosons in Glauber coherent 
states, and such states are prohibited via the global particle number SSR. 

7.1 Simple Two Mode Interferometer 

A simple description of the two mode system is provided by the Josephson 
model, where the overall Hamiltonian is of the form [^ 

H. Joseph = Hq I4o1 (209) 
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with 


Hq = huja a^a + hwb b^b 

V = A(t) exp{—iu!ot) e:xjp{i(j)) b^a + A{t) eKjp{+iu!ot) exjp{—i(f>)a^b 

Vcoi = xib'^b-a'^af ( 210 ) 

where Hq is the free boson Hamiltonian, V gives the interaction with the clas¬ 
sical field and 14oi is the collisional interaction term. For the case of quantum 
tunneling between two distinct spatial modes, the interaction term V can also 
be described in the Josephson model (see [B] for details), in which case the fac¬ 
tors multiplying b^a or a^b involve the trapping potential and the two spatial 
mode functions. A time dependent amplitude and phase might be obtained via 
adding a suitable time dependent field to the trapping potential - this would be 
experimentally difficult. The Hamiltonian can also be written in terms of spin 
operators as 

Ho = l/2{hlUa + huJb)N - huJabSz 

V = A{t) exp(—iwot) exp{i(j)) {S^ + iSy) + H.C 

Koi = 4x^2 (211) 

The coupling effect in a simple two mode interferometer can be described via 
the classical interaction term V, where now the amplitude A(t) is only non-zero 
over a short time interval. The pulsed classical field is associated with an area 
variable s, defined by 

t 

s = y dtiA{ti)/h (212) 

*0 

the integral eventually being over the pulse’s duration. The interferometer fre¬ 
quency loq is assumed for simplicity to be in resonance with the transition 
frequency ujba = ^b — ^a- The classical field is also associated with a phase 
(f, so the simple two mode interferometer is described by two interferometric 
variables 2s = 9 giving the pulse area and (j) specifying the phase. Changing 
these variables leads to a range of differing applications of the interferometer. 
When acting as a beam splitter (BS) a 2s = 7r/2 pulse is involved and (j) is 
variable, but for a phase changer a 2s = tt pulse is involved is arbitrary). 
For state tomography in the yz plane we choose 2s = 9 (variable) and ^ = 0 or 
TT. The beam splitter enables state tomography in the xy plane to be carried 
out. Generally speaking the effect of collisions can be ignored during the short 
classical pulse and we will do so here. 

7.2 General Two Mode Interferometers 

More complex two mode bosonic interferometers applied to a specific input 
quantum state will involve specific arrangements of simple two mode interferom¬ 
eters such as beam splitters, phase changers and free evolution intervals, followed 
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by final measurement of the mean population diference between the modes and 
its variance. Ramsey interferometry involves two beam splitters separated by a 
controllable free evolution time interval T. During such an interval in which free 
evolution occurs, the interaction of the classical beam splitter field with the two 
mode system can be ignored, but the effect of collisions and coupling to external 
systems may be important if collision parameters are to be measured using the 
interferometer. The overall behaviour of such multi-element interferometers will 
also depend on the initial two mode quantum input state that acts as the input 
state for the interferometer, as well as important variables such as the phase (j>, 
the centre frequency ojq, the area variable s for the classical pulses used, and 
also the the free evolution intervals (if any). The behaviour also will depend 
on the characteristic parameters such as the transition frequency Wab, collision 
parameter y and total boson number N for the two mode system used in the 
interferometry. The variables that describe the interaction with other systems 
whose properties are to be measured using the interferometer must also affect its 
behaviour if the interferometer is to be useful. Finally, a choice must be made 
for the interferometer physical quantity whose mean value and quantum fluctu¬ 
ation is to be measured - referred to as the measurable. The outcome of such 
measurements can be studied as a function of one or more of the variables on 
which the interferometer behaviour depends - referred to as the interferometric 
variable. There are obviously a wide range of possible two mode interferometers 
types that could be studied, depending on the application envisaged. Interfer¬ 
ometers also have a wide range of uses, including determining the properties of 
the input two mode state - such as squeezing or entanglement. For a suitable 
known input state they can be used to measure interferometric variables - such 
as the classical phase f of the pulsed field associated with a beam splitter or 
a parameter associated with an external system coupled to the interferometer. 
On the other hand, in a Ramsey interferometer the interferometric could be the 
collision parameter y, obtainable if the free evolution period T is known. No 
attempt to be comprehensive will be made here. 

The Ramsey interferometer is illustrated in Figure 7. 
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Figure 7. Ramsey Interferometry. Two 7r/2 beam splitters separated by a free 

evolution region. 
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For the purpose of considering entanglement tests a simple two mode in¬ 
terferometer operating under conditions of exact resonance loq = ooab will be 
treated, and its behaviour for N large when the phase </> is changed and for differ¬ 
ent choices of the input state p will be examined. Measurements appropriate to 
detecting entanglement via spin squeezing and correlation will be discussed. The 
measurable chosen will initially be half the population difference {lf^b — d)d)/2 - 
which equals Sz - generally measured after the two mode system has interacted 
with the simple interferometer, but also without this interaction. The phase </> 
will act as the interferometric variable, as will the pulse area 2s = 9. As we will 
see, different choices of input state ranging from separable to entangled states 
lead to markedly different behaviours. In particular, the behaviour of relative 
phase eigenstates as input states will be examined. Later we will also consider 
measurements involving the square of Sz- 

7.3 Measurements in Simple Two Mode Interferometer 

As discussed in the previous paragraph, the initial choice of measurable is 

M =^{M-d'<d) = Sz (213) 

and we will determine its mean and variance for the state given by 

'p* = UpU-^ (214) 

where the output state p"^ has evolved from the initial input state p due to the 
effect of the simple two mode interferometer. U is the unitary evolution operator 
describing evolution during the time the short classical pulse is applied. Collision 
terms and interactions with other systems will be ignored during the short time 
interval involved. 

We note that for an N boson state the eigenvalues of M range from —N/2 
to -|-iV/2 in integer steps. For more general states the possible values for M 
are any integer or half integer. When M is measured the result will be one of 
these eigenvalues, but the average of repeated measurements will be which 
must also lie in the range —N/2 to -|-A^/2. The variance of the results for the 
repeated measurements of M is also experimentally measureable and will not 
exceed (iV/2)^, and apart from NOON states will be much less than this. The 
experimentally determinable results for both and will depend on 

p and on the interferometer variables such as the phase </> and the pulse area 
2s = 9. 

The Hamiltonian governing the evolution in the simple two mode interfer¬ 
ometer will be Hq V. For the output state the mean value and variance are 

(m'^ = TrlM'p*) 

= Tr{^M - (m'jY' p*) (215) 
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These will be evaluated at the end of the pulse. If the input state is measured 
directly without applying the interferometer, then the mean value and variance 
are as in the last equations but with 'p^ replaced by p. 

The derivation of th e results is set out in Appendix and are given by 
the same form as (215), but with replaced by p and with M replaced by 


the interaction picture Heisenberg operator at the end of the pulse, 

which is given by 


Mh{2s,4>) = i(6|^(s,(/>)6i/(s,(/))(?())) 

= sin2s (sin(/)S'a; + cosc^S'y) + cos2s^^ 


(216) 


with 


6/r(s, 4)) = cos sb — i exp(j^) sin s a anis, (j)) = exp(—sin sb + cos so 

(217) 

The versatility of the measurement follows from the range of possible choices of 
the pulse area 2s = 9 and the phase (j). These results are valid for both bosonic 
and fermionic modes. 

We then Hnd that the general result for the mean value is 


M ) = sin 0 sin (() ( S'a;) + sin 9 cos (j) \^Sy^ + cos 9 (^Sz 


and for the variance is 

(1 — cos 29) (1 — cos 2(j)) 


(218) 


2 2 
(1 + cos 29) 


mo (1-C0S26I) (l + cos2(;i)^^- 

O -r ^ ^ Lj\Dy^Dy} 


C{Sz,Sz) 


(1 — cos 20) 


sin20(7(53;, Sy) + sin20 cos0(7(5^, 5^) + sin20 sin0(7(52, S^f) 


(219) 


where the mean value of the spin operators are = Tr{Sa'p) and the 

covariance matrix elements are given by C{SaiSp) = 1/2 ({SaSp + SpSa)J ~ 

■ The diagonal elements (7(5 q, 5a) = ^5 q^^ — l^Sa'^ = ^A5a^ 

is the variance. By making appropriate choises of the interferometer variables 0 
(half the the pulse area) and 0 (the phase) the mean values of all the spin oper¬ 
ators and all elements of the covariance matrrix can be measured. Tomography 
for the spin operators in any selected plane can be carried out. 
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7.3.1 Tomography in xy Plane - Beam Splitter 

In the beam splitter case (for state tomography in the xy plane) we choose 
‘Is = 'Kl‘l and (p (variable) giving 

Sx + coscj) Sy ( 220 ) 

and we find that for the output state of the BS interferometer the mean value 
and variance of M are given by 

siiKj) + cos (j) (^Sy'^ (221) 

i(l - cos2(^) C{S^,Sj;) + ^(1 + cos2(^) C{Sy,Sy) + sin2(() C{Sa:,Sy) 

( 222 ) 

showing the mean value for the measureable M depends sinusoidally on the 
phase (j) and the mean values of the spin operators Sx,Sy. The variance for 
the measurable depends sinusoidally on 2(j) and on the covariance matrix of the 
same spin operators. Both the means and covariances of the spin operators 
Sx, Sy now depend on the input state p for the interferometer rather than the 
output state 

7.3.2 Tomography in yz Plane 

For state tomography in the yz plane we obtain the means and covariances of 
the spin operators SyfSz- To do this we choose 2s = 9 (variable) and ^ = 0 so 
that 

MHi9,0) = sine Sy +cose Sz (223) 

and find that for the output state the mean value and variance of M are given 

by 

sin9 (^Sy'j +cos9(^Sz'j (224) 

i(l - cos26») C{Sy, 5y) + ^(1 + cos 20) C{Sz,Sz) + sin 20 C{Sy, Sz) 

(225) 

A single measurement does not of course determine the mean value 
An average over a large number of independent repetitions of the measurement is 
needed to accurately determine which can then be compared to theoretical 
predictions. This is a well-known practical issue for the experimenter that we 
need not dwell on here. A brief account of the issues involved is included in 
Appendix 
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( 226 ) 


7.3.3 Phase Changer 

In this case we choose 2s = 9 = tt and (p (arbitrary) giving 

= -Sz 

and for the output state the mean value and variance of M are given by 

(m) = - (227) 

(228) 

so the phase changer measures the negative of the population difference. Effec¬ 
tively the phase changer interchanges the modes a —>■ 5 and 6 —>■ a and this is 
its role rather than being directly involved in a measurement. Phase changers 
are often included in complex interferometers at the midpoint of free evolution 
regions to cancel out unwanted causes of phase change. 


7.3.4 Other Measurements in Simple Two Mode Interferometer 

Another useful choice of measureable is the square of the population difference 


M2= Q(St6-ata)^ 


(229) 


For the beam splitter case with 2s = 7r/2 and p (variable), we can easily show 
(see Appendix 171 that the mean value of M 2 for the output state is given by 

M 2 ^ =sin^(() ((SxY'^ +cos^ (j) ((SyY'^ + sin (j) cos (j) ((S^Sy + SyS^)^ (230) 

showing that the mean for the new observable M 2 is a sinusoidal function of the 
BS interferometer variable (p with coefficients that depend on the means of 5"^, 
Si and Sjy + SyS^. 

Choosing special cases for the interferometer variable yields the following 
useful results 


M; 




M2) ={{S,Y 

cj)—'Kj2 \ Ip 


(M 2 ) -(M 2 

\ / (f)—Tz{A \ 


= aS^Sy + SyS,)] 

b=-7r/4 \ /p 


(231) 


Hence all three quantities (( [Sy)^ ) and ( (SxSy -I- SySx )) can be mea¬ 


sured. We note that just measuring (M 2 j does not add to the results ob¬ 
tained by measuring the mean and variance of the original measureable M, 
since (mY^ = (iS.M'^') + variance ^AM|^ does of course de¬ 
pend on higher moments, for example with d) = 0 (AM^) = [A \ and 
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7.4 Multi-Mode Interferometers 


For the multi-mode case we consider two sets of modes di and hi as descibed 
in Subsection |2.2| These may be different modes associated with two different 
hyperfine states or they may be modes associated with two separated potential 
wells. The Hamiltonian analogous to that in (210) for the two mode case is 


Hq = h{uja + uji)a\ai + h{u;b + 

i i 

V = A{t) exp{—iuJot) exp{i(j>)'^^bla.i + A{t) exp{+iuJot) exp{—i(p) 

i i 

(232) 


where the collision terms are ignored since we are only considering the effect of 
the short interferometer coupling pulse. Here we have assumed that the energy 
for the mode d^ is h{u}a+uji), which is the sum of a basic energy for all a modes - 
hu!a, and an energy term huji that distinguishes differing modes di (and similarly 
for the mode bi). In addition, we assume selection rules lead to pairwise coupling 
at bi. In the case where coupling is due to pulsed external fields (microwave 
and RF) we can assume that the momenta (~ mhujtrap) associated with 
trapped modes d^ and hi are the same, since the momenta associated with the 
low frequency photons (~ hujRF jc) involved can be ignored. The spin operators 
for the multi-mode system are set out in (|^ in terms of the mode operators. 

As in SubSection |7.1| the interferometer frequency wq is assumed for simplic¬ 
ity to be in resonance with the transition frequency ujba = aib — oJa- Following 
the treatment in SubSection 7.3 the choice of measurable is the half the total 


population difference between the two sets of modes 

1 


M=- \^(b]h - aJdO = S., 


(233) 


and we will determine its mean and variance for the state given by 

p# = [/pt7-i (234) 

where the output state has evolved from the initial input state p due to 
the effect of the multi- mode interferometer. U is the unitary evolution operator 
describing evolution during the time the short classical pulse is applied. Collision 
terms and interactions with other systems will be ignored during the short time 
interval involved. 

^s in the two mode interferometer case, the results for the mean and variance 
of M depend on the pulse area 2s = 9 and the phase (f>.of the interferometer 
coupling pulse. They have the same dependence on the mean values and co- 
variance matrix for the multi-mode spin operators S^, Sy and for the input 
state p as in the two mode interferometer. Thus we then find that the general 
result for the mean value is 

I^m'^ = sm.9 sin<f sin9 cos 4> (^Sy'^ -\-cos9(^S^ (235) 
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and for the variance is 


AM^ 

(1 — cos 26) (1 — cos 2(p) 


r<to (1-COS26I) (l + cos2(;i)^^j; ^ 

O -r Lj\Dy^Dy) 


(1 + cos 26) 


C{S,,S,) 


sin2(/)C'(S'a;, Sy) + sin20 cos(l)C{Sy, Sz) + sin20 sin0C'(5'2, Sx) 

(236) 


Details of the derivation are set out in Appendix [T7| 


7.5 Application to Spin Squeezing Tests for Entanglement 


Unless stated otherwise, we now focus on spin squeezing tests for SSR compliant 
entanglement based on the beam splitter measurements (the simple two mode 
interferometer with 2s = 0 = 7r/2). By choosing the phase (j) = 0 we see 

that = (^Sy'^ and = C{Sy,Sy) = (^{Sy — (^Sy'^ giving the 

mean and variance for the spin operator Sy. By choosing the phase (j) = n 12 we 
see that and = C{Sx,Sx) = (^{Sx - (Sx'^ giving 

the mean and variance for the spin operator Sx ■ If the measurement of ( M \ 


were carried out without the beam splitter being present then ■ 

Combining all these results then enables us to see whether or not Sx is squeezed 
with respect to Sy or vice versa. This illustrates the use of the interferometer 
in seeing if a state p is squeezed. Squeezing in Sz with respect to Sy (or Sx) 
or vice versa also demonstrates entanglement and again the simple two mode 
interferometer with appropriate choices of 0 and (j) can be used to measure the 
means and variances of the relevant spin operators. 

As the presence of spin squeezing shows that the state must be entangled 
[T] the use of the interferometer for squeezing tests is important. Note that we 
still need to consider whether fluctuations due to a finite measurement sample 
could mask the test. However, as spin squeezing has been demonstrated in two 
mode systems of bosonic atoms this approach to demonstrating entanglement 
is clearly useful. 
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7.5.1 Spin Squeezing in Sx, Sy 

To demonstrate spin squeezing in Sx with respect to Sy we need to measure the 
variances (^ASx"^^ and (^ASy"^^ and the mean (^S^ 


and show that 




<2' 

P ^ 




(237) 


As we have seen, the variances in Sy^ Sx are obtained by measuring the fluctua¬ 
tion in the measureable M after applying the interferometer to the state p, with 
the interferometer phase set to <() = 0 or ()) = 7r/2 for the two cases respectively. 
The mean is obtained by a direct measurement of the measureable M 

without applying the interferometer to the state p. 


7.5.2 Spin Squeezing in xy Plane 

As shown in Section (see [1] ) squeezing in Sx compared to Sy or vice versa 
is sufficient to show that the state is entangled. However, as the interferometer 
measures the variance for the state p in the quantity 


Stt 


M//(7r/2, ())) = sm(j)Sx + cos(pSy = Sf{— + ip) 


(238) 


corresponding to the x component of the spin vector operator (ig^ after it has 
been rotated about the z axis through an angle ^ it is desirable to extend 
the entanglement test to consider the squeezing of Sf (^ + 0) with respect to 
the corresponding y component Sf (f>) - and vice versa, where 


-S'f (-^ +<(') = -cos4>Sx +sm(j)Sy 


(239) 


The variance in S^{^ + </>) can be obtained by changing the interferometer 
phase to ()) + f. Clearly [Sf{^ + <()), Sf{^ +<(')] = * >5'z, as before. 

The question is - does squeezing in either Sf (^ + <()) or S^+4>) demon¬ 
strate entanglement of the modes a and b ? The answer is that it does. 

For a separable state we have 

+ + </>)) =0 (240) 


as before, since ( Sft (^ +<(')) are just linear combinations of the zero ( Sx 

^ ^ ' P ^ ' ' p 

Since [ Sfiftf- -\- cj)), S^{^ + (jfl] = iS^ the Heisenberg uncertainty principle 

shows that (^ASf{^ (^ASf{^ > || (^Sz'^ P so spin squeez¬ 

ing in sf (^ -I- (/>) with respect to Sf (j)) or vice versa requires us to show 
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that 


(241) 

Since we measure (^AS^ (^ + 0)^ ^ the spin squeezing test is (A5#(f+ 0)2) < 


^1(5. 


Now for Sf{^ + (j)) we have 

As*(^+0f) = (stif+4'r) 


p \ Ip 

2 ^ / C 2 \ I „„„2 


= sin 0 (<52 ) + cos 0 ( S'y ) + sin 0 cos 0 \ SxSy + SySx 


and for a separable state we have from SubSection 2.4 


52 ) = {ASl)>\\{S. 


Si) = {^ASi)>-\{S. 


whilst for the remaining term 

1 
p 


s.s, + §,s.) = 1({(S')"(S)*-+)"(+}) 




^2\ 


(242) 


(243) 


= 0 


(244) 


using the local particle number SSR. 

Thus as sin^ 0+cos2 0=1 and applying similar considerations to + 0)^ 

(a 4 *(|++)^ > i|(a); 

^A§#(^+0)2^ > (245) 

showing that for a separable state there is no squeezing for S'|‘(^ + 0) compared 
to Sf (^ + 0) or vice versa. Hence squeezing in either Sf {^ + (j)) or Sf {^ + (j)) 
demonstrates entanglement of the modes a and b. 


7.5.3 Measurement of ( S 


The question remaining is whether the mean value \^Sz) can be measured ac¬ 
curately enough to apply the test for entanglement. With an infinite number 
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of repeated measurements this is always possible, since then both the variances 
ASf{^ + 4>Y) and the mean (Sz) ■ would become error free. For a finite 


p ' ' p 

number of measurements R the measurement of requires a considera¬ 

tion of the variance in Sz. For general entangled states general considerations 
indicate that this mean will be of order N and the variance will be at worst 

of order N^, Hence the variance /a(^Sz^\ in determining the mean 

Sz'j for R repetitions of the measurement would be ~ N'^/R, giving a fluctu¬ 
ation of ~ N/'/R. For this to be small compared to ~ fV we merely require 
i? ^ 1, which is not unexpected. This result indicates that the application of 
the spin squeezing test via interferometric measurement of both the variances 
ASf{^ + and the mean (looks feasible. 


7.5.4 Spin Squeezing in Sz, Sy 

To demonstrate spin squeezing in Sz with respect to Sy we need to measure the 
v„i.„ces and (aS/)^ and the mean (S.) and show that 


AS. 


1 , 




(a§,“ 




(246) 


As we have seen, the variances in Sz , Sy are obtained by measuring the fluctua¬ 
tion in the measureable M after applying the interferometer to the state p, with 
the interf erom eter phase set to (/) = 0 and the pulse area 2s = 9 made variable. 
From Eq.(225) we see that choosing 9 = 0 gives (ASz^) and choosing 9 = 


gives (ASy'^) . From Eq.(221) the mean (Sx) is obtained by a measurement 


of the mean in the measureable M after applying the interferometer to the state 
p, with the interferometer phase set to (f) = it/2 and the pulse area 2s = 7r/2. 


7.6 Application to Correlation Tests for Entanglement 

7.6.1 First Order Correlation Test 

Unless stated otherwise, we again focus on correlation tests for SSR compliant 
entanglement For the beam splitter case and by choosing the phase </> = 0 
we see that = 

(m) = 


and by choosing the phase ^ = 7r/2 we see that 


Sx) ■ The simplest form of the correlation test with n = m = 1 

p 


requires 


5 .) ^0 

p 


7^0 


(247) 
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for establishing that the state is entangled. For the separable state with Mh = 
sin (j)Sx+ cos (j)Sy = Sf (^ + </>) 




(248) 


so that the mean value of the measureable is zero and independent of the beam 
splitter phase (/).for all (j). Finding any non-zero value for (^Mh'^ would then 
show that the state p is entangled. More importantly from the general result, 
Mh) would show a sinusoidal dependence on the interferometer phase </>, so 

' p 

the appearance of such a dependence would be an indication that the state was 
entangled. 

The question remaining is whether the mean values (^Sx,y ^ can be measured 

accurately enough to apply the test for entanglement. With an infinite number 
of repeated measurements this is always possible, since then both the variances 


and the means ( S'- 


. would become error free. For a finite number 


of measurements the measurement of ( S. 


requires a consideration of the 

variances in Sx,y (or Mh to cover both cases). For a general entangled state 
we can assume that I^Mh'^ ^ N/2 and the variance will be at worst of order 

Hence the variance /a(^Mh'^ \ in determining the mean (^Sx 

for R repetitions of the measurement would be ~ N'^/R, giving a fluctuation 
of ~ N/\/~R. For this to be small compared to ~ we merely require i? 1, 
which is not unexpected. This result indicates that the application of the simple 
correlation test via interferometric measurement of and looks 

feasible. 


7.6.2 Second Order Correlation Test 

For the second order form of the correlation test with n = m = 2 requires 

2 


(a5. 


Sx)'' -(as,, 


V / {^y) ^ 

p \ ! p \ Ip \ I p 


{CiSx,Sy) + {Sx)^{Syj) ^ 0 


(249) 


We have already shown using Eqs (221) and (222) that the variances (ASy 
and (aSx'^'^ and the means (iSyJ and can be obtained via the BS 

interferometer.from the mean (M) and the variance ( AM^ ) for the choices of 


(j) = 0 and (j) = -n 12. To obtain the covariance matrix element C{Sx, Sy) we see 
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that if we choose ^ = 7r/4 then ~ 5 { 




+ C{S,,Sy), 


(ASy'^i -r>-'yJ, 
P ' ' P 

from which the covariance can be measured. Thus the second order correlation 
test can be applied. 

Alternately, if the measurement quantity for the BS interferometer is the 


square of the population difference then we see from (231) that the mean value 

of M, for certain choices of the BS variable 0 measures + (s,)' = 

2 


Sx 


ASy^ 


Sy) = 




p \ / p \ / p \ / p 

The second order correlation test is that if 


and {C{S^,Sy)US, 




S. 




{SxSy + SySx)) ^ 0 


) = ( {SxSy + SySx) 


(250) 


then the state is entangled. 


7.7 Application to Quadrature Tests for Entanglement 

As we saw previously, no useful quadrature test for SSR compliant entangle¬ 
ment in two mode systems of identical bosons of the form (A(a:^ AxbY)^ + 

{A{pa T PbY) < 2 -b 2 ( A^\ results if the overall system state is globally SSR 


compliant. However, if the separable states are non-compliant then showing 
that 

(A(7a±Sb)^)p+(A(p_aTPb)^)p < 2-b2(|iV^ +2{(a'bj +(a''V'^ )-2\{a)p+(V' 

(251) 

would demonstrate entanglement. This test requires measuring together 

with , (fl)p and (h^'^ . Although the first can be measured using the BS 

interferometer the quantities j (fl)p and (b^'^ cannot. Another technique 

involving a measuring system where there is a well-defined phase reference is 
therefore required if quadrature tests for SSR neglected entanglement are to be 
undertaken. Furthermore, the overall state must still be globally SSR compli¬ 
ant, and hence (^b^ , (d)^ and ^6'^^ are all zero, even for entangled states, so 

the test reduces to (A(x^ ± Js)^) + (^A(j)aT Pb)'^) <2-|-2/a^\ . Since for 

all states {A(xa ± xbY) + {A{pa TPbY) =2 + 2 (n\ this test must fail 
anyway. We have also seen that finding the correlation coefficient - defined in 


terms of generalised quadrature operators (1941 in Eq. (196) - to be non-zero 


does not lead to a new test for SSR compliant entanglement. The tests involv¬ 
ing two mode quadrature squeezing look more promising, assuming the relevant 
quadrature variances can be measured. 
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8 Experiments on Spin Squeezing 


We now examine a number of recent experimental papers involving squeezing 
and entanglement in BEC with large numbers of identical bosons. Their no¬ 
tation will be modihed to be the same as here. There are really two questions 
to consider. One is whether squeezing has been created (and of which type). 
The second is whether or not this demonstrates entanglement of the modes in¬ 
volved. Here we define entanglement for identical bosons as set out in Section 3 
of paper I. Many of these experiments involve Ramsey interferometers and the 
aim was to demonstrate spin squeezing created via the collisional interaction 
between the bosons. Obviously, in demonstrating spin squeezing they would 
hope to have created an entangled state, though in most cases an entangled 
state had already been created via the interaction with the first beam split¬ 
ter. Although the criterion for entanglement used in most cases was based on 
an experimental proposal [13) . [31j which regarded identical particles as distin¬ 


guishable sub-systems, the spin squeezing test based on Sz does turn out to be 
a.valid test for two mode entanglement, as explained in Section However, 
it should be noted that all the papers discussed have a different viewpoint re¬ 
garding what exactly is entangled - generally referring to entanglement of atoms 
or particles rather than modes. All the experiments discussed below establish 
entanglement, though often this was already created in a hrst tv/2 coupling 
pulse. Most are based on the spin squeezing test involving Sz, that of Gross et 
al [39] involved population difference squeezing rather than spin squeezing (see 
Subsection 3.5). The other experiment of Gross et al [ID] shows (see Fig 2b in 
[40) ') that the mean value of one of the two spin operators S^, Sy is non-zero, 
as measurement results such as in (224) for the simple two mode interferometer 
with 2s = 9, (j) = 0 can determine. This is sufficient to demonstrate two mode 
entanglement, as (68) shows. 

A key result of the present paper (and [I]) is that the conclusion that experi¬ 


ments which have demonstrated spin squeezing in Sz have thereby demonstrated 
two mode entanglement, no longer has to be justified on the basis of a proof 
that clearly does not apply to a system of identical bosons. 


8.1 Esteve et al. (2008) |41j 

• Stated emphasis - generation of spin squeezed states suitable for atom inter¬ 
ferometry, demonstration of particle entanglement. 

• System - Rb®^ in two hyperfine states. 

• BEG of Rb®^ trapped in optical lattice superposed on harmonic trap. 

• Occupation number per site 100 to 1100 atoms - macroscopic. 

• Situation where atoms trapped in just two sites treated - two mode entan¬ 
glement? 

• Glaimed observed (see Fig 1 in |1T|) spin squeezing based on N (^AS'^'^ + 

< 1 (see 


(1321). 
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• Claimed entanglement of identical atoms. 

• Spin squeezing test is based on assumption that Bloch vector is on Bloch 


sphere, a result not established. 

• Comment - spin squeezing in 
entanglement is established. 


(almost) demonstrated (see (132)), so 


8.2 Riedel et al. (2010) |42j 


• Stated emphasis - generation of spin squeezed states suitable for atom inter¬ 
ferometry, demonstration of multi-particle entanglement. 

• System - Rb®^ in two hyperfine states. 

• BEC of Rb®^ trapped in harmonic trap with non-zero magnetic field - 
Zeeman splitting of levels. 

• Number of atoms 1200 - macroscopic. 

• Process involves Ramsey interferometry - starts with all atoms in one 

state, 7r/2 pulse (duration 7r/2n ?) generates coherent spin state (a'^ -I- |0) 

(entangled), free evolution with collisions (causes squeezing), second pulse with 
area 2s = 0 and phase tt (or 0) followed by detection of population difference - 
associated with operator S^- 

• Evolution described using Josephson Hamiltonian H = 6Sz + xSt 

where Scf, = cos 4> Sx — sin (p Sy,Vt is Rabi frequency, p is phase of RF-microwave 
field, 5 is detuning, y describes collisions. Interaction picture and on resonance 
? 


• During free evolution including collisions spatial modes for internal states 
pushed apart so that y becomes much bigger in order to give larger squeezing. 

• Final pulse enables state tomography in the yz ^\a.iie to be carried out - 
measures spin squeezing for spin operator Sg = cos 9 Sz — sin 0 Sy in this plane 
(see (236) herein with p = n). 


ASg ) being less than standard 


• Claimed observed spin squeezing based on 
quantum limit N/4. (see Fig 2 in [12]). 

• No measurement made to show that | | « N/2 as required to justify 

spin squeezing test. Spin squeezing test is based on assumption that Bloch 
vector is on Bloch sphere, a result not established. 

• Claim that state of atoms at end of free evolution is four-partite entangled 
based on spin squeezing test is not substantiated, also an entangled state was 
already created by first n/2 pulse. 

• Comment - spin squeezing in Sz (almost) demonstrated (see (132)), so 
entanglement is established. An entangled state was of course already created 
by first n/2 pulse, and then modified via the collisional effects. 


8.3 Gross et al. (2010) | 3D j 

• Stated emphasis - generation of non-classical spin squeezed states for non¬ 
linear atom inteferometry, demonstration of entanglement between atoms. 

• System - Rb®"^ in two hyperfine states. 
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• Six independent BECs of Rb®^ trapped in six separate wells in a optical 
lattice. 

• Number of atoms 2300 - macroscopic, down to ca 170 in each well. 

• Evolution described using Josephson Hamiltonian H = AuJoSz+rtS^+xS"^ 
where (in the present notation) Sj = cos 7 Sa,+sin j Sy, H is Rabi frequency, 7 is 
phase of RF-microwave field, Awq is detuning, x describes collisions. Interaction 
picture and on resonance ? 

• During free evolution plus collision evolution Feshach resonance used so 
that X becomes much bigger in order to give larger squeezing. 

• One process involves Ramsey interferometry - starts with all atoms in one 

state, 7r/2 pulse (duration 7 r/ 2 f 2 ?) generates coherent spin state (d'^ + |0) 

(entangled) with (^Sz^ = 0 , free evolution with collisions (causes squeezing) 
and with spin echo pulse applied, second 7r/2 pulse with phase (j) followed by 
detection of population difference - associated with operator Sz- 

• Population difference measured after last 7 r /2 pulse shows a sinusoidal 

dependence on phase 4> (see Fig 2b in [10]). This shows that and 

are non-zero, thereby showing that the state created just prior to last pulse is 
entangled (see Bloch vector test ( 681 ). This does not of course show that the 
state is spin squeezed. 

• Another process involves generation of coherent spin state (d^ -I- |0) 

(entangled) with (^Sz^ = 0 , then free evolution with collisions (causes squeez¬ 
ing), followed by coupling pulse to rotate Bloch vector through angle a thereby 
crossing xy plane. The variance in Sz is then measured as a changes. 

• Claimed observed spin squeezing based on TV /{(^Sx 

being less than 1 (see Fig 3 in |iD]l. 

• Spin squeezing test is based on assumption that Bloch vector is on Bloch 

sphere, a result not established since and (^Sy'j not measured. 

• Claimed entanglement of ca 170 atoms. 

• Comment - spin squeezing in Sz (almost) demonstrated (see (132)), so 
entanglement is established. An entangled state was of course already created 
by first tt/2 pulse, and then modified via the collisional effects. 


+ (^v) ) 


8.4 Gross et al. (2011) [39] 

• Stated emphasis - continuous variable entangled twin-atom states. 

• System - Rb®"^ in several hyperfine states. 

• Independent BECs of Rb®^ trapped in separate wells in a optical lattice. 

• Number of atoms macroscopic, ca few 100 in each well. 

• Spin dynamics in Zeeman hyperfine states (2,0), (1,±1). 

• Initially have BEC in (2,0) hyperfine state - acts as pump mode. 

• Spin conserving collisional coupling to (1,±1) hyperfine states - which act 
as the two mode system. 
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• One boson created in each of (1, ±1) hyperfine states with two bosons lost 
from (2,0) hyperfine state due spin conserving collisions. 

• OPA type situation associated with spin changing collisions with (1,±1) 
hyperfine states acting as idler, signal modes . 

• Mean and variance of population difference between (1,+1) and (1,-1) 
hyperfine states measured. Total population also measured. 

• Entanglement test is that if the variance in population difference is small, 
but that in the total boson number is large then the state is entangled (see (73), 

(HI)- 

• Measurements (see Fig Ic in [33]) show noise in population difference is 
small, but that in the total boson number is large. 

• Further entanglement test is that if there is two mode quadrature squeezing 
then the state is entangled. 

• Comment - Number squeezing and two mode quadrature squeezing demon¬ 
strated and entanglement confirmed. 
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9 Discussion and Summary of Key Results 

The two accompanying papers are concerned with mode entanglement for sys¬ 
tems of identical massive bosons and the relationship to spin squeezing and 
other quantum correlation effects. These bosons may be atoms or molecules as 
in cold quantum gases. The previous paper I dealt with the general features 
of quantum entanglement and its specihc definition in the case of systems of 
identical bosons. In defining entanglement for systems of identical massive par¬ 
ticles, it was concluded that the single particle states or modes are the most 
appropriate choice for sub-systems that are distinguishable, that the general 
quantum states must comply both with the symmetrisation principle and the 
super-selection rules (SSR) that forbid quantum superpositions of states with 
differing total particle number (global SSR compliance), and that in the separa¬ 
ble states quantum superpositions of sub-system states with differing sub-system 
particle number (local SSR compliance) also do not occur [1]. The present pa¬ 
per II has examined possible tests for two mode entanglement based on the 
treatment of entanglement set out in paper I. 

The present paper first defines spin squeezing in two mode systems for the 
original spin operators Sx,Sy,Sz, which are defined in terms of the original 
mode annihilation and creation operators d, b and , 6^. Spin squeezing for the 
principal spin operators Jx,Jy,Jz for which the covariance matrix is diagonal, 
rather than via the original spin operators is then discussed. It is seen that the 
two sets of spin operators are related via a rotation operator and the principal 
spin operators are given in terms of new mode operators c, d and '(? ,d^, with c, d 
obtained as linear combinations of the original mode operators a, b and hence 
defining two new modes. Finally, we consider spin squeezing in the context of 
multi-mode systems. 

The consequence for the case of two mode systems of identical bosons of 
the present approach to defining entangled states is that spin squeezing in any 
of the spin operators Sy or Sz requires entanglement of the original modes 
d,b. Similarly, spin squeezing in any of the new spin operators Jx, Jy or Jz 
requires entanglement of the new modes c, d. The full proof of these results has 
been presented here. A typical and simple spin squeezing test for entanglement 
is (^ASx'^ < I 1/2 or (^ASy'^ < \ (^Sz^ |/2. We also found a simple Bloch 

vector test 0 or 0. It was noted that though spin squeezing 

requires entanglement, the opposite is not the case and the NOON state provided 
an example of an entangled physical state that is not spin squeezed. Also, 
the binomial state provided an example of a state that is entangled and spin 
squeezed for one choice of mode sub-systems but is non-entangled and not spin 
squeezed for another choice. The relative phase state provided an example 
that is entangled for new modes c, d and is highly spin squeezed in Jy and 
very unsqueezed in Jx- We then showed that in certain multi-mode cases, spin 
squeezing in any spin component confirmed entanglement. In the multi-mode 
case this test applied in the bipartite case ( Case 1 ) where the two sub-systems 
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each consisted of all the modes or all the modes hi or in the single modes case 
{Case S)where there were 2n sub-systems consisting of all the modes at and all 
the modes hi. For the mode pairs case {Case 3) where there were n sub-systems 
consisting of all the pairs of modes Oi and hi, a spin squeezing entanglement test 
was found in the situation where for separable states each mode pair involved 
a single hoson. The connection between spin squeezing and entanglement was 
regarded as well-known, but up to now the only existing proofs were based on 
non-entangled states that disregarded either the symmetrization principle or the 
sub-system super-selection rules, placing the connection between spin squeezing 
and entanglement on a somewhat shaky basis. On the other hand, the proof 
given here is based on a definition of non-entangled (and hence entangled) states 
that is compatible with both these requirements. 

There are several papers that have obtained different tests for whether a state 
is entangled from those involving spin operators, the proofs often being based 
on a definition of non-entangled states that ignores the sub-system SSR. Hillery 
et al |23] obtained criteria of this type, such as the spin variance entanglement 
test + (^ASy'^ < I The proof of this test has also been set out 

here, and the test is also seen to be valid if the non-entangled state definition 
is consistent with the SSR. The test -f (^ASy'^ < \ I suggested 

by the requirement that (^AS^j + (^ASy'j > \ \ for non-entangled states 

- since both (^AS'^'j > \ 1/2 and (^ASy'^ > \ 1/2 is of no use, since 

(^AS'^'^ + l^ASy'^ > I I for all states. However as previously noted, showing 

that either (^AS'lj < \ ^5'^^ |/2 or (^ASy'j < \ |/2 - or the analogous tests 

for other pairs of spin operators - already provides a test for the entanglement of 
the original modes a, h. This test is a different test for entanglement than that of 
Hillery et al [23]. In fact the case of the relative phase eigenstate is an example 
of an entangled state in which the simple spin squeezing test for entanglement 
succeeds whereas that of Hillery et al [23] fails. The consequences of applying 
both the simple spin squeezing and the Hillery spin operator test were examined 
with the aim of seeing whether the results could determine whether or not the 
local particle number SSR applied to separable states. The conclusion was 
negative as all outcomes were consistent with both possibilities. In addition, the 
Hillery spin variance test was also shown to apply to the multi-mode situation, 
in the Cases I and 2 described above, but did not apply in Case 3. Other 
entanglement tests of Benatti et al [50] involving variances of two mode spin 
operators were also found to apply for identical bosons. 

The present paper also considered correlation tests for entanglement. In¬ 
equalities found by Hillery et al [33] for non-entangled states which also do not 
depend on whether non-entangled states satisfy the super-selection rule include 
I ^(a)™ (6^)"^ P < giving a valid strong correlation test 

I ^(a)™ (6^)"^ P > for an entangled state. However, with 
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entanglement defined as in the present paper we have | P = 0 for 

a non-entangled state, so we have also proved a weak correlation test for entan¬ 
glement in the form | ^(a)™ P > 0. This test is less stringent than the 

strong correlation test of Hillery et al [33]., as | ^(a)™ P is then required 

to be larger. In all these cases For n ^ m none of these cases are of interest 
since for global SSR compliant states ^(a)™ would be zero. In the cases 

where n = m we show that all the correlation tests can be expressed in terms of 
spin operators, so they reduce to tests involving powers of spin operators. For 
the case n = m = 1 the weak correlation test is the same as the Bloch vector 
test. 

Work by other authors on bipartite entanglement tests has also been exam¬ 
ined here. He et al [33|, [33] considered a four mode system, with two modes 
localised in each well of a double well potential. If the two sub-^stems A and B 
each consist of two modes - with oi, 0,2 as sub-system A and bi , 62 as sub-system 
B, then tests of bipartite entanglement of the two sub-systems of the Hillery 
|33| type I for any i, j = 1 , 2 or 

involving local spin operators | Sf Sff'j apply. Raymer 

et al m have also considered such a four mode system and derived bipartite 
entanglement tests such as ^A(5' ^ ± S + ^A(S' y T S < \ (*^ 2 ) I 
involve local spin operators for the two sub-systems. 

We also considered the work of Sorensen et al [T3|, who showed that spin 
squeezing is a test for a state being entangled, but defined non-entangled states 
for identical particle systems (such as BECs) in a form that is inconsistent 
with the symmetrisation principle - the sub-systems being regarded as indi¬ 
vidual identical particles. However, the treatment of Sorensen et al m can 
be modified to apply to a system of identical bosons if the particle index i is 
re-interpreted as specifying diffferent modes, for example modes localised on 
optical lattice sites i = I,2,..,n or localised in momentum space. With two 
single particle states \4>ai)j \4>bi) with annihilation operators ai,bi available on 
each site, there would then be 2 n modes involved, but spin operators can still 
be defined. This is just a particular case of the multi-mode situation described 
above. If the definitions of non-entangled and entangled states in the present 
paper are applied, it can be shown that spin squeezing in either of the spin 
operators Sx or Sy requires entanglement of all the original modes ai,bi (Case 
I, above). Alternatively, if the sub-systems are pairs of modes ai,bi and the 
sub-system density operators were restricted to states with exactly one bo¬ 
son, then it can be shown that spin squeezing in Sz requires entanglement of 
all the pairs of modes (Case 3, above). With this restriction the pair of modes 
Qi, bi behave like distinguishable two state particles, which was essentially the 
case that Sorensen et al m implicitly considered. This type of entanglement 
is a multi-mode entanglement of a special type - since the modes di, bi may 
themselves be entangled there is an ’’entanglement of entanglement”. So with 


97 


either of these revisions, the work of Sorensen et al m could be said to show 
that spin squeezing requires entanglement. However, neither of these revisions 
really deals with the case of entanglement in two mode systems, and here the 
proof given in this paper showing that spin squeezing in requires entangle¬ 
ment of the two modes provides the justification of this result without treating 
identical particles as distinguishable sub-systems. Sorensen and Molmer m 
have also deduced an inequality involving ( AS ^ ) and | ( S' 2 ) I for states where 


~ \ y) ~ ^ based on just the Heisenberg uncertainty principle. This 
is useful in terms of confirming that states do exist that are spin squeezed still 
conform to this principle. 

Entanglement tests involving quadrature variables have also been published, 
so we have also examined these. Duan et al [28], Toth et al [34| devised a test for 
entanglement based on the sum of the quadrature variances (A(a;^ ± Xb)"^) + 
(^A{pa TPb)^) > 2 for separable states, which involve quadrature components 
xa,Pa,xb,Pb constructed from the original mode annihilation, creation opera¬ 
tors for modes A, B. Their conclusion that if the quadrature variances sum is 
less than 2 then the state is entangled is valid both for the present definition of 
entanglement and for that in which the application of the super-selection rule 
is ignored. However, for quantum states for systems of identical bosons that 
are global SSR compliant (^A{xa ^xb)"^) + {A{paT Pb)'^) = ‘2 + 2^-^^ for 
all such states - both separable and entangled, and although this is consistent 
with [28j . [34] we have concluded that the quadrature variance test can never 
confirm entanglement. A more general test [35] involving quadrature operators 


A® , Xg required showing that 


vB v<t> 
-^A -^B 


^ 0. This was shown to be equivalent 


to showing that yf 0 or 7^ 0, the Bloch vector or weak correlation 

test. A two mode quadrature squeezing test was also obtained, but found to be 
less useful than the Bloch vector test. 

Overall then, all of the entanglement tests (spin squeezing and other) in 
the other papers discussed here are still valid when reconsidered in accord with 
the definition of entanglement based on the symmetrisation and super-selection 
rules, though in one case Sorensen et al m a re-definition of the sub-systems 
is required to satisfy the symmetrization principle. However, further tests for 
entanglement are obtained in the present paper based on non-entangled states 
that are consistent with the symmetrization and super-selection rules. In some 


cases they are less stringent - the correlation test in Eq.(165) being easier to 
satisfy than that of Hillery et al |33j in Eq. (168). The tests introduced here 


are certainly different to others previously discovered. 

The theory for a simple two mode interferometer was then presented and 
it was shown that such an interferometer can be used to measure the mean 
values and covariance matrix for the spin operators involved in entanglement 
tests for the two mode bosonic system. The treatmebnt was also generalised to 
multi-mode interferometry. The interferometer involved a pulsed classical field 
characterised by a phase variable </> and an area variable 2s = 9 defined by the 
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time integral of the field amplitude, and leads to a coupling between the two 
modes. For simplicity the centre frequency was chosen to be resonant with the 
mode transition frequency. Measuring the mean and variance of the population 
difference between the two modes for the output state of the interferometer for 
various choices of (p and 9 enabled the mean values and covariance matrix for 
the spin operators for the input quantum state of the two mode system to be 
determined. More complex interferometers were seen to involve combinations of 
simple interferometers separated by time intervals during which free evolution 
of the two mode system can occur, including the effect of collisions. 

Experiments have been carried out demonstrating that spin squeezing oc¬ 
curs, which according to theory requires entanglement. An analysis of these 
experiments has been presented here. However, since no results for entangle¬ 
ment measures are presented or other independent tests for entanglement car¬ 
ried out, the entanglement presumably created in the experiments has not been 
independently confirmed. 
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11 Appendix A - Spin Squeezing Test for Bipar¬ 
tite Multi-Mode Case 


We now consider spin squeezing for the multi-mode spin operators given in Eqs. 
§ and (§ in Subsection |2.2| We consider separable states for Case 1, the 
density operator being given in Eq. (|62[). In this bipartite case the two subsys¬ 


tems consist of all modes ai and all modes bi. The development involves ex¬ 
pressions such as = Trci^l(r)p^) and 

c c 

= Trc($Kr)$Kr')?K), ($c(r)$c(r'))^ = rrc($c(r)$c(r')pg), 
c 


^ = Trc{^l{r)^c{r')^), where C = A,B. 


Firstly, we have 


S, 


= ^ f dr 

R 2-’ 


§t(r))^ 0 (252) 


since from the local particle number SSR for sub-systems A and B we have 


follows that 


= ^vl/^(r)^ =0. A similar result applies to so it then 


S.) = (Sy)=0 


(253) 


This immediately yields the Bloch vector entanglement test. It also leads to the 
spin squeezing in S' ^ entanglement test, namely if S is squeezed with respect 
to S X or S y (or vice versa), then the state must be entangled. The question 
then is: Does spin squeezing in S^, with respect to Sy (or vice versa) require 
the state to be entangled for the two n mode sub-systems A and Bl 
To obtain an inequality for the variance in Sa,, we see that 


S; 


^ffdrdr' X {(^§l(r)§l(r')'^^ (^$„(r)$„(r'))^ 

+ ($b(i’)$fc(i-') 


B 


A 


B 


H, \ / ix \ ' R. 

-h($b(r)$t,(r')\ ($l(r)$);(r')\ } 


$l(r)$a(r') 


R 


R 


(254) 


From the local particle number SSR for sub-systems A and B we have ^'h|)(r)4'j)(r')^ = 

^ ^ \A 

^a(r)^a(r') ) = 0, SO the first and fourth terms are zero. Using the field op- 

' 

erator commutation rules we then obtain 




+ \ldr { ^$];(r)$b(r)^^ -h (|$l(r)$a(r)^^} (255) 
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so that 



Hence from (32) 


\ II dr dr' (|$l(r')$a(r)^^ 

+ \ldr { (^I(i’)$fe(r)^^ + ^$l(r)$a(r)^^} 


(a52) 

> 12Pr{\II drdr' ^$j;(r)$b(r'))^ ($I(r')$a(r))^ 

+ \ldr { 

The same result applies to (^ASy'j. 

Now we can easily show that 

= E-Pfl^/c^r{($I(r)$6(r))^ - (|$l(r)$a(r)^^} 

so that 

l\{s.)\<EPR\ldr{ ($t(r)$,(r))^ + ($t (r)$ Jr))^} 


(256) 


(257) 


(258) 

(259) 


as ^'l'j)(r)'I'h(r)^ and (r)4'Q(r)^ are real and positive. 

Hence we find that 

> (|$];(r)$b(r'))^(|$l(r')$a(i-))^ ( 260 ) 

= 12 PrI II dr dr' TrB{I>b{r')p^ $j;(r)} rrA{$a(r) §l(r')}( 261 ) 

R ^ 

= ^//c?i’dr'rr|$a(r)$b(r')p^ep$l(i’')^b(i’)} ( 262 ) 


giving three forms that the inequality for ~ sl I satisfy in 

the case of a separable state. The last form involves a double space integral of a 
quantum correlation function. Note the order of r and r'. It is straightforward to 
show that the right side of the inequality is real, but to achieve an entanglement 
test involving spin squeezing for S' 2 , we need to show that it is non-negative. 
Identical inequalities can be found for ( AS y)~5l(5'z)|. 
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11.0.1 Mode Expansions 


If we use Eq.([^ to expand the field operators then using Eq.(2611 we have 


R ^ ij kl 

X ^TrA{ai'pRa]}TrB^k'pR^}^ 

= T.Pr\ E(^S Af^ (263) 

R 4 y 

= + (264) 

R 4 


where mode orthogonality has been used and we have introduced matrices 
and whose elements are 

Afj = TrA{a, Pr a]} = Trsib^ Pr (265) 

It is easy to show that Ai^ = {A^)* and = (B^)* showing that the matrices 

A^ and B^ are Hermitian, as is A^B^ + B^A^. The quantity E(^ii B^ + 

ij 

B^ Afj) is real. The question is: Is it also positive 1 

For the simple case where there is only one spatial mode for each component 

the right side of the inequality is just equal to a’^} 

R '■ 

^ Pr\N^ , where and give the mean numbers of bosons in sub- 
R 

■^A 

systems A and B for the states and p^- The right side of the inequality is 
positive, showing that the separable state is not spin squeezed for S x with re¬ 
spect to S'y (or vice versa), leading as before to the test that such spin squeezing 
requires entanglement. 


11.0.2 Positive Definiteness 

For the multi-mode case we now take into account that the sub-system density 
operators 'p^ and are positive-definite. Their eigenvalues and are 
real and non-negative as well as summing to unity, and we can write the density 
operators in terms of their orthonormal eigenvectors \AR,X) and \BR,fi) as 

PR = E^t^ARA){AR,X\ p^=j:7r^^\BR,p){BR,p\ (266) 

A 
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Then from (2651 


= E = E (BR, Ml bjb, \BR, /r) 

A M 

(267) 

Consider a 1 x n row matrix = {■^ 1 ,^ 2 ) - -jCn } 


= E E -^1 «]«* 

A ij 

= J2ni^ {AR,X\ h\hA\AR,X) (268) 

A 

where we have introduced the operator ^Ia = Since ^^A^^ is always 

non-negative for all this shows that A^ is a positive definite matrix. Similarly, 
considering a 1 x n row matrix 77 I = {rjl, p^T-^fin } introducing the operator 
fts = Vi bi we find that 

^ {BR, ^i\ ns \BR, p) (269) 


which is also always non-negative, showing that B^ is also a positive definite 
matrix. 

We can then express the positive definite Hermitian matrices A^ and B^ 
in terms of their normalised column eigenvectors 0^ and respectively, where 
the corresponding real, positive eigenvalues are Va and ajs Thus we have (for 
ease of notation R will be left understood) 


9i 
B^ Cf 


Then 


( 0 ^< = <5„, A^ =J2’^XXV 

a 

(Cf )^cf = = T. (Cf )^ (270) 


Br{Y^ E Oi (0f )^Cf (Cf )^} 

a /3 

+^^{E E cf (cf X)'’} 

a /3 

EE^“^/3 [(^ffcf] [iXot] 

ct /3 

+EE^“^/3 [(cf)^^f] Pffcf] 

a p 

2EE^«'^/3i[(^ffcf]i^ 

a p 


(271) 
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Hence we have using (264) 



( 272 ) 


where the right side of the inequality is non-negative. The same result applies 
to (^ASy'^ — 11 I- Thus separable states are not spin squeezed in S' 2 : or 
in Sj,. 

Thus we have established the spin squeezing test for the multi-mode Case 
1 - states that are spin squeezed in S^; compared to Sy.{or vice versa) must be 
entangled states for the two subsystems consisting of all modes and all modes 

h- 

For the other spin components, the Bloch vector result in 
S y'j = 0 for separable states enables us to show that if is squeezed com¬ 
pared to S' 2 ;.(or vice versa) or if is squeezed compared to Sy.{or vice versa) 
then the state must be entangled. Thus spin squeezing in any spin component 
requires the state to be entangled, justt as for the two mode case. 


253) that (Sx) = 
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12 Appendix B - Spin Squeezing Tests for Other 
Multi-Mode Cases 


12.1 Single Mode Sub-Systems 


We now consider spin squeezing for the multi-mode spin operators given in Eqs. 
§ and in Subsection |2.2| We consider separable states for Case 2, the 
density operator being given in Eq. (|63|). In this single mode sub-system case 


there are 2n subsystems consist of all modes dt and all modes hi. 

This case is that involved in the modihed approach to Sorensen et al and we 

ot.(aS'I) that applies 


will see that it leads to a useful inequality for {^AS 


when non-entangled states are those when all the separate modes di and bi are 
the sub-systems . We will follow the approach used for the simple two mode 
case in Section |3l 

Firstly, the variance for a Hermitian operator 12 in a mixed state 


P = '^PrPr 

R 


(273) 


is always greater than or equal to the the average of the variances for the separate 
components 

>^Pr (^Ah^^ (274) 

R ^ 

where = Tr(p AQ. ^) with AH = 12—^17^ and ^A12 = Tr(pj^ AQr ^) 

with AQr = . The proof is straight-forward and given in Ref. [20] , 

Next we calculate (^ASl'^ , (^ASy'^ and 
case where 

P = ^ Pr {pr ® ® {pr ® Pr'^ ® {pr ® Pr^ ® •••■ (275) 

R 

as is required for a general non-entangled state all 2n modes. This situation 
is that of Choice 2 for the sub-systems, as described in SubSection |3.3[ As 
the density operators for the individual modes must represent possible physical 
states for such modes, so super-selection rule for atom number applies and we 
have 

= Tr(p“*(a,)P) = 0 (idlY) =Tr{rR(alr)=0 

= Trip^Yihr) = o (^iblr)^^ = TriYYi^r) = o 

(276) 
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The Schwinger spin operators are 


i i 

Sy = ^(Sla,-al6.)/2* = E^y 

i i 

S. = Y.(^%-a%)/2 = Y,Sl 


( 277 ) 


where ai, bi and a), 6| respectively are mode annihilation, creation operators. 
Note that this expression for the spin operators is the same as §) for the multi- 
mode case treated in SubSection 2.2 From Eqs. (277) we find that 


I 


(278) 


so that on taking the trace with 'pj^ and using Eqs. (275) we get after applying 
the commutation rules [e, e^] = 1 (e = a or 6) 




As we also have 






Si 

R \ / R 


using Eqs. (275) and we see finally that the variance (ASi) is 


(As»)^=E((af)„-E(a,„ 


all the terms with i ^ j cancelling out. and therefore from Eq. (|274|) 

E f((5;)2) - ( 5 ;' 


(279) 


(280) 


(281) 


(282) 


An analogous result applies for ( AS y 


But using (276) 


{Sir = 


(Si) 


A 


R 


-(SJoi^loi + aJSiSJai -I- afbia]bi) 

+ ((2ta).)^) + i(((at«),)^ (283) 


and 


S':) =0 

R 


(284) 
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It then follows that 


R 


R 


(285) 


so that 


^Sl)>Y,PnE 

(286) 


The same result applies for ( AS y 


Now using (276) 






= 

R * 

R * 

< E^«4SI(((6%)^-((«^«)0fi))l 

R ® 

< E^« 4 S(((^%) +((«^«).)«)) (288) 


and thus 


ASi)-^\(S 


> E E ()(((*’*>•)„ + <P’s).>h) + )((P’S).>H ((*’‘).)„) 

R 

-E^«JS(((&%)^ + ((a^2).)^)) 


> 0 


(289) 


A similar proof shows that (^ASy'^ — 5 I ^«S'z^ | > 0 for the non-entangled state 
of all 2 n modes. 

This shows that for the general non-entangled state with all modes a-i and bi 
as the sub-systems, the variances for two of the spin fluctuations and 

ASy'j are both greater than i| (^Sz'j |, and hence there is no spin squeezing for 


no 



Sx or Sy. Note that as 


= 0, the quantity 




+ 1 5-L 


is the same as | |, so the alternative criterion in Eq. (16) is the same as 

that in Eq. © which is used here. 

Hence we have shown that for a non-entangled physical state for all the 2n 
modes di and k 


as: 


> 


s. 


and ( AS ' 


> 


(290) 


so that spin squeezing in either Sx or Sy requires entanglement. 

From (285) we see that x'^ = = 0 for the general separable state 

showing there is a Bloch vector test for entanglement such that if either (^S 
or y'j is non-zero, then the state must be entangled. 

Finally, if there is spin squeezing in Sz with respect to Sx or vice versa, or 
spin squeezing in Sz with respect to Sy or vice versa, it follows that one of (^S ^ 

or is non-zero. But as both these quantities are zero for a non-entangled 

state, if follows that spin squeezing in Sz also requires entanglement. 

Thus, spin squeezing in any spin operator Sx,Sy or Sz is a sufficiency test 
for entanglement of all the separate mode sub-systems. 


12.2 Two Mode Sub-Systems 

We now consider spin squeezing for the multi-mode spin operators given in Eqs. 
([^ and ^ in SubSection 2.2 We consider separable states for Case 5, the 
density operator being given in Eq. (64). In this mode pair sub-system case 


there are n subsystems consist of all pairs of modes di and bi. 

This case is also involved in a modified approach to Sorensen et al and we 
show a useful inequality for ( AS I) applies when non-entangled states are those 


when the pairs of modes di and bi are the separate sub-systems, but only in 
restricted situations. The pairs of modes corresponding to localised modes on 
different lattice sites or pairs of modes with the same momenta doe represent the 
closest way of simulating the approach used by Sorensen et al where identical 
particles i were regarded as the sub-systems. 

Now the general non-entangled state will be 




Pr 


'Pr^ 


'fR' 


(291) 


where the are now the density operators for sub-syst em i consisting of the pair 
of modes di and bi (which are of the form given in Eq. (138)) and the conditions 
in Eq. (276) no longer apply. The Fock states are of the form \Nia)®\Nn,) for the 
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pair of inodes di and bi, and for this Fock state the total occupancy of the pair 
of modes is Ni = Nia + Nn,. From the super-selection rule the density operator 
for the *th pair of modes di and bi is diagonal in the total occupancy. For 
Ni = 0 there is one non zero matrix element ((0|j^ (8) (0|-j,) (|0)0 |0)- j). For 

Ni = 1 there are four non zero matrix elements, which may be written 


((lLa®(0L6)Pfl(|l)»a®|0)*h) = Paa 

((l|*a®(0L6)Pfl(|0)»a®|l)*h) = Plb 

((0La®(lL6)Pfl(|l)»a®|0)*h) = Pla 

((0La®(lL6)Pfl(|0)»a® |l)*h) = Plb 


(292) 


For Ni = 2 there are nine non zero matrix element ((2|j^0(O|-^) p)j (|2)j^010)-j,), 
((0|ia ® (2|ib)Pfl (|0)io 0 |2)j^) and the number increases with Ni. 

If we restrict ourselves to general entangled states for one particle, where Ni 
= 1 for all pairs of modes, then the density operator p)^ is of then form 


Pfl = Paa(|l).a(lL®|O).5(O|J+P*„6(|l),,(OL0|O),,(l|J 

+Plam.a (l|,a ® |1),. Kb) + PM (0|.a ® |l).b Kt) (293) 


In addition Hermitiancy, positivity, unit trace Tr('p]^) = I and < 1 can 

be used as in Eq (3271 to parameterise the matrix elements in (292). 


Ptb = cos ai 

= sin^ Oi cos^ Ui sin^ exp(+i0j) 


Paa = Sin a, 
Pab 


pla = y sin^ Oi cos^ ai sin^ exp(— 

(294) 


The expectation values for the spin operators 5*, S'® and S® associated with 
the zth pair of modes are then 

= Tr{p^n^{V.ai + a^fy 
i^^)R ^ J-{Pab-Pba) 

(Sl)^ = l{plb-Pla) (295) 


as in the Appendix 14 


which are of exactly the same form as in Eq. 
derivation of the original Sorensen et al m results based on treating identical 
particles as the sub-systems. The proof however is now different and rests on 
restricting the states p^ to each containing exactly one boson. 

The remainder of the proof is exactly the same as in Appendix and we 
find that 

ASl)>^({s,y + {Syy) (296) 
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for non-entangled pairs of modes Oi and bi. Thus when the interpretation is 
changed so that are the separate sub-systems are these pairs of modes, it follows 
that spin squeezing in with respect to S'a, or 5" y requires entanglement of 
all the mode pairs, but only if there is one particle in each mode pair. 

In general, spin squeezing in either S ^ or Sy is not linked to entangle¬ 
ment for Case 3 sub-systems, as has been pointed out in SubSection |3.3| by a 
counter-example involving the relative phase state. Also there is no Bloch vector 
entanglement test. For we have in general 

= Tr{p^ji^(blai + afb,) 

(Sl)^ = rr(p^i(&|a,-al60 (297) 

and the local particle number SSR does not require these quantities to be zero 
for sub-systems consisting of pairs of modes and bi. Thus in general and 

can be non-zero for a separable state, so the Bloch vector entanglement 
test does not apply. 
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13 Appendix C - Hillery Spin Variance - Multi- 
Mode 


13.1 Bipartite Case 


We first consider Case 1 where there are two sub-systems each consisting of all 
the modes or all the modes bi. We use the results from (2571 to find that for 
a separable state 


(as») + (as>) 

> E^Mlfdrdr' /^(r)$b(r')\ /^l(r')^a(r)) 

\ / R \ / H 

+ \ldr { ^ (298) 

The same result would have occured if the local sub-system SSR had been 
disregarded, the terms such as ^ ff dr dr' x 
cancelling out. 

The mean number of bosons is obtained from ([^ and hence 

^ (^) = +($l(r)$a(r)) (299) 

R ^ ^ 

Thus we have 

- 2 (^) 

> (^b(j')^b(i’'))^ (^l(i'')^a(r)) (300) 

^ \ / H \ / H 

Using the mode expansion ([^ we then get 

(aU) + (as;) - i («) 

R ki \ lR\ ! R 

= EPR'^TrA{d,'pjid\}TrB^j'pRVi} (301) 

^ ij 

= EPR\Tr{A^B^ + B^A^) (302) 

R 2 


after orthogonality is used and the matrix elements Af- and Sjl are introduced 
from (265). 

Since we have shown in Appendix [H] that the right side of the last inequality 
is always non-negative, the Hillery spin variance entanglement test follows that 
if 

'asI) + (asI)<Un) (303) 
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then the quantum state must be an entangled state for the case of two sub¬ 
systems each consisting of all the modes or all the modes bi. 


13.2 Single Modes Case 


We now consider separable states for Case 2, the density operator being given in 
Eq. ( |63[ ) . In this single mode sub-system ease there are 2n s ubsy stems consisting 
of all modes Oi and all modes bi. We use the results from (286) to find that for 
a separable state 


(a5' 


{asI 




The same result would have occured if the local sub-system SSR had been 
disregarded, the terms such as I 


(305) 


A \b]b]) cancelling out. 

\ / Ri 

The mean number of bosons is obtained from ([^ 


Thus we have 


AS: 


AS: 


--2W 


R 


(306) 


which is always non-negative. 

The Hillery spin variance entanglement test follows that if the inequality in 
(303) occurs then the quantum state must be an entangled state for the case of 
2 n sub-systems consisting of all the modes ai and all the modes bi. 


13.3 Two Modes Case 

We now consider separable states for Case 3, the density operators being given 


in Eq. (64). In this two mode subsystem case there are n subsystems consisting 


of all mode pairs Oi and bi. We consider a special separable state with just one 
term where 

- ^ ^ah(2) ^ ^ ^ ^ab(n) ^3^7^ 


P sep P 


) p ' b 




We use the results from (281) to find that 


(a5' 


{AS: 


Y.{{{ASlf) + ((ASlf)) 


(308) 
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where for a = x,y. This result did not depend on applying 

the local SSR. 

Now suppose each of the two mode states is an entangled state of the 
modes at and bi in which the Hillery spin variance test is satished. Then 

((A§;f) + ((A5;f)<i(n,) 

Hence 

(as|) + (as;) 

< 

= i («) (309) 


where N — ^ ■ Ni is the total number operator and Ni = blbi + djdi 

Thus the Hillery spin variance test is satisfied even though the state (3071 
is separable, shoeing that the test cannot be applied for multi-mode Case 3. 
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14 Appendix D - Derivation of Sorensen et al 
Results 


Sorensen et al |13j derive a number of ineqalities from which they deduce a 
further inequality for the spin squeezing parameter in the case of a non-entangled 
state. From this result they conclude that spin squeezing implies entanglement. 
The final inequality they obtain for a non-entangled state is 


A50>i 


2 /^ \ 2 
5,) +(S 


(310) 


Their approach is based on writing the density operator for a non-entangled 

(311) 


state of N identical particles as in Eq. (1301 


P — ^ Pr Pr® Pr® P^r' 

R 


■■■ = J2^rpr 

R 


The spin operators are defined as 

* * 

* 

Sz = = ^(|(/>6(i)) ((/'b(i)| - |(/'a(0) (<(>a(*)l)/2 


(312) 


where the sum i is over the identical atoms and each atom is associated with 
two states \4>g} and {(j)},). Clearly, the spin operators satisfy the standard com¬ 
mutation rules for agular momentum operators. 

Sorensen et al m state that the variance for Sz satisfies the result 


(Asb = f-EPiiE(S:)' _ 

4 R i ' ' R R 


EPR{Szy -{Sz'^ 


(313) 


To prove this we have 

’Si) = Y.PRTrCpj,Y.i:SlSi) 


E^«(e((* 1)’) +E(*l),(Si) 


R 

R VPl 


R 

M) 

R \ / R , 


(314) 
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where we have used 

^(l</'6(0) (</>b(*)l<^b(*)) - (!</'&(*)) (0b(*)ka(O) (</>a(*)l) 

+ ^(-(l<^a(0) (</>a(*)l'?^’b(*)) ('/'bWI + (l</>a(*)) ('/'a W l</>a(*)) ('/'a(*)l) 

('^bWI + (l</>aW)(<^a(*)l) 

it. (315) 

a result based on the orthogonality, normalisation and completeness of the states 
l</'a(0) >kb(0)- Also 

(S.) = TtCprY^SI) 

\ t H 2 

= ?(®)« 



so eliminating the term ( S (*^ 2 ) (^z 

Sion for - (s^'j . 

Next, Sorensen et al m state that 


gives the required expres- 


(^'<nj:prE(si)1 

\ / R i \ / R 

To prove this we have 


(5,)"<]v^p„El(5;) p 

R * 


R * 

R * 

R * 


(317) 


(318) 


since the modulus of a sum is less than or equal to the sum of the moduli. Now 



( 319 ) 
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using the general result that ( Pr V^r I <J2Pr Cfl, where ^ = 1 with 

\ R. Jr r 

here \/Cr = X) I I- Next consider 

y = "£ 1 ( 4 )^“ 

y-z = Ei\{si)\-\{si)\r> 0 (320) 

■ic^j ' 'IX ' 'IX 


so that 


<nJ2prE\{si 


Sy) <nJ2prE\{K 


R 


R 


(321) 


which is the required result. The inequality for (Sy) is proved similarly. 


Another inequality is stated [13] for . This is 

[s.y 


To show this we have 


S: 


R 


< 


R 

R * 

R 


(322) 


(323) 


so that 

{s.y = i(s.)t< 

R 

= ( 324 ) 
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using the general result that ( ^ P/? Cn j < Y^PrCr^ where ^ = 1 with 

\R Jr r 

herev^=|(§,)^|. 

Finally, we find that 

^ i" 

fl * ' ^ 

( 325 ) 

To show this we use the properties of the density operator for the ith particle 
of Hermitiancy, positiveness, unit trace Tr(p)^) = 1 and Tr(p]J)^ < 1. In terms 
of matrix elements of the density operator between the two states \4>ai'>‘))j 

the quantities , (^Sy'^ and £^re 

(Sl)^ = Tr{pn\m^)){Ui)\ + \U^)){M^m 

= 2 Pba) 

{^y}^ = J-iPab-Pba) 

(si)^ = l{plb-pL) ( 326 ) 

where p*^ = {4>ci'^)\'p‘R l^d(*))- The Hermitiancy and positiveness of p)j show 
that Pit, and pl^ are real and positive, pit, = (plJ* and p^Plfc- IPahP > 0- The 
condition Trilpp) = 1 leads to p),^ + pit, = 1, from which Tr(ffj^)‘^ < 1 follows 
using the previous positivity results. Taken together these conditions lead to 
the following useful parametrisation of the density matrix elements 

Pla = sin^ a, pit, = cos^ a, 

pit, = \J sin^ at cos^ at sin^ (3^ exp{+i(j)J) p(,^ = sin^ at cos^ at sin^ /?j exp(—i^j) 

( 327 ) 

where at, (3^ and (p^ are real. In terms of these quantities we then have 


(-5*)^ = ^sin2a, sin^/?. cos(/)i 

~ ^ sin 2Q;i sin^/3j sin (()j 

{Sl)^ = i cos2a, (328) 
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It is then easy to show that 


Si 


SI 


2 11 

= T — 7 sin^ (1 — sin"^ /S,) 

fi 4 4 


< 


( 329 ) 


and the final inequality (3251 then follows by taking the sum over particles i 


and then using — 1- If only Hi® Schwarz inequality is used instead of 

the more detailed consequences of Hermtiancy, positiveness etc it can be shown 

that (si\ + (si\ + 


R 


/ 


R 


SI) < |, which though correct is not useful. 


Combining the inequalities in Eqs. (317), (322) and (325) into Eq. (313) 
shows that 


AS”) 



N 


~4 

> 

N 

T 

> 

N 

T 

> 

1 

N 


R 


-EPrE(sI 

R i ' 

- -N 
4 


5, 


EPr{SzY 


(a 


(330) 


for the case of a non-entangled state. This result is that in Sorensen et al. m- 
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15 Appendix E - Heisenberg Uncertainty Prin¬ 
ciple Results 


Here we derive the results in SubSection |4.6| leading to inequalities for the vari¬ 
ance (AJU considered as a function of | ( J^) | for states where the spin op¬ 


erators are chosen such that \^Jxj = yj = 0- 

From the Schwarz inequality (j^ < so that 


+ + + + (331) 


giving Eq. (157|. Subtracting \^Jxj = vj — 0 from each side gives 


+ ('^•^/ Ad(J-l-l) 


(332) 


Substituting for from the Heisenberg uncertainty principle result in 


Eq. (158) gives 


ajiY - ("j(j + i)-(U\'Kaj? 


ii(d)Ao 


(333) 


The left side is a parabolic function of (^AJ^J and for this to be negative 
requires (AjU to lie between the two roots of this function, giving 


AJl) > 


X — 


1 <) 




(334) 


AJU < 




(335) 


which are the required inequalities in Eq. (159) and (160). 
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16 Appendix F - ’’Separable but Non-Local” States 
IS 

It is instructive to apply the various entanglement tests to the so-called 
separable but non-local states considered in Refs. 0, m. for which the sub¬ 
system states are definitely not SSR compliant. These states should not pass 
the the Hillery tests [23], [33] for SSR neglected entanglement, but they may 
pass the entanglement tests in this paper and in Ref. [Tj since these states would 
be regarded as SSR compliant entangled. Note that these states are consistent 
with the global particle number SSR, so there is no dispute about whether they 
are possible two mode quantum states. The issue is rather whether they should 
be categorised as separable or entangled, and that depends on how separable 
(and hence entangled) states are first defined. As discussed previously, the 
interferometric measurements discussed here do not enable us to choose one 
definition over the other - that is an issue involved swhat types of quantum 
states would be allowed in the separate sub-systems. 

The first example of such states is the mixture of two mode coherent states 
is represented by the two mode density operator 

P = I^\a,a){a,a\ 

= (l«)H)a®(H(«l). (336) 

where \oi)(^ is a one mode coherent state for mode c = a,b with a = |a| exp(—*0), 
and modes a, b are associated with bosonic annihilation operators a, b. The 
magnitude |a| is fixed. This state gobally but not locally SSR compliant. 

Now 


But 


(a'fb^ = TrJ — a'fb{\a){a\)^(g){\a){a\)i, 

= Tr J ^ {\a) {a\d^)^® (b\a) {a\^^ 


(337) 



(a^a|a) (a|)a® («l)^ 

((a|a^a|a))^(8) ((a|6^S|a))^ 

|a|^ (338) 


Hence we have | P > 0 and | P = (a’^db^by This shows the state 

is SSR compliant entangled. However it fails the Hillery test for SSR neglected 
entanglement which is consistent with being a SSR neglected separable state 
from the i, |43| viewpoint. 
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The second example of such states has an overall density operator which is 
a statistical mixture given by 

P = ^(l^l) (V'll)a «) IV'l) (^ll)b + (V-iDa «) IV'i) (^J)6 

+ ^(|V'-l) (V'-l|)a 0 IV'-l) (V'-Ipfc + (V'-i|)a ® \^-i) (V'-^|)6 

(339) 

where = (|0) +w|l))/-\/2, with w = The Itp^) are superposi¬ 

tions of zero and one boson states and consequently the local particle number 
SSR is violated by each of the sub-system density operators {tpi^\)a and 

Now using 6 IV’^) = (w |0 ))/a/2, = ((0|a;*)/y2 and jwp = 1 

(s.i) = 


But 

(a'^ab%'^ = 


Hence we have | P > 0 and | P < (^^ab^'^. This shows the state 

is SSR compliant entangled. However it fails the Hillery test for entanglement, 
so is consistent with being a SSR neglected separable state 0, m viewpoint. 
It should be noted however that the density operator can also be written as 

P = J(|0)(0|)^®|0)(0|)B + i(|l)(l|)A®|l)(l|)B 

+ ^(I^+>(^+IUb (342) 

where |'I'+)^^ = (|0)^ 11)^ -I- |1)^ |0)^)/-\/2. In this form the terms correspond 
to a statistical mixture of states with 0,1, 2 bosons. The first two terms corre¬ 
spond to separable states, in which the sub-system density operators are SSR 


LO 

LO 

1 1 * 1 


(340) 


® (b^b\^J 


Ujlb^ 


T X I b^b 


iyi|^|2l|^|2 

4 2^ 2' ' 2 ' ' 

LO 

I 

4 


(341) 
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compliant. The final term however is a one boson Bell state which is generally 
regarded as the paradigm of a two mode entangled state. Hence regarding the 
overall state as separable is highly questionable. 
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17 Appendix G - Derivation of Interferometer 
Results 

17.1 General Theory - Two Mode Interferometer 

Introducing the free and interaction evolution operators via 

U = UoU.nt 

Uq = exp{—iHot/h) (343) 

it is straightforward to show that for 

M=^{Vb-a^a) (344) 

we have 

(m) = Tr{MHp) 

= Tr{{MH - {MH)y p) (345) 

giving the mean and variance in terms of the input density operator and inter¬ 
action picture Heisenberg operators 

Mh = — o^h^h) 

bn = U^l b Umt an = U^] a Uint (346) 

where we have used the results bUo = exp{—iujbt) b and Uq^ aUo = exp(—iWat) a- 

The interaction picture Heisenberg operators satisfy 

d d 

ifi-^bn = [bHjVn] ib,—aH = [aH,VH] (347) 

where 

Vh = A{t) exp{—iu}ot) exp{i(j))VfjaH exp{+iu!bat) 

-|-.4(t) exp(-l-iwot) exp{—i(f))a\jbH exp{—iujbat) 

= A{t) exp{i<j)) + A{t) exp{—i(j))a^^bH (348) 

for resonance. 

We then Hnd that the Heisenberg picture operators satisfy coupled linear 
equations 

a- _ a_ 

ih—bn = A{t) exp{+i(j))aH ~ -^(0 &xp{—i<p) bn (349) 

which after replacing the time t by the area variable s then involve time inde¬ 
pendent coefficients 

^ d ^ 

^-^bH{s) = exp{+i(j))aH{s) i — anis) = exp{-i<j)) bnis) 
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(350) 



The equations can then be solved via Laplace transforms giving 

bnis, (f)) = cos sb — i exp(z^) sin s a anis, (j)) = exp(—sin sb + cos s a 

(351) 

where now 2s is the area for the classical pulse. 

Hence we have in general 

MH{2s,(j)) = i(6|^(s,(())6//(s,(()) - a|^(s,(/))a//(s») 

= sin2s (sin^S'a; + cos(/>5y) + cos2s^^ (352) 


The versatility of the measurement follows from the range of possible choices of 
the pulse area 2s and the phase (j). 

Writing 2s = 9 we can then substitute into EqJ ^5 ) to ob tain result s for 
A'l'j and These are set out in SubSection 


7.3 


in Eqs. (218) and (219) 


in terms of the mean values of the spin operators and the matrix elements of 
the covariance matrix for the spin operators, all for the quantum statep. 


17.2 Beam Splitter and Phase Changer 

For the beam splitter we have 2s = 7r/2 and (j) (variable) so that 

(/)) = sin^Sa; + cos^S'y (353) 

whilst for the phase changer we have 2s = tt and (j) (arbitrary) so that 

Mh{tt,(1>) =-S^ (354) 


17.3 Other Measureables 

We can also consider other choices for the measureable, which then enable us 
to determine other moments of the spin operators. A case of particular interest 
is the square of the population difference 

M 2 = ^-(6^6 —a^a)^ (355) 

It is then straightforward to show for the beam splitter case with 2s = 7r/2 and 
(j) (variable) 

M2hC^, 4>) = {sin(j)Sx +cos (j)S^ 

= sin^ (p + cos^ (j) (S'y)^ + sin p cos p {SxSy + SySx) 

(356) 


Hence 

= sin^p (|(5a;)^^ +cos^ p ({Syf'^ + sinp cosp ((SxSy + SySx)'j (357) 
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showing that the mean for the new observable M 2 is a sinusoidal function of the 
BS interferometer variable </> with coefficients that depend on the means of 3“^, 
Si and Sjy + SyS^. 

17.4 General Theory - Multi-Mode Interferometer 

The derivation follows the same steps as in SubSection |17.1[ However here we 
have the results bi Uq = exp{—i{u}b + 0Ji)t) hi and ^ a-i Uq = exp{—i{(jja + 
uii)t)ai. The factors involving exp[—iuji)t cancel out in the derivation of the 
Heisenberg equations, which here are 


o ^ ^ o ^ ^ 

= A{t) exp{+i4>)aiH ih—a-iH = A{t) exp{-i(j)) b.^H (358) 

and the solutions are 

biH^s^/p) = cos s bi—i exp{i(j)) sinsdi aiHis,4>) = —jexp(—sin s 5^+cos sdi 

(359) 

where now 2s is the area for the classical pulse. 

Hence we have in general 




= sin2s (sin^S'j; + cosc/iS'y) + cos 2 s.S 2 


(360) 


This leads to the same formal results (2351 and ( |236 ) for the mean and variance. 
The versatility of the measurement again follows from the range of possible 
choices of the pulse area 2s and the phase p. 
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18 Appendix H - Limits on Interferometry Tests 


The tests for entanglement in a particular quantum state are given in terms of 
the mean value and variance for certain physical quantities. Interferometers are 
used to enable these means and variances to be determined from measurements 
on another physical quantity when either the state being tested is acted upon 
by the interferometer or it is being unaffected. Quantum theory enables us to 
predict two things. Firstly, for any physical quantity M we can predict the 
possible values that measurements could result in. Results from a succession of 
measurements would confirm what these values are. Secondly, for any quantum 
state, we can predict the probability that measurement leads to a specific value. 
A single measurement only yields one of the possible values, so independent 
repetitions of such measurements is needed to confirm what the probabilities 
for measuring particular values are - ideally an infinite number of repeated 
measurements would be required. If this was possible, the computed mean 

and variance of the measurements for the physical quantity M would 

confirm the quantum theory predictions for any quantum state. A finite but 
large number of independent measurements - each based on the same probabiity 
distribution for the possible results, would enable us to estimate the actual mean 
and variance of the measured values from the sample measurements. These 
estimates would not be precisely accurate. The question is - how big would 
the sample of repeated measurements need to be for the purpose of using the 
estimated mean and variance in the tests for entanglement ? 

Statistical theory in the form of the central limit theorem |44) can be applied 
here. This tells us if the number R of repeated measurements is large, then the 
mean of the sample measurements approaches the true mean and the variance 
in the sample estimation of the mean is given by the true variance divided by 
R 


(M) 

\ ! sample 




2 

a/m 


sample 


—)• 



(361) 


We can use our theoretical estimate of the variance to get an idea 

of how large the sample of measurements must be in order that the standard 
deviation of the sample estimate for the mean is small enough that the mean 
can confidently be stated to exceed or be less than the quantity on the other 
side of the inequality in the entanglement test. 
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19 Appendix I - Relative Phase State 


The relative phase eigenstate (see i, m) for an N boson two mode system 
has provided an important example of different outcomes for the simple spin 
squeezing and Hillery spin squeezing tests, so here its properties are set out in 
more detail. The results for interferometric measurements on the relative phase 
state are also presented. 

The relative phase state is a globally compliant entangled state of the sub¬ 
systems a and b and is defined by 

N/2 

= ^ eM^kO^)\N/2-k)JN/2 + k), (362) 

k=-N/2 


where 9p = p{2TTj {N + 1)), p = —N/2, —N/2 + 1, -I-7V/2 is a quasi-continuum 
of TV -I- 1 equispaced phase eigenvalues, and |iV/2 — k)^, \N/2 + /c)^ are Fock 
states for sub-systems a and b. The Hermitian relative phase operator 0jv for 
N boson states is then defined as 

= ( 363 ) 

p 


and \N,9p) is an eigenvector with eigenvalue 9p . 

Since these states are entangled with maximum mode entropy, are spin 
squeezed and are fragmented BEC (two modes have macroscopic occupancy) 
it is of some interest to examine their interferometric properties for the simple 
beam splitter interfometer. As shown in [5] the relative phase state has the 
following mean values for the spin operators when /> = \N, 9p) {N, 9p\ 


= ^cos9p (^Sy'^ =-^sm9p (s^'^ = 0 (364) 

so that for the measurable 

^ sm{cj) - 9p) (365) 

We thus have a large amplitude - proportional to V - sinusoidal dependence 
for the mean value of the measureable on the interferometer phase detuning 
((/> — 9p), and which goes to zero when cf = 9p. Since we never have both 

and equal to zero the simple correlation test confirms that the relative 

phase eigenstate is entangled. 

As mentioned above, the relative phase state is highly spin squeezed. To 
describe this it is convenient to introduce rotated spin operators J^, Jy and 
given by (see Ref [5], Eqn. 179) 



Jy = sin 9p Sx + cos 9p Sy 
Jz = — cos 9p Sx + sin 9p Sy 


(366) 
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The new spin operators are Schwinger spin operators for new modes c, d where 


d = — exp(ii0p) — d'j /v^ b = — exp{—^idp) + d'j 


(367) 


and the relative phase state also an entangled state for new modes. This can 
be shown by substituting for the |A^/2 — and \N/2 + k)f^ in terms of Fock 
states for the new modes c, d. 

These new angular momentum operators are principal spin operators for 
which the covariance matrix is diagonal. For the mean values 


Jx) =0 


Jy) =0 


P) = — 


Ntt 


(368) 


In terms of spin operators discussed above (see Eqs. (2381 and (239)) we have 


Jx = Sz, Jy = Sfpf + 6p) and + dp) so the variances for Jy and 

Jz can be measured using the simple BS interferometer, and the mean for is 
also measureable by simply measuring the mean population difference without 
subjecting the relative phase eigenstate to the BS interaction. 

Inverting these expressions and substituting gives the measureable in terms 
of the new spin operators 


Mh 


cos{(j) — Op) Jy — sin(()) — Op) Jz 


(369) 


Hence we find for the variance of the measureable 

= cos'^ P - 0 p) C {Jy Jy)+siP{ct>-6 p)C{JzJ.) 

—2 siii{4> — Op) cos{4> — Op) C{Jy, Jz) (370) 

As Jx,Jy and Jz are principal spin operators C{Jy, Jz) = 0 and substituting 
for the variances C{Jy, Jy) = 1/4 + 1/8 InN and C{Jz, Jz) = (1/6 — 7r^/64) 

(see i) we get for the variance of the measureable for an input relative phase 
eigenstate 

= COS^(0-^p)(i + ilniV)+ sin^(0-^p)(i - 

“ ) + ( 371 ) 

for (j) «. Op. The other variance is C{ Jxt Jx) = (1/12) iV^. The variance for 
the measurable depends sinusoidally on 2{(j) — Op). Thus the quantum noise in 
the measureable also goes to essentially zero at 4> = Op, when the mean value 

also goes to zero. The width for this low noise window scales as 
1/A - which corresponds to the Heisenberg limit. At the zero of the mean 
value, the relative fluctuation varies as 1/A as in the Heisenberg limit. Since 
for (j) = Op we have Mh = Jy = Sf{^ + Op) and (^AM^^ = (1 + |lnA) 
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whilst ~ Thus the spin squeezing test in Eq.(|24l|) is 

satished, confirming again that the relative phase eigenstate is an ent 
state of modes a and b. 


angled 
and 


In regard to particle entanglement [45], [46] with p = \N,9p) {N,0p 
with Ua = {N/2 — k), Ub = {N/2 + k), the quantities in Eqs. ( 132 ) and ( 133 ) 
of paper 1 are given by 


P 


{riant,) _ 


1 


Pr, 


+ 1 
1 

fvTT 


|fV/2 - fc), {N/2 - ® \N/2 + k), {N/2 + k\, (372) 

(373) 


and since is a separable state, it follows that Ep{ p) = 0. Thus the 

measure of particle entanglement is zero for what is clearly an entangled state. 
Hence the particle entanglement measure has not detected entanglement in this 
example. 

The relative phase state is therefore a promising candidate for use as an 
input state in two mode interferometry. More elaborate interferometers where 
the interferometric variable is associated with other systems whose parameters 
are to be measured might be developed. The main issue would be whether 
such a relative phase state could be prepared. This is an issue being dealt with 
elsewhere [47]. 
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